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ogy, with particular emphasis on applications to the problem of post-inflation reheat- 
ing. The Schwinger-Keldysh closed-time-path (CTP) formalism is utilized along with 
the two-particle-irreducible (2PI) effective action in order to obtain coupled, nonper- 
turbative equations for the mean field and variance in a general curved background 
spacetime, both as a closed system in the case of a self-interacting inflaton field, and 
as an open system in the case of coarse-grained dynamics of the inflaton field interact- 
ing with fermions. For a model consisting of a quartically self-interacting O(A^) field 
theory (with unbroken symmetry) in spatially fiat FRW spacetime, the dynamics of 
the mean field is studied numerically, at leading order in the large- expansion, with 
initial conditions appropriate to the end state of slow roll in chaotic inflation scenarios. 
The time evolution of the scale factor is determined self-consistently using the semi- 
classical Einstein equation. It is found that cosmic expansion can dramatically affect 



the efficiency of parametric resonance-induced particle production. The production of 
fermions due to the oscillating inflaton mean field is studied for the case of a scalar 
inflaton coupled to a fermion field via a Yukawa coupling /. The dissipation and noise 
kernels appearing at 0(/^) in the one-loop CTP effective action are shown to sat- 
isfy a zero-temperature fluctuation-dissipation relation (FDR). The normal-threshold 
0(/^) parts of the one-loop CTP effective action are also shown to satisfy a FDR. 
The effective stochastic equation obeyed by the inflaton zero mode at 0(/^) contains 
multiplicative noise. It is shown that stochasticity becomes important to the dynamics 
of the inflaton zero mode before the end of reheating. The thermalization problem 
is discussed, and a strategy is presented for obtaining time-local equations for equal- 
time correlation functions which goes beyond the Hartree-Fock approximation. For the 
field theory, the correlation entropy associated with a particular coarse graining 
scheme consisting of slaving the three-point function to the mean field and two-point 
function is computed, and found not to be conserved. 
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CHAPTER 1 



Introduction 



1.1 Background 

The inflationary Universe [1-12] has for over a decade been the new paradigm for 
addressing many basic issues in cosmology such as the spatial flatness-oldness prob- 
lem, the large-scale homogeneity (horizon) problem, and the small-scale inhomogene- 
ity (structure formation) problem. The linkage between observations, especially those 
from the recent Cosmic Background Explorer (COBE) data, and theory, based on 
grand unified theories (GUT's) and Friedmann-Robertson- Walker- (FRW-)de Sitter 
models, has been pursued in earnest, but most theoretical discussions to date are 
largely phenomenological and somewhat utilitarian in nature [13,14]. This lack of 
rigor and precision is understandable for at least two reasons: the precise physical 
conditions between the Planck and GUT scales (when the most cosmologically signif- 
icant inflationary evolutions are believed to have taken place) have not been clearly 
understood, and the theoretical framework for the treatment of processes affecting the 
inception and completion of inflation were not well developed. As stressed earlier by 
[15, 16], the important physical processes which can determine whether inflation can 
occur, sustain, and flnish with the necessary features are affected by at least three 
aspects: the geometry, topology, and dynamics of the spacetime [17], the quantum 
field theory aspects pertaining to the analysis of infrared behavior, and the statistical 
mechanical aspects pertaining to nonequilibrium processes. These quantum and sta- 
tistical processes include phase transition, particle creation, entropy generation, fluc- 
tuation or stochastic dynamics, and structure formation [18, 19]. Most of these invoke 
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the quantum field and statistical mechanical aspects, and for processes occurring at 
the Planck scale (which are instrumental in starting certain models of inflation, such as 
proposed in [5,20,21]), also the geometry and topology of spacetime. Two important 
problems involving field theory in curved spacetime [17], namely, the back reaction 
of cosmological particle creation [22-26] on the structure and dynamics of spacetimes 
[20, 25, 27-34], and the effects of geometry and topology of spacetime on cosmological 
phase transitions [16, 21, 35-38], were investigated systematically and comprehensively 
in the 1970s and 1980s. The statistical mechanical aspect has not been considered with 
equal mastery. 

The statistical mechanical aspect enters into all three stages of inflationary cos- 
mology: (i) At the inception: What conditions would be most conducive to starting 
inflation? Do there exist metastable states for the Higgs boson field which can gener- 
ate inflation [39]? Can thermal or quantum fluctuations assist the inflaton in hopping 
or tunneling out of the potential barrier in the spinodal or nucleation pictures? Most 
depictions so far have been based on the finite temperature effective potential, which 
assumes an unrealistic equilibrium condition and a constant background field. How- 
ever, when asking such questions in critical dynamics one should be using a Langevin 
or Fokker-Planck equation (a generalized time-dependent Landau-Ginzberg equation 
[40]) incorporating dynamic dissipation and intrinsic noise consistently, (ii) During 
inflation, the dynamics of the inflaton field can be more easily understood in terms 
of a Kadanoff-Migdal exponential scaling transform [41]. The reason why the infla- 
ton evolves as a classical stochastic field [42-44] at late times involves the process 
of decoherence, caused by noise and fluctuations from environmental fields [45]; this 
necessitates statistical mechanical considerations. The evolution of the classical den- 
sity contrast (containing the seedings of structures) from quantum fluctuations of the 
inflaton also requires both quantum and stochastic field theory considerations [46-55] . 
(iii) In the reheating epoch, particle creation induces dissipation of the inflaton field. 
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and the interaction of quantum fields is the source for reheating the Universe. This 
last epoch is the focus of this dissertation, as we shall detail below. 

The construction of a viable theoretical framework for treating quantum statistical 
processes in the early Universe has been underway for the past decade (for a review, 
see [56]). This framework has now been successfully established, and its application to 
the problems mentioned above has just begun. The cornerstones are the Schwinger- 
Keldysh closed-time-path (CTP) [57-69] effective action and the Feynman- Vernon 
influence functional [70-76] formalisms. They are useful for treating particle creation 
back reaction [67, 69], fluctuation or noise, and dissipation or entropy problems [53, 54, 
76] . Other essential ingredients include the Wigner function [77, 78] , the n-particle- 
irreducible (nPI) effective action [37,68,79,80], and the correlation hierarchy [81,82] 
for treating kinetic theory processes [68, 83] and phase transition problems [54, 84]. In 
this dissertation we apply these techniques to the problems of inflaton damping due 
to back reaction from parametric particle creation (Chapter 3) and dissipation due to 
particle creation (Chapter 4), which are relevant in in the third epoch depicted above. 
In parallel, these newly developed methods in statistical field theory are now being 
applied to derive the classical stochastic dynamics of the inflaton (in the second epoch) 
[45] , and the statistical field theory of spinodal decomposition (in the first epoch) [40] . 

1.2 Issues 

Most all inflationary cosmologies share the feature of a period of cosmic expansion 
driven by a nearly constant vacuum energy density p (a "vacuum-dominated" era 
with effective equation of state, p = —p): In a Priedmann- Robertson- Walker (FRW) 
spacetime, the scale factor expands exponentially in cosmic time, resulting in extreme 
redshifting of the energy density of all other forms of matter and fields. As long as the 
interaction time scale of any physical process involving given fields is longer than the 
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cosmic expansion time , the fields will remain in disequilibrium. This condition 
can prevail in all three stages of inflation, and one should use a fully nonequilibrium, 
nonperturbative treatment of the dynamics of the inflation field. The physics of the 
reheating epoch is important because it directly determines several important cos- 
mological parameters which are relevant to later evolution of the Universe, and in 
principle verifiable by observational data. For example, the reheating temperature is 
a vital link between the inflationary Universe scenario and GUT scale baryogenesis 
[86] , and may provide a mechanism to explain the origin of dark matter [87, 88] . 

It is generally believed that at the end of inflation, the state of the inflaton field can 
be approximately described by a condensate of zero-momentum particles undergoing 
coherent quasioscillations about the true minimum of the effective potential [10, 11, 89]. 
The reheating problem involves describing the processes by which the many light fields 
coupled to the inflaton become populated with quanta, and eventually thermalize. It 
is commonly believed that if the flelds interact sufficiently rapidly and strongly, the 
Universe thermalizes and turns into the radiation-dominated condition described by 
the standard Priedmann solution, but this has not been proven satisfactorily. 

There has been a great deal of work over the past 15 years on the reheating problem, 
and in attempting to understand reheating, a wealth of interesting physics has been 
revealed (see, e.g., [90]). To date, the work on particle production during reheating 
largely follows two distinct approaches, each pursued in two stages. 

In the first stage of work on the reheating problem (group lA, [91-94]), time- 
dependent perturbation theory was used to compute the rate of particle production 
into light fields (usually fermions) coupled to the inflaton. Particle production rates 
were computed in flat space assuming an eternally sinusoidally oscillating inflaton 
field. The inflaton evolution in FRW spacetime was modeled with a phenomenolog- 
ical c-number equation involving the Hubble parameter H and the classical inflaton 
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amplitude ^, 

4> + m^(t)+{r^ + 3H)^ = 0, (1.1) 

where T^, given by the imaginary part of the self-energy of (/), is the total perturbative 
decay rate, and (p = d(j)/dt. Bose enhancement of particle production into the spatial 
Fourier modes of the inflaton fluctuation field tp (and light Bose fields coupled to the 
inflaton) was not taken into account. 

In the second stage of this first approach to the reheating problem (group IB, 
[87,88,95,96]) Eq. (1.1) was still utilized to model the mean-field dynamics, but 
with computed beyond first-order in perturbation theory. In the work of Shtanov, 
Traschen, and Brandenberger [95] and Kofman, Linde, and Starobinsky (KLS) [87], 

was computed for a real self-interacting scalar inflaton field 4> which was both 
Yukawa-coupled to a spinor field ip, and bi-quadratically coupled to a scalar field x 
(KLS studied both the cf) ^ —(f) symmetry-breaking and unbroken symmetry cases). 
Prom the one-loop equations for the quantum modes of the x, tp, and fields (in which 
the mean field appears quadratically as an effective mass), approximate expressions 
for the growth rate of occupation numbers were derived, assuming a quasi-oscillatory 
mean field 4>- For bosonic decay-product fields, it was found that first-order time- 
dependent perturbation theory drastically underestimates the particle production rate 
for modes which are in an instability band for parametric resonance.^ Parametric am- 
plification of quantum fluctuations^ in Bose decay-product fields can result in rapid 
out-of-equilibrium transfer of energy from the inflaton mean field to the (spatially) 
inhomogeneous inflaton modes and light Bose fields coupled to the inflaton. This 
phenomenon was called preheating by KLS. It has been suggested that exponential 

^Thc first study to point out tiiat parametric resonance effects can dramatically effect particle 
production in an out-of-equilibrium phase transition was [89] . 

^Parametric amplification of quantum fluctuations refers to the increase in expectation values of 
occupation numbers for parametric oscillators, due to a time-dependent perturbing frequency. 
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growth of quantum fluctuations can in some cases lead to out-of-equilibrium (nonther- 
mal) symmetry restoration in the "new" inflation models with a spontaneously broken 
symmetry [97,98]. (See, however, the work of Boyanovsky et al, which reached a dif- 
ferent conclusion on the possibility of nonequilibrium symmetry restoration [99].) This 
may have interesting implications for baryogcnesis, defect formation, and generation 
of primordial density perturbations [87,90,98]. 

In both stages of this first approach, the back reaction of the variance of the inflaton 
on the mean-field dynamics, and of the variance on the quantum mode functions, were 
not treated self-consistently. The effect of spacetime dynamics was either excluded 
entirely, or not included self-consistently using the semiclassical Einstein equation. Due 
to the potentially large initial inflaton amplitude at the onset of reheating, particularly 
in the case of chaotic inflation [11], the effect of cosmic expansion on quantum particle 
production needs to be included. Since the mean field and variance (mean-squared 
fluctuations) are coupled, the back reaction of particle production on the mean-field 
dynamics must be accounted for in a self-consistent manner. 

In the decade before the advent of inflationary cosmology, there was active research 
on quantum processes in curved spacetimes. An important class of problems is vacuum 
particle creation [22-26] and its effect on the dynamics and structure of the early 
Universe [25, 27-34] at the Planck time. The effect of spacetime dynamics and the 
importance of parametric amplification on cosmological particle creation were realized 
very early [22,24,26]. Most of the effort in the latter part of the 1970s was focused 
on obtaining both a regularized energy-momentum tensor and a viable formalism for 
the treatment of back reaction effects. The wisdom gained from work in that period 
before the infiationary cosmology program was initiated is particularly relevant to the 
reheating problem. Simply put, for obtaining a finite energy- momentum tensor for a 
quantum field in a cosmological spacetime, the adiabatic [26, 100-102] and dimensional 
[103] regularization methods are the most useful. For studying the back reaction of 
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particle creation, the Schwinger-Keldysh (CTP, "in-in") effective action formalism 
[57,59,63-67,69] is more appropriate than the usual Schwinger-DeWitt ("in-out") 
method [104, 105]. 

The second approach to the post-inflationary reheating problem is built upon the 
body of earlier work on cosmological particle creation. Following the application of 
closed-time-path techniques to nonequilibrium relativistic field theory problems [67, 
68] , several authors (which we call group 2A) derived perturbative mean-field equations 
for a scalar inflaton with cubic [69] and quartic [106] self-couplings, as well as for 
a scalar inflaton Yukawa-coupled to fermions [107]. The closed-time-path method 
yields a real and causal mean-field equation with back reaction from quantum particle 
creation taken into account. For the case of Bose particle production, perturbation 
theory in the coupling constant is known to break down for sufficiently large occupation 
numbers, which occurs on the time scale ri for parametric resonance effects to become 
important [108, 109]. It is, therefore, necessary to employ nonperturbative techniques 
in order to study reheating in most inflationary models. 

The second stage of work in this second approach to the reheating problem used 
the closed-time-path method to derive self-consistent mean-field equations for an in- 
flaton coupled to lighter quantum fields (group 2B, [110-116]). In the first of these 
studies [110-113], the coupled one- loop mean-field and mode-function equations were 
solved numerically in Minkowski space, implicitly carrying out an ad hoc nonpertur- 
bative resummation in h. In the one- loop equations, the variances for the inflaton 
((^^) and light Bose fields (x^) do not back-react on the mode functions directly. How- 
ever, mean-field equations were derived for an 0(Ar)-invariant linear a model (with a 
self-interaction) at leading order in the large- AT approximation by Boyanovsky et 
al. [108]. In this approximation, the variance docs back-react on the quantum mode 
functions. At leading order in the 1/iV expansion, the unbroken symmetry dynamical 
equations for the quartic 0{N) model are formally similar to the dynamical equations 
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for a single A<I>'^ field theory in the time-dependent Hartree-Fock approximation [79]. 
The nonequilibrium dynamics of the quartically self-interacting 0(A'^) field theory in 
Minkowski space has been numerically studied at leading order in the 1/iV expan- 
sion in both the unbroken symmetry [99, 108, 117] and symmetry-broken [99, 108, 118] 
cases. Some analytic work has been done on the self-consistent Hartree-Fock mean- 
field equations for a quartic scalar field in Minkowski space [109]. In addition, the 
Hartree-Fock equations for a field in the slow-roll regime have been studied nu- 
merically in Minkowski space [119] and in FRW spacetime [120]. However, the effect of 
spacetime dynamics on reheating in the O(A^) field theory has not (to our knowledge) 
been studied using the coupled, self-consistent semiclassical Einstein equation and 
matter-field dynamical equations, though some simple analytic work has been done 
on curvature effects in reheating [96,114]. The semiclassical equations for one-loop 
reheating in FRW spacetime were derived in [115]. The (t>^X^ theory has been studied 
in FRW spacetime by [53, 116, 121]. In addition, numerical work has been done on 
symmetry-breaking phase transitions in both a scalar field in de Sitter spacetime 
[122], and an 0{N) theory in FRW spacetime [123, 124]. 

1.3 Organization 

This dissertation is organized as follows. In Chapter 2, which describes work published 
in Ref. [80], we construct the two-particle- irreducible (2PI), closed-time-path (CTP) 
effective action for the 0{N) field theory in a general curved spacetime. Prom this 
we derive a set of coupled equations for the mean field and the variance, which are 
useful for studying the nonperturbative, nonequilibrium dynamics of a quantum field 
when full back reactions of the quantum field on the curved spacetime, as well as the 
fluctuations on the mean field, are required. Renormalization of the effective action 
at leading order in the 1/N expansion is then discussed. 
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In Chapter 3, which describes work pubhshed in Ref [125], we study the non- 
perturbative, nonequilibrium dynamics of a quantum field in the preheating phase 
of inflationary cosmology, including full back reactions of the quantum field on the 
curved spacetime, as well as the fluctuations on the mean field. We use the 0{N) 
field theory with unbroken symmetry in a spatially flat FRW spacetime to study the 
dynamics of the inflaton in the post-inflation, preheating stage. Oscillations of the in- 
flaton's zero mode induce parametric amplification of quantum fluctuations, resulting 
in a rapid transfer of energy to the inhomogeneous modes of the inflaton field. The 
large-amplitude oscillations of the mean field, as well as stimulated emission effects 
require a nonperturbative formulation of the quantum dynamics, while the nonequi- 
librium evolution requires a statistical field theory treatment. We adopt the coupled 
nonperturbative equations for the mean field and variance derived in Chapter 2 while 
specialized to a dynamical FRW background, up to leading order in the 1/A^ expan- 
sion. Adiabatic regularization is employed. The renormalized dynamical equations are 
evolved numerically from initial data which are generic to the end state of slow roll 
in many inflationary cosmological scenarios. We find that for sufficiently large initial 
mean- field amplitudes > Mp/300 (where Mp is the Planck mass) in this model, the 
parametric resonance effect alone (in a collisionless approximation) is not an efficient 
mechanism of energy transfer from the mean field to the inhomogeneous modes of the 
quantum field. For small initial mean-field amplitude, damping of the mean field due to 
particle creation is seen to occur, and in this case can be adequately described by prior 
analytic studies with approximations based on field theory in Minkowski spacetime. 

In Chapter 4, which describes work to be published [126], we present a detailed 
and systematic analysis of the coarse-grained, nonequilibrium dynamics of a scalar 
inflaton field coupled to a fermion field in the late stages (dominated by fcrmion particle 
production) of the reheating period of inflationary cosmology with unbroken symmetry. 
We derive coupled nonperturbative equations for the inflaton mean field and variance 
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at two loops in a general curved spacetime, and show that the equations of motion 
are real and causal, and that the gap equation for the two-point function is dissipative 
due to fermion particle production. We then specialize to the case of Minkowski 
space and small-amplitude inflaton oscillations, and derive the perturbative one-loop 
dissipation and noise kernels to fourth order in the Yukawa coupling constant; the 
normal-threshold dissipation and noise kernels are shown to satisfy a zero-temperature 
fluctuation-dissipation relation. We derive a Langevin equation for the dynamics of 
the inflaton zero mode. We then show that the variance of the inflaton zero moe can 
be non-negligible during reheating, which is the primary physical result of the chapter. 

In Chapter 5, which describes work to be published [127], we set the stage for a 
study of the thermalization process in reheating by investigating how entropy can be 
defined for an interacting quantum field. We discuss various definitions of entropy 
but focus our attention on the correlation entropy. We discuss how an effectively 
open system arises when hierarchy of correlation functions is truncated and one of 
the higher correlation functions is slaved to the lower correlation functions. We show 
how the dynamics of a nonperturbative truncation of the Schwinger-Dyson equations 
can be reduced to coupled equations for the equal-time correlation functions. We 
then compute the correlation entropy for the case of the truncated at third order 
in the correlation hierarchy, where the three-point function is slaved to the mean 
field and two-point function, for the case of a (spatially) translationally invariant 
Gaussian density matrix. We then discuss the possible benefits of this approach to the 
thermalization problem. 

1.4 Notation 

Throughout this dissertation we use units in which c = = 1. Planck's constant h 
is shown explicitly (i.e., not set equal to 1) except in Chapter 5. In relativistic units 
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where c = 1, Newton's constant is G = HMp^'^, where Mp is the Planck mass. We 
work with a four-dimensional spacetime manifold, and follow the sign conventions^ of 
Birrell and Davies [17] for the metric tensor Qf^i,, the Riemann curvature tensor Rf^uap, 
and the Einstein tensor G^^. We use greek letters to denote spacetime indices. The 
beginning latin letters a, b, c, d, e, f are used as time branch indices (see Sec. 2.2), and 
in Chapters 2 and 3, the middle latin letters i,j,k,l,m,n are used as indices in the 
0{N) space (see Sec. 2.5). In Section 5.3 the middle latin letters are used as indices to 
indicate spatial coordinate. The Einstein summation convention over repeated indices 
is employed. Covariant differentiation is denoted with a nabla or a semicolon. 



^In the classification scheme of Misner, Thome and Wheeler [128], the sign convention of Birrell 

and Davies [17] is classified as (+, +, +). 
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CHAPTER 2 



O(A^) Quantum fields in curved spacetime 
2.1 Introduction 

One major direction of research on quantum field theory in curved spacetime [17, 129, 
130] since the 1980s has been the apphcation of interacting quantum fields to the 
consideration of symmetry breaking and phase transitions in the early Universe, from 
the Planck to the grand unified energy scales [123, 131-138]. In a series of work, Hu, 
O'Connor, Shen, Sinha, and Stylianopoulos [16,21,35-37,39,139-141] systematically 
investigated the effect of spacetime curvature, dynamics, and finite temperature in 
causing a symmetry restoration of interacting quantum fields in curved spacetime. In 
general one wants to see how quantum fluctuations around a mean field (j) change 
as a function of these parameters. For this purpose, the two-particle-irreducible (2PI) 
effective action was constructed for an A'^-component scalar 0{N) model with quartic 
interaction [35, 37, 137]. Hu and O'Connor [37] found that the spectrum of the small- 
fluctuation operator contains interesting information concerning how infrared behavior 
of the system depends on the geometry and topology. The equation for cf) containing 
contributions from the variance of the fluctuation field [ip^) depicts how the mean field 
evolves in time. This program explored two of the three essential elements of an inves- 
tigation of a phase transition [16] , the geometry and topology and the field theory and 
infrared behavior aspects, but not the nonequilibrium statistical-mechanical aspect. 

For this and other reasons, Calzetta and Hu [67] started exploring the closed- 
time-path (CTP) or Schwinger-Keldysh formalism [57-59,63], which is formulated 
with an "in-in" boundary condition. Because the CTP effective action produces a 
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real and causal equation of motion [65,66], it is well suited for particle production 
back-reaction problems [31,32,69]. Use of the CTP formalism in conjunction with 
the 2PI effective action [79] and the Wigner function [77] enabled Calzetta and Hu to 
construct a quantum kinetic field theory (in flat spacetime), deriving the Boltzmann 
field equation from first principles [68] . The necessary ingredients were then in place for 
an analysis of nonequilibrium phase transitions [84]. In recent years these tools (CTP, 
2PI) have indeed been applied to the problems of heavy-ion collisions, pair production 
in strong electric fields [117], disoriented chiral condensates [142, 143], and reheating 
in inflationary cosmology [144]. However, none of these recent works has included 
curved spacetime effects in a self-consistent manner, where the spacetime governs the 
evolution of a quantum field and is, in turn, governed by the quantum field dynamics. 
This is especially important for Planck scale processes involving quantum fluctuations 
with back reaction, such as particle creation [76], galaxy formation [54], preheating, 
and thermalization in chaotic inflation [122, 124]. 

In this Chapter we return to the problems begun by Calzetta, Hu and O'Connor 
a decade ago. We wish to derive the coupled equations for the evolution of the mean 
field and the variance for the 0{N) model in curved spacetime, which should provide a 
solid and versatile platform for studies of phase transitions in the early Universe. The 
first order of business is to construct the CTP-2PI effective action in a general curved 
spacetime. The evolution equations are derived from it. We must also deal with the 
divergences arising in it. From the vantage point of the correlation hierarchy (and the 
associated master effective action) as applied to a nonequilibrium quantum field [82], 
there is a priori no reason why one should stop at the 2PI effective action. Indeed, 
the 2PI effective action corresponds to a further approximation from the two-loop 
truncation of the master effective action constructed from the full Schwinger-Dyson 
hierarchy [81, 82]. For problems where the mean field and the two-point function give 
an adequate description (which is not the case near the critical point, where higher- 
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order correlation functions become important [145]), the CTP-2PI effective action is 
sufficient. In particular, the 2PI effective action contains the commonly used leading- 
order laige-N, time-dependent Hartree-Fock, and one-loop approximations. 

The 0{N) model has been usefully applied to a great variety of problems in field 
theory and statistical mechanics [146]. At leading order in the large- expansion, the 
0{N) field theory yields nonperturbative,^ local evolution equations for the mean field 
and the modes of the fluctuation field, which are valid in the regime of strong mean field 
[79, 117]. This approximation has recently been applied to problems of nonequilibrium 
phase transitions [99, 108, 118]. In the "preheating" problem studied in Chapter 3, we 
shall see that it is particularly useful for describing the nonperturbative dynamics 
of the infiaton field in chaotic infiation scenarios [5], where the infiaton mean-field 
amplitude can be on the order of the Planck mass at the end of the slow roll period 
[11, 147]. The 1/N expansion has many attractive features, as it is known to preserve 
the Ward identities for the 0(iV) theory [148] and to yield a covariantly conserved 
energy-momentum tensor [149]. Furthermore, in the limit of large N, the quantum 
effective action for the matter fields can be interpreted as a leading-order term in the 
expansion of the full (matter plus gravity) quantum effective action [149]. 

Mazzitelli and Paz [150] have studied the and 0{N) field theories in a general 
curved spacetime in the Gaussian and large-AT approximations, respectively. Their 
approach differs from ours in that it is based on a Gaussian factorization which does 
not permit systematic improvement either in the loop expansion or in the 1/N ap- 
proximation. In contrast, our approach is based on a closed-time-path formulation of 
the correlation dynamics. The evolution equation we obtain for the two-point function 
contains a two-loop dissipative contribution (due to multiparticle production) which is 

^By this, we mean that the solution to the coupled equations for the mean field and inhomogeneous 
modes represents a nonperturbative resummation of an infinite subclass of diagrams in the ordinary 
IPI effective action, which is a functional of the mean field only. 
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not present in the large- approximation. At leading order in the large- A'^ approxima- 
tion, our results agree with theirs, so that their renormalization counterterms can be 
directly applied to the leading-order-large- AT limit of the mean field and gap equations 
derived here. 

This chapter is organized as follows. In Sees. 2.2 and 2.3 we present self-contained 
summaries of the two essential theoretical methodologies employed in this study, the 
closed-time-path formalism and the two-particle-irreducible effective action. The adap- 
tation of these tools to the quantum dynamics of a field theory in curved spacetime 
is presented in Sec. 2.4. The 0{N) scalar field theory is treated in Sec. 2.5, where we 
study the two-loop truncation of the 2PI effective action. 

2.2 Schwinger-Keldysh formalism 

The Schwinger-Keldysh or "closed-time-path" (CTP) formalism is a powerful method 
for deriving real and causal evolution equations for expectation values of operators 
for quantum fields in disequilibrium. A quantum field may be defined to be out 
of equilibrium whenever its density matrix p and Hamiltonian H fail to commute, 
i.e., [H,p\ ^ 0. Such conditions can occur, for example, in a field theory quantized 
on a dynamical background spacetime, and also in an interacting field theory with 
nonequilibrium initial conditions. Although in Chapters 2 and 3 we shall be concerned 
with closed-system, unitary dynamics of a single self-interacting quantum field, the 
methods discussed here are also well suited to studying the dynamics of an open 
quantum system, as shown in Chapter 4 below. Excellent reviews of the Schwinger- 
Keldysh method are Chou et al. [63] as applied to nonequilibrium quantum field theory 
and Calzetta and Hu [67] as applied to the back reaction problem in semiclassical 
gravity. 

Let us briefiy review the Schwinger-Keldysh method as applied to the effective 
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mean-field dynamics of a self-interacting scalar field theory in Minkowski space. The 
classical action for a scalar theory in Minkowski space is 



(2.1) 



where A is a coupling constant with dimensions of 1/h (inverse mass times inverse 
length), m is the "mass" with dimensions of inverse length, and r]'^'^ = diag(-|-, —,—,—) 
is the Minkowski space metric tensor. The Euler-Lagrange equations are obtained by 
variation of the action, where it is understood that the variations of must satisfy 
boundary conditions in order that surface terms can be discarded. 

Let us denote the Heisenberg field operator for the canonically quantized theory 
with classical action (2.1) by ^h{x)- By "effective mean-field dynamics" we mean that 
we seek a dynamical equation for the mean field (f), which is the expectation value of 

^(x) = {n,m\^n{x)\n,m), (2.2) 

in a quantum state |r2, in) for which (f) is initially displaced from zero. In what follows, 
we shall assume that the quantum state initially corresponds to the vacuum state for 
the fluctuation field, defined as the difference between the Heisenberg field operator 
$H and the mean field, 

ip^{x) = ^nix)-4>ix). (2.3) 

It is important to note that even if the quantum state \Q, in) initially corresponds 
to the vacuum state for the fluctuation field, in a nonequilibrium setting (e.g., time- 
dependent background field (/>), it will not remain so. At later times, |il,in) will not 
correspond to the no-particle state for the fluctuation field. In what follows, we will 
simply refer to \n, in) as the initial "vacuum state," though it should be understood 
as the vacuum state for modes of the fluctuation field, and not of the field operator 
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An initial quantum state |0,in) with nonvanishing mean field such as described 
above may arise in the following way. Let us suppose that coupled to the scalar field 
(f) there is an external source^ which is nonvanishing for t < to, and which is removed 
at to, 

F{t) = F9{to-t). (2.4) 

Let us denote the quantum state for t < to by \Q, in), and suppose that in this state, 
for t <to, the expectation value of $h is given by a constant ^q. We will assume that 
the constant F satisfies 

F^-'-^lU (2.5) 

where Vgs is the vacuum cfi'cctive potential [151] for the theory with classical action 

(2.1), and that for t < to, \Q,in) corresponds to the vacuum state for the fiuctuation 

field $H — ^0- Then the expectation state ^ is equal to the constant stable equilibrium 

configuration (f)Q for t < to. The expectation value ^ is spatially homogeneous for all 

times due to the spatial translation invariancc of the fluctuation-field vacuum state 

in) and the action (2.1). Because of the instantaneous change in the external source 

at t = to, we may use the sudden approximation, in which |fi,in) is taken to be the 

initial quantum state for the t > to evolution. The physical picture here is that the 

expectation value of the scalar field operator is like a classical field initially held fixed 

at (j) = (j)Q for t < to and which is suddenly released at t = to- Let us ask whether 

there is an action r[(j)] whose variation gives the dynamical equation governing the 

subsequent evolution of the mean field 4>, including quantum corrections. As stressed 

above, because of the time-dependence of background mean field 4>, the condition 

[H,p] ^ for t > to, and in this case, the conventional "in-out" generating functional 

will not yi(^ld ( he coi rcct d\"iiainics of the mean field (f). Nevertheless, it is instructive 

^We assume the source is sufficiently weak tfiat we do not need to take into account nonperturbative 
Scliwinger-type particle production effects. 
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to see why this is so. 

In the conventional Schwinger-DeWitt or "in-out" approach [104, 129], one couples 
an arbitrary c-number source J (which is a function on M^) to the field ^ and computes 
the vacuum persistence amplitude in the presence of the source J. This amplitude has 

a path integral representation 



Z[J] = exp (^j^W[J]^ = J D(t>e^Y> ^ (s^[(p] + J d^xJ{x)(t){x) 



(2.6) 



where the functional integral is a sum over classical histories of the (f) field for which 
(f) — (j) is pure negative frequency [i.e., all spatial Fourier modes of 4> — ((> have a time 

dependence like exp(iujt), a; > 0] for t < to and pure positive frequency [~ exp{—iujt)] 
in the asymptotic future.^ The generating functional of normalized amplitudes is given 

by 

W[J] = -imnZ[J], (2.7) 

and it is well known that W[J] is the sum of all connected diagrams of the field theory 
in the presence of the source J [152]. The "classical" field 4>j = is a function of 
J defined by 

Assuming the functional relation (2.8) can be inverted to yield J in terms of ^, one 
can define an effective action [whose variation gives the inverse of Eq. (2.8)] as the 
Legendre transform of W, 

r[4>] = W[J] - J d^xJ{x)4>{x), (2.9) 

where we have dropped the J subscript on (f) for simplicity of notation. By differenti- 
ating Eq. (2.9), we find that 

ST 

= --J, (2.10) 
ocp 



^These boundary conditions on the functional integral are equivalent (up to an overall normalizar- 
tion) to adding a small imaginary term —itfjP' to the classical action, where e > 0. 
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where J is understood to be a functional of (p. The functional relation (j) = 8W/5J 
can indeed be inverted if V is well defined. The field equation satisfied by if) in the 
physical case of interest, J = F{t), is given by^ 

F(t) = -I34. (2.11) 

d(p{x) 

It should be emphasized that the T[(j)] appearing in Eq. (2.11) follows directly, by 
way of the Legendre transform, from the "in-out" generating functional defined in 
Eq. (2.6). The problem with Eq. (2.11) is that its solution, ^, is not the time-dependent 
expectation value of the Heisenberg field operator for the quantum field. Let us see 
why this is so, and what the interpretation of (p, as defined by Eq. 2.8, should be. 

At one loop (in Minkowski space) , the vacuum state is determined by an expansion 
of the fluctuation field operator in terms of spatial eigenmodes of the Klein-Gordon 
operator with time-dependent frequency 

a;| = P + m2 + ^0'. (2.12) 

In a nonequilibrium setting, such as in a curved or dynamical spacetime (where the 

scale factor is an additional time-dependent parameter in the effective frequency), or 

when (f) is time-dependent, the notion of positive frequency in the asymptotic past is in 

general different from that in the asymptotic future, in the present case because (p{t = 

00) 7^ (f)Q. Hence, the "in" vacuum state for the fluctuation field, |ri,in), and the "out" 

vacuum state for the fluctuation field, |0,out) [defined as the vacuum state for $h — 

(/)(t = 00)], are not equivalent. It is useful at this point to go over to an "interaction 

picture," where the "interaction" is the coupling Jcp to the external source. In this 

representation, the evolution of the field operator $h is just the Heisenberg evolution 

for the theory without the external source J (hence we retain the H subscript), and 

*While it is correct to include F{t) in the equation of motion as we have shown, one may simply 
set F{t) = in the equation of motion, since F{t) = ior t > to, and just incorporate F{to) = F into 
the initial data for S. Henceforth, this is what we shall do. 
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the evolution of a state vector is due to the coupling J'&h- In particular, the quantum 
state |r2, in) defined for t <tQ will evolve to a state \il,m)j{t) at time t, given by 

\n,in)j{t) = Texp(^ d^xJ{x)^u{x)] \n,m), (2.13) 
\riJxO<t J 

where T denotes temporal ordering. For convenience we will denote the t = oo limit 

of this state by 

|J^,in)j(+oo) = |l^,in)j. (2.14) 

In this interaction picture, the "in-out" generating functional takes a particularly 
simple form, 

Z[J] = {n,out\n,m)j{+oo) = (Q,out|Texp j d'^xJ{x)^^{x)^ (2.15) 

where T denotes temporal ordering. This amplitude is in general complex, even in the 
J = limit. The imaginary part of W[Q] gives the integrated probability P to produce 
a particle pair [31, 32] over the entire time range integrated in the classical action S, 

P = 2ImVF[0]. (2.16) 

It follows that the "classical field" (jyj defined as the functional derivative of W[J\ with 
respect to J, 

(n,out|$H(x)|0,in) = ^^ (2.17) 

is a matrix element which will in general be complex. In addition, the dependence of 
(0,out|$H(a^)|^^,in) on J will not, in general, be causal [65,66]. In curved spacetime, 
the quantum expectation value of the energy- momentum tensor, (T^i,), is obtained 
by functional differentiation of W with respect to g^'^ , which yields a complex matrix 
element of r^iy($H) between the "in" and "out" vacua, where Tjj^iy (6) is the clas- 
sical energy-momentum tensor for the field [17,67,129]. This is problematic because 
the quantum-corrected energy-momentum tensor, suitably regularized, constitutes the 
right-hand side of the semiclassical Einstein equation. 
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If instead of the above "in-out" approach one uses the closcd-timc-path formalism, 
one can obtain a real and causal evolution equation for the mean field ^ (the expecta- 
tion value of $h), as well as for the expectation value of the energy-momentum tensor. 
Here we briefly illustrate the procedure for the case of the theory in Minkowski 
space; the generalizations to curved spacetime and to higher correlation functions will 
be discussed in Sees. 2.4 and 2.3, respectively. Let x'^ = x° be far to the future of 
any dynamics we wish to study. It is not necessary to assume that A = or that 
the Hamiltonian is time independent at = x^. Here, we will specify initial condi- 
tions at x^ = — oo, for simplicity. As in the previous "in-out" approach, suppose we 
wish to compute the quantum-corrected equation governing the classical field (j). Let 
M = {{x^, x) \ — oo < x^ < x^} be the portion of Minkowski space to the past of time 
x^. We start by defining a new manifold as a quotient space 

X = (MX {+,-})/-, (2.18) 

where ~ is an equivalence relation defined by 

{x,+) {x',+) iffx = a;', 

(x,-) ~ (x',-) iSx = x', (2.19) 

{x, +) ^ {x' , — ) iff X = x' and x^ = 

It is straightforward to define an orientation on A4, provided we reverse the sign of 
the volume form between the -|- and — pieces of the manifold. Choosing an overall 
sign, we take the volume form on the -|- branch to be 

e = da;° A dx^ A dx'^ A dx^, (2.20) 

and the volume form on the — branch to be — e. The next step is to generalize the 
usual effective action construction of Eqs. (2.6)-(2.9) to the new manifold M. With 
the volume form on M thus defined, we can generalize the classical action (which 
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is a functional on M) to a functional on Ai as follows, 

S^[cf>+,<f>_] = S-[cf>+]-S-[cf>_], (2.21) 

where ^+ and ^_ denote the ^ field on the + and — branches of M , respectively. For 
a function (f) on M., let us define the restrictions of ^ to M x {+} and M x {— } by 
0+ and (p-, respectively. In order for to be a function on Ai, the restrictions must 
satisfy 

</.+ (x)Uo = </.-(x)Uo. (2.22) 

The generating functional of n-point functions (i.e., expectation values in the |i7,in) 
quantum state) for this theory is then defined by 



Z[J+, J_] = ^ Dct)+D<j>_ exp ^ (s^[ct)+,ct)-] + J^d\{J+ct)+ - J_(/._)) 



, (2.23) 



where J+ and J_ are c-number sources on the + and — branches of M, respectively. 
The designation "ctp" on the functional integral in Eq. (2.23) indicates that one in- 
tegrates over all field configurations (^+ , ^_ ) such that (i) ^+ = 0_ at the a;° = 
hypersurface and (ii) 0+ (<?-) consists of pure negative (positive) frequency modes at 
= —oo. It is not necessary for the normal derivatives of and 0_ to be equal at 
= x^ [67]. Because the theory is free in the asymptotic past,^ a positive frequency 
mode^ is a solution to the spatial-Fourier transformed Euler-Lagrange equation for 
whose asymptotic behavior at = -co is exp{—iiox^), for a; > 0. 

In the interaction picture where the sources J+ and J_ govern the state vector 
evolution, the expression Z[J+, J_] in (2.23) is seen to be the amplitude for the quan- 

®The vacuum boundary conditions for a theory which is not free in the asymptotic past are more 
comphcated, but can be treated by methods discussed in Sec. 2.5. 

®Here, the choice of vacuum boundary conditions corresponds to adding a small imaginary part 
ie{(f>+^ — cj>- ^) to the classical action . Alternatively, the boundary conditions correspond to the usual 
prescription — > — ie in S^[(p], but with <S^ now redefined as S^[4>+,4>-] = S^[4i+] — S^[(f)-]*, 
where ★ denotes complex conjugation [67]. 
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turn state to evolve forward in time under the source J-l from |i7,in) at x' 







-oo, to 



some arbitrary state IV') at x^ = x^, times the amphtude for the state at x° = x 



to evolve backwards in time under the source J_ to the state |0, in) at = — oo. The 
state must be summed over a complete, orthonormal basis of the Hilbert space of 
the (p field. In this picture, the CTP generating functional takes the form 



ZlJ+,J^ 



E 



{ft, in|rexp 



iTexp - 



i 








h 




J — oo 






I 


dx 




-oo 



dx° / d^xJ-^n{x) 



d'^xJ+^ii{x) |0,in) 



(2.24) 



where T and T denote temporal and anti-temporal ordering, respectively. The gener- 
ating functional for connected diagrams is then defined by 



W[J+,J. 



-zMnZ[J_|_, J_ 



(2.25) 



^aia2---an 



Classical fields on both -|- and — branches are then defined as 

where o, b are time branch indices with index set {-|-, — }. The matrix c"^ is the n = 2 
case of the n-index symbol c"i«2 - an (Jefined by 

1 if ai = 02 = • ■ • = o„ = +, 

-1 if oi = 02 = • • • = a„ = — , (2.27) 
otherwise. 

The functional differentiation in Eq. (2.26) is carried out with variations 6J+ and SJ- 
which satisfy the constraint that = 6J- on the x^ = x^ hypersurface. The J± 
subscript in Eq. (2.26) indicates the functional dependence on J±, which has been 
shown to be causal [65,66]. In the limit Jj^ = J- = J, the classical fields on the + 
and — time branches become equal. 



J+=J-=J 



J+=J-=J 



ct>{x)j= j(0,in|$H(x)|J^,in)^, (2.28) 
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where |r2,in)j, defined above in Eq. (2.14), is the state which has evolved from the 
vacuum at x° = — oo under the interaction ^h-'^, and Eq. (2.28) becomes the expecta- 
tion value ($h) in the limit J = 0. The effective action is again defined via a Legendre 
transform, with c"^ now acting as a "metric" on the internal 2x2 "CTP" space {+, — }, 

r[<^+, 0_] = W[J+, J_] - c"^ / d^xJa{x)^i,{x), (2.29) 

JM 

where the J subscripts on are suppressed for simplicity of notation. The func- 
tional dependence of J± on 4> via inversion of Eq. (2.26) is understood. By direct 

computation, tlic inverse of Eq. (2.26) is found to be 

^.(x),, = -J^^, (2.30) 

where we have indicated the explicit functional dependence of J± on with a sub- 
script, and Cah is the inverse of the matrix c"^ defined above. In the limit 0+ = = 0, 
Eq. (2.30) yields the evolution equation for the expectation value j($h)j = in the 
state which has evolved from |f2,in) under the source interaction J$h- The evolution 
equation for 0, the vacuum expectation value (fi, in|$H|^^5 in), is therefore 



5r[<^+,, 



= JI1 

4'+=4>-=<t> 



= 0. (2.31) 

<)>+=<t>-=4' 



Using Eqs. (2.30) and (2.29), an integro-differential equation for V can be derived [66], 
r[^^,^_]=-zMn|£ Z^<^^i^0_e^(^^'^-^-l-^-^^^'^(^"-^"))|, (2.32) 

in which the functional differentiations of V with respect to 0± are carried out with the 
constraint that the variations of 0± satisfy (5^+ = (50_ when x° = x^. The difference 

(pa — (pa is naturally interpreted as the fluctuations of a particular history (pa about 
the "classical" field configuration Let us, therefore, define the fluctuation field 

'^a^4>a- 4>a (2-33) 
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in analogy with Eq. (2.3. Performing the change of variables (pa functional 
integrand, we obtain 

T[hA-] = -^n^r.^l Z^^^Z^^_e*(^''[^-^-'^-^-'-^-'^^^'^-)|. (2.34) 
This functional integro-differential equation has a formal solution [152] 

r[<^+,<^_]=5^[,^+,^_]-^lndet(^-i) + ri[0+,0_], (2.35) 



where A^^{x,x') is the second functional derivative of the classical action with respect 
to the field 

iA'^'ix^x') = ff' [h,4>-]- (2.36) 

6(pa{x)(t)h{x ) 

The inverse of^»^is the one-loop propagator for the fluctuation field (p. The functional 
Fi in Eq. (2.34) is defined as —ih times the sum of all one-particle-irreducible vacuum- 
to- vacuum graphs with propagator given by A~^{x^ x') and vertices given by a shifted 
action 5;^^, defined by 

-\j d'x [ dV ( [^i]) ^aix)ip,ix'). (2.37) 

For simplicity of notation, we do not explicitly indicate the functional dependence of 
on (f)±. Figure 2.1 shows the diagrammatic expansion for Fi, where lines represent 
the propagator A~^{x,x'), and vertices terminating three lines are proportional to (j) 
[152]. Each vertex carries a spacetime label in M and a time branch label in {-|-, — }. 
The lowest-order contribution to Fi is of order h'^, i.e., a two loop graph. The one- 
loop propagator does not depend on h. The In(det^) term in Eq. (2.35) is the 
one-loop [order h] term in the CTP effective action. The CTP effective action, as a 
functional of (f)±, can be computed to any desired order in the loop expansion using 
Eq. (2.35). As with the ordinary in-out effective action [153], the CTP effective action 
contains divergences at each order in the loop expansion, which need to be regularized. 
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Figure 2.1: Diagrammatic expansion for the Fi part of the IPFCTP effective action 

The renormahzabiUty of the field theory in the "in-out" representation is sufficient to 
guarantee renormahzabiUty of the "in-in" equations of motion for expectation values 
[66-68]. 

Functionally differentiating F[0_|_,0_] with respect to either or and mak- 
ing the identification 4>+ = 4'- = 4' shown in Eq. (2.31)] yields a real and causal 
dynamical equation for the mean field 4>- Thus the IPI effective action F[^-|-] yields 
mean-field dynamics for the theory, which is a lowest-order truncation of the corre- 
lation hierarchy [81,82]. However, for a detailed study of nonpcrturbative growth of 
quantum fluctuations relevant to nonequilibrium mean-field dynamics (or a symmetry- 
breaking phase transition), it is also necessary to obtain dynamical equations for the 
variance of the quantum field, 

($1) - ($h)' = ($1) -4>' = i^i) = hG++{x,x), (2.38) 

where hG++{x,x') is the time-ordered Green function for the fluctuation field (fj^, 
{T{ip{x)iiip{x')ji)) . A higher-order truncation of the correlation hierarchy is needed 
in order to explicitly follow the growth of quantum fluctuations; the two-particle- 
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irreducible (2PI) effective action formalism, to which we now turn, serves this purpose. 

2.3 Two-Particle-Irreducible formalism 

In a nonperturbative study of nonequilibrium field dynamics in the regime where 
quantum fluctuations are significant, the IPI effective action is inadequate because 
it does not permit a derivation of the evolution equations for the mean field (<1>h) 
and variance {(f'^), at a consistent order in a nonperturbative expansion scheme. In 
addition, the initial data for the mean field (p do not contain any information about the 
quantum state for fluctuations around the mean field. The two-particle-irreducible 
(2PI) formalism can be used to obtain nonperturbative dynamical equations for both 
the mean field (p{x) and two-point function G{x,y), which contains the variance, as 
shown in Eq. (2.38). The 2PI method generalizes the IPI effective action to an action 
T[(f>, G] which is a functional of possible histories for both cf) and G. Alternatively, the 
2PI effective action can be viewed as a truncation of the master effective action to 
second order in the correlation hierarchy [82]. In this section we briefly review how 
the 2PI method works; more thorough presentations can be found in [68, 79]. 

Generalizing the IPI method where the mean field is fixed to be (j), the 2PI method 
fixes the mean field to be (p and the sum of all self-energy diagrams to be G. This leads 
to a drastic compression of the diagrammatic expansion containing the full information 
of the field theory [81]. Coupled dynamical equations for cj) and G are obtained by 
separately varying r[^, G] with respect to G and (f>. Imposing SF/S^ = yields an 
equation for the mean field </>, and setting 5T/6G = yields an equation for G, the 
"gap" equation. The variance (y?!) is the coincidence limit of the two-point function 
HG, as seen from Eq. (2.38). In a nonequilibrium setting, the closed-time-path method 
should be used in conjunction with the 2PI formalism in order to obtain real and causal 
dynamical equations for cp and G [68, 81, 84]. 
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Let us apply the 2PI method to a scalar theory in Minkowski space, with vac- 
uum initial conditions. In a direct generalization of Sec. 2.2, we now couple both a local 
source Ja{x) and nonlocal source Kab{x,x') (which are c-number functions on M) to 
the scalar field via interactions of the form hc°^Ja<i>a and hc"'^c'^Kac{x,x')<^b{x)<pd,{x')- 
Following Eq. (2.23), the CTP generating functional is defined as a vacuum persis- 
tence amplitude in the presence of the sources J and K, which has the path integral 
representation 

Z[J,K]= f D(j>+D(t)_e^p Us^[(j)+, f d^xd'^Jah 

+ 1 f d^X f d\'c''''c'"^Kacix,x')cl>b{x)Mx') 

(2.39) 

Here, is as defined in Eq. (2.21), and we are using Z[J,K] as a shorthand for 

Z[J^,J-;K^^,K ,K |_]. The generating functional for normalized n-point 

functions (connected diagrams) is defined by 

W[J, K] = -ihhiZ[J, K\. (2.40) 

The "classical" field 0„(x)j^ and two-point function Gab{x,x')jii are then given by 

2 , , 5W[J,K] 

hGab{x, x')jK = ^^<^^bdj^^^^^ - 4>a{x)jK(i>b{x')jK, (2.42) 

where we use the subscript JK to indicate that ^„ and Gab are functions of the sources 
J and K. 

In the limit K = J = 0, the classical field 0„ satisfies^ 

(f^+),/=A=o = {^-)j=K=o = {n\^n\n) = h (2.43) 



'^For simplicity, the "in" vacuum for the fluctuation field is henceforth denoted by \Q) instead of 
|f2,in), since we will not have any need to refer to the "out" vacuum in what follows. 
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i.e., (p^ becomes the expectation value of the Heisenberg field operator <I>h in the 
quantum state i.e., the mean field. In the same limit, the two-point function 
Gab is the CTP propagator for the fluctuation field defined by Eq. (2.3). The four 
components of the CTP propagator are, for J = = 0, 

nG++(x,x')|j=x=o = {^\T K(^)<^h(^')) I^^), (2-44) 

nG__(x,x')|j=i,=o = mf ^{x)^^{x')) (2.45) 

hG+.{x,x')\j=K=o = {9,\^^{x')^^{x)\a), (2.46) 

nG_+(x,x')|j=x=o = mv^{x)^^{x'm, (2.47) 

in the Heisenberg picture. In the coincidence limit x' = x, all four components above 
are equivalent to the variance ((^|) defined in Eq. (2.38), 

{<pl{x)) = hGab{x,x). (2.48) 

Provided we can invert Eqs. (2.41) and (2.42) to obtain J and K in terms of and 
G, the 2PI effective action can be defined as the double Legendre transform (in both 
J and K) ofW[J,K], 

r[^,G] = W[J,K] - f d^xc''^Ja{x)$b{x) 
Jm 

-\ I d\ I d^x'c''''c''''Kac{x,x')[hGu{x,x') + 4>^{x)^a{x)\. (2.49) 

As with W[J,K] above, r[0,G] denotes r[^+, 0_; The JK 
subscripting of (p and G is suppressed, but the functional dependence of and G on J 
and K through inversion of Eqs. (2.41) and (2.42) is understood. By direct functional 
differentiation of Eq. (2.49), the inverses of Eqs. (2.41) and (2.42) are found to be 



(5^„(x) 



(2.50) 

^111^ (2.51, 



29 



where the subscript "ijiG" indicates that K and J are functionals of and G. Once 
r[0, G\ has been calculated, the evolution equations for ^ and G are given by 

5T[lG] 

4.+-- 

= 0. (2.53) 



sniG] 



= 0, 



(2.52) 



SGabix,y) 

Of course, the two equations contained in Eq. (2.52) (corresponding to a = + and 
a = —, respectively) are not independent, just as in Eq. (2.31). In addition, only two 
of equations (2.53) are independent, one on the diagonal and one off diagonal in the 
"internal" CTP space. Using both Eq. (2.49) and Eq. (2.39), an equation for T[^,G] 
in terms of the sources K and J can be derived. 



r[(f),G] = -imn{ I D(P+D(t)_exp 

J ctp 



/< 

JM 



-(5^[0+,0_]+c'^'' / d''xJa{x)[<t>^{x)-(t^,{x)] 



f d^x [ .^.X.(.,.O[^e(.)0.(.O - 4(.)^.(.') - nGUx,x')^) 



(2.54) 

The sources K and J in the right-hand side of Eq. (2.54) are functionals of 0, through 
Eqs. (2.50) and (2.51). Expressing this functional dependence, we obtain a functional 
integro-differential equation for F, 



mG] = f d'x f d'x' /^^^'^\ Gab{x,x') 
Jm Jm SGba[x,x) 



SGbaix ,X 

— ihln-i I D(j)-^-D(j)- exp 



ctp 



4 sn<P,a] 

X 

Ha 



'-'a - VaJ 



\l d'xj d'x' [Ux) - k{x)][Mx') - k{x')]) 
nJM Jm (>Gba{x ,x) / 



(2.55) 



We have dropped the JK subscripting because the functional derivatives in the equa- 
tion are only with respect to and G. As in Sec. 2.2, a change of variables D4>± — Dip^ 
is carried out in the functional integral, with the fluctuation field defined by 
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Eq. (2.33). The resulting equation 

Jm J Ad ^Gba{x',x) 



I 

*f ctp 



iYihy\ \ Difj^Dif - exp 



S''[^+ + (t)+,^- +0-] 



Jm 56„ nJM Jm oGbaix ,x) 



(2.56) 



has the formal solution [79] 



T[4>,G\=S^[4>\-'hudei{Gab) 

ih 
+ 2 



d^x I d^x'A''\x',x)Gab{x,x') + r2[^,G], (2.57) 
Jm 



IM 

where is the second functional derivative of the classical action 5^, evaluated at 0^ 
[as defined in Eq. (2.36)]. The functional r2 is —ih times the sum of all two-particle- 
irreducible vacuum-to- vacuum diagrams with lines given by Gab and vertices given by 
a shifted action S^^^ defined by Eq. (2.37). The shifted action for the scalar field 
theory is 

SLM = - SU<p^], (2.58) 

in terms of 

SfM = Ij'x (^\^' + , (2.59) 

where the functional dependence of 5^,^^ on (f)± is not shown explicitly. Two types 
of vertices appear in Eq. (2.59): a vertex which terminates four lines and a vertex 
terminating three lines which is proportional to the mean field (j). The expansion 
for r2 in terms of G and is depicted graphically up to three-loop order in Fig. 2.2, 
where lines represent the propagator G and vertices are given by 5;^^ [79] . The vertices 
terminating three lines are proportional to (p. Each vertex carries a spacetime label 
in M and a time branch label in {-|-, — }. In general, the 2PI effective action contains 
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Figure 2.2: Diagrammatic expansion for r2 part of the 2PI-CTP effective action 

divergences at each order in the loop expansion. It has been shown that if a field 
theory is renormalizable in the "in-out" formulation, then the "in-in" equations of 
motion derived from the 2PI effective action are renormalizable [68]. We will discuss 
the renormalization of the 2PI effective action below in Sec. 2.5.4. 

Various approximations to the full quantum dynamics can be obtained by trun- 
cating the diagrammatic expansion for Throwing away r2 in its entirety yields 
the one-loop approximation. In Fig. 2.2, there are two two-loop diagrams, the "set- 
ting sun" and the "double-bubble." Retaining just the double-bubble diagram yields 
the time-dcpcndcnt Hartree-Fock approximation [79]. Retaining both diagrams gives 
the two-loop truncation of the theory.® This approximation will yield a non-time- 
reversal-invariant mean-field equation above threshold, due to the setting sun diagram 
[82]. 

The time-rcvcrsal noninvariance of the mean-field and gap equations generated 

by the 2PI effective action is a consequence of the fact that the 2PI effective action 

really corresponds to a further approximation from the two-loop truncation (in the 

sense of topology of vacuum graphs) of the master effective action [82] . The two-loop 

truncation of the master effective action is a functional r;=2[0)G, C3] which depends 

on the mean field <j), the two-point function G, and the three-point function C3. The 

*A different approximation, the 1/N expansion, is used in Sec. 2.5 to study the nonequihbrium 
dynamics of the 0(iV) field theory. 
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four-point function C4 also appears, but is not dynamical due to a constraint. The 
full set of equations, 



Ha 
SGab 

dri=2[lG,Cs] 



0, (2.60) 

= 0, (2.61) 
= 0, (2.62) 



is time-reversal invariant. However, the 2PI effective action is obtained by solving 
Eq. (2.62) with a given choice of causal boundary conditions and substituting the 
resulting C3 into Ti=2, to obtain r2[</>,G]. This "slaving" of C3 to (p and G with a 
particular choice of boundary conditions is what breaks the time-reversal invariance of 
the theory [82]. This is shown explicitly below in Chapter 5. In Chapter 3 where we 
discuss the preheating dynamics of the inflaton field, we work with further approxima- 
tions which discard the setting sun diagram, and thus regain time-reversal invariance 
of the dynamical equations. 

2.4 A^"* field theory in curved spacetime 

In this section the quantum dynamics of a scalar field theory is formulated in 
semiclassical gravity, which means that the matter fields are quantized on a curved 
classical background spacetime.^ The two-particle-irreducible effective action is used in 
conjunction with the CTP formalism to obtain manifestly covariant, coupled evolution 
equations for the mean field ($h) and variance ($|) — ($h)^ in the model. 

Let us consider the scalar field theory in a globally hyperbolic, curved back- 
ground spacetime with metric tensor Qfj^i,. The diffeomorphism-invariant classical ac- 

®The semiclassical approximation is consistent with a truncation of the quantum effective action 
for matter fields and gravity perturbations at one loop [i.e., at order 0{h)] [31] or [in the case of the 
0{N) field theory studied here] at leading order in the 1/N expansion [37, 149]. 
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tion for this system is 

5[<^,5n = '5°[5n+5'['^,5n, (2.63) 

where g*^" is the contravariant metric tensor, and 5° and 5"^ are the classical actions 
of the gravity and scalar field sectors of the theory, respectively. For the scalar field 
action, we have 



S^[ct>,gn = -\ jd'x^g 



(2.64) 



where ^ is the (dimensionless) coupling to gravity (necessary in order for the field 
theory to be renormalizable [154]), □ is the Laplace-Beltrami operator defined in 
terms of the covariant derivative V^, and R is the scalar curvature. The constant 
m has units of inverse length, and the self-coupling A has units of 1/h. Following 
standard procedure in semiclassical gravity [17], we define the semiclassical action for 
gravity to be 

S"" [^'^l = J d^x^ [r-2K + cR" + bR'^^Rap + aR'^f^^^Ra^^s] , (2.65) 

where a, b, and c are constants with dimensions of length squared, RajS-fS is the Rie- 
mann tensor, R^^p is the Ricci tensor, A is the "cosmological constant" (with units of 
inverse length-squared) , is the square root of the negative of the determinant of 
g^j,i,, and G (with units of length divided by mass) is Newton's constant. As a result 
of the generalized Gauss-Bonnet theorem [155], the constants a, 6, and c arc not all 
independent in four spacetime dimensions; let us, therefore, set a = 0. Classical Ein- 
stein gravity is obtained by setting 6 = and c = 0. Minimal and conformal coupling 
(for the ^ field to gravity) correspond to setting ^ = and ^ = 1/6, respectively. 

The motivation for including the arbitrary coupling ^ and the higher-order cur- 
vature terms R^ and R^^Ra/s in the classical action S is that we wish to study the 
semiclassical dynamics of the theory. In the semiclassical gravity field equation and 
matter field equations, divergences arise which require a renormalization of b, c, G, A, 
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m, ^, and A [17]. These quantities, as they appear in Eqs. 3.1-2.65, are understood to 
be bare; their observable counterparts are renormahzed. 

The classical Euler-Lagrange equation for ^ is obtained by functionally differenti- 
ating S^[(f),g'^''] with respect to (f), and setting SS^ /S(f) = 0, 

n + m'^ + ^R + ^(t)'^^ = 0. (2.66) 

The Euler-Lagrange equation for the metric g^i, is obtained by functional differen- 
tiation of S with respect to gi^" (it is assumed that the variations 5(f) and 5g'^'^ are 
restricted so that no boundary terms arise) , 

G^,, + Ag^u + c + b = -SttGT^,, (2.67) 
where the tensors Gj^n,, ^^^Hj^^i,, and ^"^^Hj^ii, are defined by [101, 102] 

= 2R,^, + 2g^,nR - ^g^^R^ + 2RRf,,, (2.68) 

= 2i?/,„ - Di?^, - ^g^^^DR + 2R^^Ra, - ^g^^uR^^R^fS, (2-69) 

Gij,i, = Rfiv — -jRgixv (2.70) 
In Eq. 2.67, T^^ is the classical energy-momentum tensor, 

T^, = (1 - 20<p,^<p,,+ (2i - 5m.5^"0;p0;. - 2e0;^.0 + 2ig^,4>U(t> 



(2.71) 



We are interested in the dynamics of the expectation value of the scalar field operator 
and its higher moments, which in nonequilibrium field theory does not follow directly 
from functional differentiation of the usual Schwinger-DeWitt or "in-out" effective 
action. Instead, as discussed above in Sec. 2.2, the Schwinger-Keldysh formalism 
should be used. Here we discuss the implementation of the Schwinger-Keldysh method 
in curved spacetime. 

The first step is to generalize the closed-time-path (CTP) manifold 7W, defined in 
Eq. (2.18), to curved spacetime. Let be a Cauchy hypersurface chosen so that its 



35 



past domain of dependence [156], Z)_(S^), contains all of the dynamics we wish to 
study. Let us then define the manifold (with boundary) 

M = D_(S^). (2.72) 

The CTP manifold M is defined following Eq. (2.18) as a quotient space constructed 
by identification on the hypersurface C dM as in Eq. (2.18) 

M = {Mx {+,-})/ ^, (2.73) 

where the equivalence relation is the same as Eq. (2.19) except that the matching of 
+ and — time branches is now done on S^. M x {+, — } defined by 

(x,+) ~ (a:',+) iffx = a;' 

(x,-)-(x',-) iSx = x' (2.74) 
(x, +) ~ {x', — ) iff a; = x' and x G S*. 

We construct an orientation on M. using the canonical volume form from M, em, 

Cm = dx'^ A dx^ A dx^ A dx^ , (2-75) 

and define the volume form on A4 to be 

{Cm on M X {+}, 
(2.76) 
-CM on M X {-}. 

Finally, we let (p and g'^'^ be independent on the + and — branches of M, provided 
that g'^ = g'^ and ^+ = 0- on S^. In other words, <p and gi^" must be a scalar and 
a tensor, respectively, on M. In terms of the volume form em, we can write a scalar 
field action on M, 

<S^[</.±,5r] = S^^+,g^''] - S^^^,gn, (2-77) 

where S^[(l)\ is given by Eq. (2.64), and g'^ is the metric tensor on the + and — 
branches of M. Using Eq. (2.65) we can similarly define the gravity action S'^ on M, 

,5^] = S^lgl""] - S^lg^, (2-78) 
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where it is understood that only configurations of g!^ satisfying the constraint = 
g^ on Syr are permitted. 

In semiclassical gravity the scalar field theory (with action S^) is quantized on a 
classical background spacetime, with metric g^^, whose dynamics is determined self- 
consistently by the semiclassical gcomctrodynamical field equation. Let us denote the 
Heisenberg-picture field operator for the canonically quantized ^ field by $h- We wish 
to compute the quantum effective action F for this scalar field theory, using the two- 
particle-irreducible (2PI) method described in Sec. 2.3. In terms of (now defined 
on the curved manifold A^), we define the 2PI, CTP generating functional Z[J, g^^\ 
as follows: 

+ \l d^xV^ [ d^x'J^^,c-^'''c'^'''Kac{x,x')ct>,{x)Ux') 

(2.79) 

where we have written Z[J,K,g^^^] as a shorthand for Z[J±^K±±,g'^]. The three- 
index symbol c"**^ is the n = 3 case of the n-index symbol -an defined by Eq. (2.27). 
The boundary conditions on the functional integral define the initial quantum state 
(assumed here to be pure). In this and in Chapters 3-4, we are interested in the case 
of a quantum state corresponding to a nonzero mean field 0, with vacuum fluctuations 
around the mean field. This entails a definition of the vacuum state for the fluctuation 
field, defined in Eq. (2.3). In a general curved spacetime, there does not exist a unique 
Poincare-invariant vacuum state for a quantum field [129, 130]. For an asymptotically 
free field theory, a choice of "in" vacuum state corresponds to a choice of a particular 
orthonormal basis of solutions of the covariant Klein-Gordon equation with which to 
canonically quantize the field. 

From Eq. (2.79), we can derive the two-particle-irreducible (2PI) effective action 
T[(f)^G,g^^] following the method of Sec. 2.3, with the understanding that F now 
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depends functionally on the metric g^^ on the + and — time branches. The functional 
differentiations should be carried out using a covariant generalization of the Dirac 5 
function to the manifold M [17]. 

1 



Jm ^/-g{x') 



= 1. 



(2.80) 



The functional integro-differential equation (2.56) for the CTP-2PI effective action can 
then be generalized to the curved spacetime M in a straightforward fashion, modulo 
the curved-spacetime ambiguities in the boundary conditions of the functional integral 
(2.79). 

The (bare) semiclassical field equations for the variance, mean field, and metric can 
then be obtained by functional differentiation of S°[g^'^] + T[(p,G, g^'^] with respect 
to G±±, (j)±, and g^^ , respectively, followed by metric and mean-field identifications 
between the + and — time branches, 

S{S°[g^'']+T[4>,G,gn) 



Sgii" 
5r[4>,G,gf^'^] 



Ha 

5V[lG,gn 



SGab 



0, 



0, 



0. 



(2.81) 
(2.82) 
(2.83) 



As above, CTP indices are suppressed inside the functional arguments. Eqs. (2.81), 
(2.82), and (2.83) constitute the semiclassical approximation to the full quantum dy- 
namics for the system described by the classical action 5° [g'^'^\ + g'^'^] . The 
semiclassical field equation (with bare parameters) for g^^'^ is obtained directly from 
Eq. (2.81), 

G^,u + Ag^, + c + b = -87rG(T^,), (2.84) 

where (T^i/) is the (unrenormalized) energy-momentum tensor defined by 

2 (5T[lG,gn 



'-9 



6g 



(2.85) 



<!>+-- 



9+ =91 
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Equation (2.84) gives the spacetime dynamics; the dynamics of (/) and G arc given 
by the mean-field (2.83) and gap (2.82) equations. In Eq. (2.85), the angle brackets 
denote an expectation value of the energy-momentum tensor (with the Heisenberg 
field operator $h substituted for cf) in the classical theory) with respect to a quantum 
state \Q) defined by the boundary conditions on the functional integral in Eq. (2.79). 
In four spacetime dimensions the unrenormalized (T^i/) has divergences which can be 
absorbed by the renormalization of G, A, b, and c [17]. It is often useful to carry out 
the renormalization in the evolution equations for expectation values rather than in 
the CTP effective action [66]. 

The energy-momentum tensor (T^i/) is obtained by variation of the 2PI effective 
action T, which is a functional of the metric on both the -|- and — time branches. 
Prom Eq. (2.79), it is possible to derive r[(j),G, gi^"] as an arbitrary functional of g'^ 
and g'^'^ . However, in practice it is often easier to fix the metric to be the same on 
both the -|- and — time branches, i.e., 

5f = 5^^^^^'^^ (2.86) 

in the computation of r[(/), G, g^'^]. Once T[(/), G, g^^\ (or some consistent truncation of 
the full quantum effective action for S^) has been computed, it is then straightforward 
to determine how V[^, G, g'^"] should be generalized to the case of an arbitrary metric 
on M, for which g^ and g^ are independent. The bare energy-momentum tensor 
(T^u) can then be computed using Eq. (2.85). Accordingly, in Sec. 2.5 below, we fix 
g^l" = = g""" in the calculation of r[4>,G, g'"'']. 

The semiclassical Einstein equation is a subcase of the general geometrodynamical 
field equation (2.84), obtained (after renormalization) by setting the renormalized 
6 = c = A = [17]: 

G^, = -87rG(r^,), (2.87) 
where we assume that the cosmological constant vanishes. Having shown how to 
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derive coupled evolution equations for the mean field, variance, and metric tensor in 
semiclassical gravity, we now turn our attention to the scalar 0{N) model in curved 
spacetime. 

2.5 O(A^) field theory in curved spacetime 

In this section we derive coupled nonperturbative dynamical equations for the mean 
field (p and variance ((/?|) for the minimally coupled 0{N) scalar field theory with quar- 
tic self-interaction and unbroken symmetry. The background spacetime dynamics is 
given by the semiclassical Einstein equation. These equations self-consistently account 
for the back reaction of quantum particle production on the mean field, and quantum 
fields on the dynamical spacetime. In our model the Heisenberg-picture quantum state 
\4>) is a coherent state for the field $h at the initial time r/o, in which the expectation 
value ($h) is spatially homogeneous. The coherent state is defined with respect to the 
adiabatic vacuum constructed via matching of WKB and exact mode functions for the 
fluctuation field in some asymptotic region of spacetime. 

The 0(iV) field theory has the property that a systematic expansion in powers 
of 1/7V yields a nonperturbative reorganization of the diagrammatic hierarchy which 
preserves the Ward identities order by order [148]. Unlike perturbation theory in the 
coupling constant, which is an expansion of the theory around the vacuum config- 
uration, the 1/N expansion entails an enhancement of the mean field by ^/N■, this 
corresponds to the limit of strong mean field. (This is precisely the situation which 
can arise in chaotic inflation at the end of the slow-roll period, where the inflaton mean 
field amplitude can be as large as Mp/3 [125].) As discussed in Sees. 2.2 and 2.3, the 
nonequilibrium initial conditions for the mean field as well as the nonperturbative as- 
pect of the dynamics requires use of both closed-time-path and two-particle- irreducible 
methods. The 1/N expansion can be achieved as a further approximation from the 
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two-loop, two-particle-irreducible truncation of the Schwinger-Dyson equations. 

Although in this study we assume a pure state, the 2PI formalism is also useful 
for an open system calculation, in which the mean field is defined as the trace of 
the product of the reduced density matrix p and the Heisenberg field operator $h, 
Tr(p<I>H). When the position-basis matrix clement {(j)i\p{riQ)\(j)2) can be expressed as 
a Gaussian functional of and ^2, the nonlocal source K can encompass the initial 
conditions coming from p(fo) in a natural way [68]. In order to incorporate a density 
matrix whose initial condition is beyond Gaussian order in the position basis, one 
should work with a higher-order truncation of the master effective action [82]. This 
we will do for the theory in Chapter 5. The leading-order 1/N approximation is 
equivalent to assuming a Gaussian initial density matrix; therefore, the 2PI effective 
action is adequate for our purposes. 

2.5.1 Classical action for the 0{N) theory 

The 0(iV) field theory consists of N spinless fields (f) = {(f)^}, i = 1, . . . ,N, with an 
action which is invariant under the A?^-dimensional real orthogonal group. The 
generally covariant classical action for the 0(iV) theory (with quartic self- interaction) 
plus gravity is given by 



where S'^[g^'^] is defined in Eq. (2.65) for the spacetime manifold M with metric g^i,, 
and the matter field action S^[<p^,g'^^] is given by 



(2.88) 



M 



(2.89) 
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The 0(iV) inner product is defined by^° 

$.$=(l)'(fPSij. (2.90) 

In Eq. (2.89), A is a (bare) coupling constant with dimensions of 1/h, and ^ is the 
(bare) dimensionless couphng to gravity. The classical Euler-Lagrange equations are 
obtained by variation of the action S separately with respect to the metric tensor 
g^iy and the matter fields (f)^. In the classical action (2.89), the 0(iV) symmetry is 
unbroken. However, the 0{N) symmetry can be spontaneously broken, for example, 
by changing tv? to —tv? in . In the symmetry-breaking case with tachyonic mass, 
the stable, static equilibrium configuration is found to be 

f.f^^Jfl^^^, (2.91, 

which is a constant. If we wish to study the action for small oscillations about the 
symmetry-broken equilibrium configuration, the 0(N) invariance of Eq. (2.89) implies 
that we can choose the minimum to be in any direction; we choose it to be in the first, 
i.e., {4>^)'^ = v^. In terms of the shifted field a = 4?- — v and the unshifted fields (the 
"pions") tt' = i = 1, . . . , A'' — 1, the action becomes 



C7(n + + iR)a + 7? • (□ + + ^E)?? 



(2.92) 



One can show that the effective mass of each of the "pions" vf (defined as the second 
derivative of the potential) is zero, due to Goldstone's theorem. The theorem holds 
for the quantum-corrected effective potential as well [153]. In this dissertation we 

^"in our index notation, the latin letters i,j,k,l,m,n are used to designate 0(JV) indices (with 
index set {1, ...,7V}), while the latin letters a,b,c,d,e, f are used below to designate OTP indices 
(with index set {-|-, — }). 
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study the unbroken symmetry case, in order to focus on the parametric amplifica- 
tion of quantum fluctuations; this avoids the additional complications which arise in 
spontaneous symmetry breaking, e.g., infrared divergences [37,40, 84]. 



2.5.2 Quantum generating functional 

In this section we derive the mean-field and gap equations at two-loop order. The 
2PI generating functional for the 0{N) theory in curved spacetime is defined using 
the closed-time-path method in terms of c-number sources J* and nonlocal c-number 

sources K^Jj^ on the CTP manifold Ai, 

Z[Jl K%, g^,] =11 [ DcPl exp ^ f 5^ [</>, + / d'^xV^c'^'Ja ■ h 

+ \f d*x^ [ d^x'^c-'c'=''K^J,{x,x')4{x)cp'a{x%kSj) 
^ Jm Jm j 

(2.93) 

where we have (as discussed above) suppressed time branch indices on the metric 
tensor, and the CTP classical action is defined as in Eq. (2.77), 

The sources are coupled to the field by the 0{N) vector inner product 

Ja-$b = Ji'PiSij. (2.95) 

For simplicity, we shall suppress all indices inside functional arguments. In addition, 
the boundary conditions on the asymptotic past field configurations for cp'tj. in the 
functional integral correspond to a choice of "in" quantum state |(?!)) for the system. 
The generating functional for normalized expectation values is given by 

W[J, K, g^,] = -ihlnZlJ, K, g^,] , (2.96) 
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with the additional functional dependence of both W and Z on g^'^ understood. In 
terms of this functional, we can define the "classical" field (p and two-point function 
G by functional differentiation, 

= ^J}^5^J ^ (2.97) 



U^mx') + hG^,{x,x') = 2^-^—,^^5V. (2.98) 



In the zero-source limit iC*], = J* = 0, the classical field 0„ satisfies 

{4^+)j=K=o = {t)j=K=o = mV\(t>) ^ f (2.99) 

as an expectation value of the Heisenberg field operator $g in the quantum state 
The fluctuation field is defined [following Eq. (2.3)] in terms of the Heisenberg field 
operator $h and the mean field <j) (times the identity operator), 

(^j, = - ^\ (2.100) 

In the same limit J = K = 0, the two-point function G^^^ becomes the CTP propagator 
for the fluctuation field. The four components of the CTP propagator are (for J* = 

(2.101) 
(2.102) 
(2.103) 
(2.104) 

In the coincidence limit x' = x, all four components (2.101)-(2.104) are equivalent to 

the mean-squared fluctuations (variance) about the mean field 4> , 

n&l^{x,x\j=K=o = = («)')• (2.105) 





^')\J- 


=K=0 


= {n\T{^i{x)<fl,{x'))\n), 


nG'i_{x, 




=K=0 


= mf{if\{x)ipi{x')m 


hG%_{x 


,x')j= 


=K=0 


= {n\ipiix')ipiix)\n), 


hG%{x 




=K=0 


= {n\^i,{x)^Ux'm)- 
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Provided that Eqs. (2.97)-(2.98) can be inverted to give J* and in terms of (f)^ 



and G^-J,, we can define the 2PI effective action as a double Legendre transform of 



^ JM Jm 



/< 



-g'cVK^J^ix, x') \hGt!i{x, x') + SikSji, 



(2.106) 

where J* and K^J^ above denote the inverses of Eqs. (2.97) and (2.98). Prom this 
equation, it is clear that the inverses of Eqs. (2.97) and (2.98) can be obtained by 
straightforward functional differentiation of P, 



(5P 



„ab 



(2.107) 

j_ .^^ ^ = --c'^c'^^Kl^ix', x). (2.108) 

Performing the usual field shifting involved in the background field approach [152], it 
can be shown that the 2PI effective action which satisfies Eqs. (2.106), (2.107), and 
(2.108) can be written 



r[,^,G,5n=5"[<^,5n-ylndet 
ih 



+ - / d^x^g dV^^«j'(x',x)G^4(x,x')+r2[0,G,5n, (2.109) 

^ JM JM 

where the kernel A is the second functional derivative of the classical action with 
respect to the field cf), 

1 / S'^S" ..,.\ 1 



^•^ij {XiX ) 



,9' 



(2.110) 



v^\5<i)i{x)<fiix'y '■ -jv-9' 

and r2 is a functional to be defined below. Evaluating Afj by differentiation of 
Eq. (2.89), we find 



iAf^{x,x') = -|(5ijc"^p^ + + ^R{x) 



+ 2iV^' 



flbcd 



[^tix)^Ux)]SijdM + 24>lix)^'aix)5ikSji\ y\x - x')^, (2-111) 
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where the four-index symbol c"*"^"^ is defined by Eq. (2.27). In Eq. (2.109), the symbol 
r2 is defined as —ih times the sum of all two-particle-irreducible vacuum-to-vacuum 
graphs with propagator G and vertices given by the shifted action S^^^, which takes 
the form 

Jm \sk J 

1 /■ f / r2cF \ {Z.LIZ) 

-- d^x d^x'l U,gn]^i{x)ipi{x'). 

•^Jm Jm V(^<^i(x)</>^,(x') V ' 

The expansion of r2 in terms of G and (f) is depicted graphically in Fig. 2.2 for the 

theory. In the present case of the 0{N) field theory, each vertex now also carries 

an 0(N) label. From Eqs. (2.112) and (2.89), S^^^ is easily evaluated, and we find 

SL[v,9n = SL[v^+,gn - SfJ^.,gn (2.113) 

"1 



(2.114) 



The two types of vertices in Fig. 2.2 are readily apparent in Eq. (2.114). The first term 
corresponds to the vertex which terminates four lines; the second term corresponds to 
the vertex which terminates three lines and is proportional to 4>- 

The action T including the full diagrammatic series for r2 gives the full dynamics 
for (f) and G in the O(A^) theory. It is of course not feasible to obtain an exact, 
closed-form expression for r2 in this model. Instead, various approximations to the 
full 2PI effective action can be obtained by truncating the diagrammatic expansion 
for r2. Which approximation is most appropriate depends on the physical problem 
under consideration. 

1. Retaining both the "setting-sun" and the "double-bubble" diagrams of Fig. 2.2 
corresponds to the two- loop, two-particle-irreducible approximation [82]. This 
approximation contains two-particle scattering through the setting-sun diagram. 

2. A truncation of r2 retaining only the "double-bubble" diagram of Fig. 2.2 yields 

equations for (p and G which correspond to the time-dependent Hartree-Fock 
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approximation to the full quantum dynamics [79, 117]. This approximation does 
not preserve Goldstone's theorem in the symmetry-breaking case, but it is energy 
conserving (i.e., leads to a covariantly conserved energy-momentum tensor) [117]. 

3. Retaining only the (TrG^-'^)^ piece of the double-bubble diagram corresponds to 
taking the leading order 1/N approximation, as will be shown below in Sec. 2.5.3. 

4. A much simpler approximation consists of discarding r2 altogether. This yields 
the one-loop approximation. The limitations of the one-loop approximation for 
nonequilibrium quantum field dynamics have been extensively documented in 

the literature [68, 81, 82, 84]. 

Let us first evaluate the 2PI effective action at two loops [37, 54] . This is the most 
general of the various approximations described above, i.e., both two- loop diagrams 
in Fig. 2.2 are retained. The 2PI effective action is given by Eq. (2.109), and in this 
approximation, r2 is given by 



field and gap equations, respectively. The two-loop gap equation is given by 




\c^^'^'' [ d^'xV^iGlix, x)G,^^(x, x) + 2Gf,ix, x)G^{x, x)\5,,6ki 
d'xV^ I dV^<^:(x)<^:',(a:')[G^i:.(a:,x')Gf;(x,x')G^^:(:c,x') 

JM 

+ 2G;{, (x, x')g2' (x, x')G% (x, x')\5,k5rk' ■ 



(2.115) 



Functional differentiation of T[(l),G,g^'^] with respect to (j) and G leads to the mean- 



(G-i)fj'(x, x') = AfAx, x') + —c^^^'S^x - x') x) + 2G;^,(x, x) 
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+ 20:(x)4(a;')G^i (x, rcOG^^,, (x, x') 

+ 24>%)4>Ux')G2{x,x')GfA^,x') 

+ 20'(x)<^',,(x')G^,i (x, x')GiUx, x') 



(2.116) 



The mean-field equation is found to be 



4iV2 ./M 



+ 



2N 



5,j<f>, GZ{x, x) + 5^i5TU(^l,{x, x) + x)] = 0, (2.117) 



where the nonlocal function S'^(a:,x') is defined by 

S^-(x, xO = c^^'-'^c"'^'^''^' <^:, (xO [G^f (x, x')G^i; (x, x')G^^, (x, x') 

+ 2G^f(x,x')G2'(x,x')G^^(x,x') 

+ G^:^(x',x)G^X(^''^)^^d'd(^''^) 

+ 2G%{x', x)Gf;{x', x)G^^(x', x)] SjkSj^k'Sw 

(2.118) 

Taking the limit = 4>^- = 4>' i^i Eqs. (2.116) and (2.117) yields coupled equations 
for the mean field 4>^ and the CTP propagators G^-j,, on the fixed background space- 
time g^'^ . The equations, as well as the semiclassical Einstein equation obtained by 
differentiating r[(f),G, g'^'^] with respect to g'^'^, are real and causal, and correspond 
to expectation values in the limits (f>+ = (f>- = ^. The O(A^) parts of the above 
equations are nonlocal and dissipative. The nonlocal aspect makes numerical solution 
difficult; the dissipative aspect will be addressed below in Chapter 4. One can regain 
the perturbative (amplitude) expansion for the CTP effective action at two loops by 
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expanding the one-loop CTP propagators in Eq. (2.117) in a functional power series 
in the mean field (f). 

2.5.3 Large- approximation 

We now carry out the 1/A^ expansion to obtain local, covariant, nonperturbative mean- 
field and gap equations for the 0(A'^) field theory in a general curved spacetime. The 
\/N expansion is a controlled nonperturbative approximation scheme which can be 
used to study nonequilibrium quantum field dynamics in the regime of strong quasiclas- 
sical field amplitude [99, 108, 117, 118]. In the large-A^ approach, the large-amplitude 
quasiclassical field is modeled by A'" fields, and the quantum-field-theoretic generating 
functional is expanded in powers of 1/N. In this sense the method is a controlled 
expansion in a small parameter. Unlike perturbation theory in the coupling con- 
stant A, which corresponds to an expansion of the theory around the vacuum, the 
large- A'^ approximation corresponds to an expansion of the field theory about a strong 
quasiclassical field configuration [117]. At a particular order in the l/N expansion, 
the approximation yields truncated Schwinger-Dyson equations which are 0{N)- and 
renormalization-group-invariant, unitary, and (in Minkowski space) energy conserving 
[117]. In contrast, the Hartree-Fock approximation cannot be systematically improved 
beyond leading order, and in the case of spontaneous symmetry breaking, it violates 
Goldstone's theorem and incorrectly predicts the order of the phase transition [118]. 

Let us implement the leading order large- A^ approximation in the two- loop, 2PI 
mean-field and gap equations (2.117) and (2.116) derived above. This amounts to 
computing the leading-order part of F in the limit of large N , which is 0{N). In the 
unbroken symmetry case, this is easily carried out by scaling by \/iV and leaving G 
unsealed [79], 

<^l(x) ^ ViV0„(x), (2.119) 
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a'l^{x,x')^Gab{x,x')5'^, 



(2.120) 
(2.121) 
(2.122) 



for all The Heisenberg field operator (fl^ scales like (fl^ in Eq. (2.122). In the above 
equations, the connection between the large-A'^ limit and the strong mean-field limit 
is clear. 

The truncation of the 1/N expansion should be carried out in the 2PI effective 
action, where it can be shown that the three-loop and higher-order diagrams do not 
contribute (at leading order in the 1/N expansion). Let us now also allow the metric 
5^1^ to be specified independently on the -|- and — time branches. We find, for the 
classical action, 

where 



(2.123) 



AT 



/ 

Jm 



2 Jm 

The inverse of the one-loop propagator is^^ 



(2.124) 



iA''\x,x') 



(5^(x — x)- 



1 



-9 b 
(2.125) 



Finally, for the CTP-2PI effective action at leading order in the 1/N expansion, we 
obtain 



ihN 



Indet [Gab] 



+ 



ihN 



[ d^xy/^, [ d^x'^,A''\x',x)Gab{x,x') (2.126) 
Jm Jm 



\h^N 



flbcde 



/ d'^xyf^^Gab{x,x)Gcd{x,x) + 0{1). 
Jm 



^^Note that the index 6 is not to be summed in the right-hand side of Eq. (2.125), and the c subscript 

on □ and 7? is a CTP index. 



50 



Applying Eq. (2.82) and taking the limits (p^ = (p- = (p and = g'^ = g^^ , we 
obtain the gap equation for Gab at leading order in the 1/N expansion, 

{G-^r\x,x') = A^\x,x') + '^c''^-'^G,d{x,x)5\x - x')^ + ^ (^) ' ^^"^^^^ 



where 
iA''\x,x') = - 



c«^'[n + m'' + ei?(x)] + ^c"^-^4(x)<^,(x) 



5^{x-x')^=. (2.128) 



-9' 

Similarly, we obtain the mean-field equation for cp at leading order in the 1/A^ expan- 
sion, 

U + m'' + ^R+^^^ + ^G{x, x)^ ^{x) + O (^1^ = 0, (2.129) 

where G{x,x) = Gab{x,x); this definition is the same for all a, 6 € {+,—}, which can 
be seen from Eq. (2.127). To get a consistent set of dynamical equations, we need only 
consider the -|— |- component of Eq. (2.127). It should also be noted that Gab{x,x) is 
formally divergent. Regularization of the coincidence limit of the two-point function 
and the energy-momentum tensor is necessary, and will be carried out in Chapter 3 
below. Multiplying Eq. (2.127) by G and integrating over spacetime, we obtain a 
differential equation for the -|— |- CTP Green function. 



n, + m'' + ^R{x) + ^4>\x) + ^G{x,x)]Gix,x') + o(^^ - '^'^ 



x — x') 



V-9' 

(2.130) 

where appropriate boundary conditions must be specified for to obtain the time- 
ordered propagator as the solution to Eq. (2.130). 

Equations (2.129) and (2.130) are the covariant evolution equations for the mean 
field 4> and the two-point function G++ at leading order in the 1/N expansion. Follow- 
ing Eq. (2.105), we denote the coincidence limit hG{x,x) by {(p"^). With the inclusion 
of the semiclassical gravity field equation (2.84), these equations form a consistent, 
closed set of dynamical equations for the mean field (f>, the time-ordered fiuctuation- 
field Green function G++, and the metric g^i,. 
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The one-loop equations for (p and G can be obtained from the leading-order equa- 
tions by dropping the hG{x,x) term from Eq. (2.130), while leaving the mean-field 
equation (2.129) unchanged. In the Hartree approximation, the gap equation is un- 
changed from Eq. (2.127), and the mean-field equation is obtained from Eq. (2.129) 
by changing h ^ 3h [108]. The principal limitation of the leading-order large- 
approximation is that it neglects the setting-sun diagram which is the lowest-order 
contribution to collisional thermalization of the system [68]. The system, therefore, 
does not thermalize at leading order in the 1/N expansion, and the approximation 
breaks down on a time scale T2 which is on the order of the mean free time for binary 
scattering [118]. 

Let us now use Eq. (2.81) to derive the bare semiclassical Einstein equation for 
the 0{N) theory at leading order in 1/N. This equation contains two parts SS'^/Sg'^ 
and 6T/5g^ . The latter part is related to the bare energy-momentum tensor (T^;y) 
by Eq. (2.85). At leading order in 1/N, {T^j^i,) is given by a sum of "classical" and 
"quantum" parts (distinguished by the latter having an overall factor of h), 

XN 



(2.131) 



where we define the classical part of by 



= N 



(2.132) 



and the quantum part of (T^y) by 

(1 - 20V^V', + (2^ - ^ g^^g'^V^V, - 2^V^V, + 2^g^,gf"'VpV, 

G++(x,x')l+0(l). 



= Nh lim 



(2.133) 
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The above expression for is divergent in four spacetime dimensions, and needs 
to be regularized or renormalized. The energy-momentum tensor in the one-loop 
approximation is obtained by neglecting the 0{h^) terms in Eq. (2.133). It can be 
shown using (2.130) that the energy-momentum tensor at leading order in the 1/N 
expansion is covariantly conserved, up to terms of order 0(1) (next-to- leading-order) . 
The bare semiclassical Einstein equation is then given (in terms of (1)^;/) shown above) 
by Eq. (2.84). 

At this point we formally set N = 1 since we are not including next-to-leading- 
order diagrams in the 1/A^ expansion. This can be envisioned as a simple rescaling of 
the Planck mass by \/iV, since the matter field effective action T is entirely of order 
0{N). We now turn to the issue of renormalization. 

2.5.4 Renormalization 

To renormalize the leading-order, large-A'', CTP effective action in a general curved 
spacetime, one can use dimensional regularization [157], which requires formulating 
effective action in n spacetime dimensions. This necessitates the introduction of a 
length parameter into the classical action, A A/i"^"", in order for the classical 
action to have consistent units. As before, we maintain the restriction g'^ = g'^ = g^^^ , 
and we suppress indices inside functional arguments. 

Making a substitution of the gap equation (2.130) into the leading-order-large- AT, 
2PI effective action (2.126), we obtain 

r[<^,5n = S^[lgn + ^Trln [{0-^ _ 



/ d^xV^c''''^[Gab{x,x)Gcd{x,x)], (2.134) 



in terms of the CTP propagator Gab{x,x'). For convenience we shall rewrite the gap 
equation as 

{G-Y = i + X"'(^)) S{x - x')-^, (2.135) 
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in terms of a four-component "effective mass 
X"''(a;) = (m2 + ei?)c«^ + 



abed r 



.<f>d + hGcdix,x)]. 



(2.136) 



The divergences in the effective action can now be made explicit with the use of a CTP 
generahzation of the heat kernel [105, 157, 158]. Let us define a function K\{x,y,s) 
which satisfies the equation 



ds 



+ / d^z^,c^{G-^r%x,z)K\{z,y-s) = Q, 

JM 



(2.137) 



(where s is a real parameter) with boundary conditions 



K\{x,y;0) = S\S{x-y) 



(2.138) 



at s = [154]. Prom Eqs. (2.137) and (2.134) it follows that K+_ = K'^ = for all 
X, y, and s, and that K\ [K^J) is a functional of (0_) only. The CTP effective 
action can then be expressed as 



(2.139) 



in terms of a functional Psd on M defined by 



+ 



ihN 



, -5 / —K\{x,x]s) 

^ JM Jo 

8 Jm 



xv-g 



f 



dsK'^^{x, x; s) 



1 2 



, (2.140) 



and similarly for P^d. It follows from Eq. (2.137) that K\{x,x; s)[(f>^] is exactly the 
same functional of as K~_{x, x; s)[0_] is of we denote it by K{x, x; s)[0], where 
(/) is a function on M. 

The divergences in the effective action Pgo arise in the small-s part of the integra- 
tions, so that only the first term on the right-hand side of the equation 

r°° ds /"*° ds f°° ds 

/ —K{x,x;s)= —K{x,x;s)+ —K{x,x;s) (2.141) 

Jo S Jo S Jso 
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is divergent. Using the s — > 0"*" asymptotic expansion for K{x,x; s) [154], one has (for 
a scalar field, such as the unbroken symmetry, large- A'" limit of the 0(A'") model) 

oo 

K{x, x; s) ~ (47rs)-t ^ s"*a^(x), (2.142) 

m=0 

where the an{x) are the well-known "HaMiDeW coefficients" made up of scalar invari- 
ants of the spacetime curvature [105, 129, 130]. The divergences then show up as poles 
in l/(n — 4) after the s integrations are performed. They have been evaluated for the 
theory in a general spacetime by many authors (see, e.g., [154, 159, 160]) and in 
the large- iV limit of the O(A^) model [150]. At leading order in the 1/A^ expansion, 
the renormalization of A, ^, m, G, A, b, and c is required, but no field amplitude 
renormalization is required [117,150]. In Chapter 3 below we carry out an explicit 
renormalization of the large- AT dynamics in spatially flat FRW spacetime. 

2.6 Summary 

In this chapter, we have presented a method for deriving causal, coupled equations of 
motion for the mean field and two-point function for an interacting quantum field in 
an arbitrary, classical background spacetime. We derived the equations at two loops at 
leading order in the 1/A^ expansion, and in the latter case showed that renormalization 
counterterms for the "in-out" formulation of the theory are all that is necessary to 
renormalize the effective action. In Chapter 3, we apply these equations to the study 
of inflaton dynamics during the reheating period of inflationary cosmology. The mean 
field and gap equations derived here are also useful, by changing — > —mP, for 
describing the dynamics of symmetry-breaking phase transition [99, 118]. 
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CHAPTER 3 



Parametric particle creation and curved spacetime effects 
3.1 Introduction 

In this Chapter we study the nonperturbative, out-of-equihbrium dynamics of a mini- 
mally coupled scalar 0{N) field theory, with quartic self-interaction, in a spatially flat 
FRW spacetime whose dynamics is given self-consistently by the semiclassical Ein- 
stein equation. The purpose of this study is to understand the preheating period in 
inflationary cosmology, with particular emphasis on the effect of spacetime dynamics 
on the phenomenon of particle production via parametric ampliflcation of quantum 
fluctuations. Of primary interest is obtaining the dynamics of the inflaton (includ- 
ing back reaction from created particles) using rigorous methods of nonequilibrium 
fleld theory in curved spacetime [68, 83]. We have chosen to focus in this Chapter on 
parametric amplification of quantum fluctuations because this phenomenon can be the 
dominant effect in the preheating stage of unbroken symmetry inflationary scenarios, 
among which the chaotic inflation scenarios most directly necessitate (through initial 
conditions) considerations of Planck-scale physics. "New" inflationary scenarios which 
involve a spontaneously broken symmetry often contain additional subtleties (e.g., in- 
frared divergences, spinodal instabilities), and are the subject of ongoing investigation 
[40]. The results of our work are, therefore, particularly relevant to chaotic inflation 
scenarios [161]. The additional interactions which should be included to treat the 
broken-symmetry case are discussed in [80]. 

In Chapter 2 we derived the evolution equations for the mean fleld ($h) (subscript 
H denotes the Heisenberg fleld operator) and mean-squared fluctuations (variance) 
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($1) — (^>h)^ using the closed-time-path (CTP), two-particle-irreducible (2PI) effective 
action [79] in a fully covariant form. Here we use these results for the case of spatially 
fiat FRW spacetime. The quantum state for the field theory (in the case of FRW 
spacetime) consists of a coherent state for the spatially homogeneous field mode, and 
the adiabatic vacuum state for the spatially inhomogeneous modes. At conformal past 
infinity, the spacetime is assumed to be asymptotically de Sitter, and the mean field 
is chosen to be asymptotically constant. 

In this Chapter we study the 0{N) field theory using the 1/N expansion, which 
yields nonperturbative dynamics in the regime of strong mean field. This is particularly 
important for chaotic inflation scenarios [5] , in which the inflaton mean-field amplitude 
can be as large as Mp/3 at the end of the slow-roll period [11, 147]. Treatments of the 
reheating problem which rely on time-dependent perturbation theory do not apply to 
such cases where the inflaton undergoes large-amplitude oscillations, in contrast with 
nonperturbative methods such as large N. 

We now summarize the principal distinctions between our work and previous treat- 
ments of preheating in inflationary cosmology, which were summarized in Section 1.2. 
Our work improves on the group lA methodology by including parametric resonance 
effects. As it is based on first-order, time-dependent perturbation theory, the group 
lA approach cannot correctly describe the inflaton dynamics with large initial mean- 
field amplitude. In addition, our work improves on both the group lA and group IB 
studies by including the effect of back reaction from quantum particle creation on both 
the mean field and the inhomogeneous modes. We are treating the inflaton dynamics 
from first principles, without assuming a phenomenological equation (with a damping 
term T(f) put in by hand) for the mean field. In our approach the damping of the mean 
field is due to back reaction from quantum particle production in the self-consistent 
equations for the mean field and its variance. In contrast, the analytic results of the 
group lA and IB work are based on the assumption of either large-amplitude mean- 
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field oscillations (A^ /2 w?) or harmonic oscillations {w? » cf) /2), and, therefore, 
cannot describe the interesting case of inflaton dynamics in which neither term dom- 
inates the tree-level effective mass, i.e., m ~ . Furthermore, our work improves 
on group lA, group IB, and group 2A in that the closed-time-path effective action 
is computed in curved spacetime without assuming that ^ tq (where tq is the 

period of mean-field oscillations). In our work, the dynamics of the two-point func- 
tion (which reflects quantum particle production) is formulated in curved spacetime 
assuming only that semiclassical gravity is valid, i.e., Mp » H. 

Most significantly, our work improves on all the previous treatments in that it in- 
cludes curved spacetime effects systematically using the coupled, self-consistent semi- 
classical Einstein equation and matter field equations. Among the group 2B studies 
of preheating dynamics, inflaton dynamics has been studied primarily in fixed back- 
ground spacetimes: Minkowski space [99,108], de Sitter space [122,124,162], and in 
radiation-dominated, spatially fiat FRW spacetime [124]. In the present work, the 
spacetime is dynamical, with the renormalized trace of the semiclassical Einstein equa- 
tion governing the dynamics of the scale factor a. This permits quantitative study of 
the transition of the spacetime from the (slow-roll) period of vacuum-dominated ex- 
pansion to the radiation-dominated ("standard") FRW cosmology. In particular, our 
method yields the spacetime dynamics naturally, without making reference to an "ef- 
fective Hubble constant" (which has been used in calculations on a fixed background 
spacetime [162]). 

With additional couplings (see [80]), our method may also be used to study pre- 
heating in "new" inflationary scenarios [4]. In new inflation, the vacuum-dominated 
expansion of the Universe is typically driven by the classical potential energy of the 
mean field as it rolls towards the symmetry-broken ground state. In one of the group 
2B studies (Boyanovsky et al. [162]), a quench-induced phase transition is studied 
with small initial mean-field amplitude, in which the classical terms in the mean-field 
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equation are dominated by spinodal fluctuations. As a result, the mean field in their 
model does not oscillate about the symmetry-broken ground state as is generally ex- 
pected in a new inflation preheating scenario (this point was emphasized in [163]). The 
initial conditions studied in [162] are more appropriate to a study of defect formation 
in a quench-induced phase transition than preheating dynamics in new inflation. 

In addition, the renormalization scheme employed in [162] is not generally covari- 
ant (as can be seen by comparing it with [164]), and covariant conservation of the 
renormalized energy-momentum tensor is put in by hand. The regularization scheme 
employed here is the well-tested adiabatic regularization [26, 100-102], which is simple 
to use and physically intuitive. It also, in the one-loop case, ensures both covariant 
conservation of the regularized energy-momentum tensor and agreement with mani- 
festly covariant regularization procedures such as point splitting [164]. 

A related difference between our approach and that of some of the group 2B stud- 
ies is the choice of vacuum state. The choice of initial conditions for the quantum 
mode functions in most studies of reheating in FRW spacetime [116, 120, 165] has 
been to instantaneously diagonalize the matter-field Hamiltonian at the initial-data 
hypersurface. However, as has been pointed out long ago [130], this method does not 
correspond to the vacuum state which registers the least particle flux on a comoving 
detector. In our work we use the de Sitter-invariant (or Bunch-Davies) vacuum, ob- 
tained via the adiabatic construction; the adiabatic vacuum most closely aligns with 
an intuitive notion of vacuum state in a cosmological spacetime [17]. 

In many of the group 2B studies [108, 124, 162] the large- AT equations for the mean 
field and variance are derived using a factorization method which does not readily 
generalize to next-to-leading order in the 1 /N expansion. In all of the group 2B studies 
of nonpcrturbative inflaton dynamics of which wc arc aware, the equations for the 
mean field and variance are not derived using methods which encompass higher-order 
correlations in the Schwinger-Dyson hierarchy. As found in earlier studies of phase 
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transitions [37,84], this is necessary in order to derive the correct infrared behavior 
of a quantum field in a study of critical phenomena. In the present work, we use 
the result of Chapter 2 in which the CTP two-particle-irreducible (2PI) formalism 
is derived. It has a direct generalization in terms of the n-particle-irreducible (nPI) 
"master" effective action [82]. The master effective action can be used to derive a 
self-consistent truncation of the Schwinger-Dyson equations to arbitrary order in the 
correlation hierarchy [82]. The techniques employed here are, therefore, most readily 
generalized to the study of phase transitions in curved spacetime, where higher-order 
correlation functions can become important [40]. 

In summary, our approach to the inflaton dynamics problem has the following 
advantages: it is nonperturbative and fully covariant; it is based on rigorous meth- 
ods of nonequilibrium field theory in curved spacetime; we use the correct adiabatic 
vacuum construction; and we employ an approximation scheme which can be system- 
atically generalized beyond leading order, within a fully covariant and self-consistent 
theoretical framework. 

Our results are obtained by solving the mean-field and spacetime dynamics self- 
consistently using the coupled matter-field and semiclassical Einstein equations in a 
FRW spacetime, including the effect of back reaction of the variance on the mean field. 
Within the leading-order, large N approximation used here, we find that (using the 
conventional value for the self-coupling, A = 10~^^) for sufficiently large initial mean- 
field amplitude, parametric amplification of quantum fluctuations is not an efficient 
mechanism of energy transfer from the mean field to the inhomogeneous field modes. 
In this case the energy density of the inhomogeneous modes remains negligible in 
comparison to the mean-field energy density for all times. This can be understood 
from the time scales for the competing processes of parametric resonance and cosmic 
expansion. When the time scale for parametric amplification of quantum fluctuations 
Ti is of the same order as (or greater than) the time scale for cosmic expansion H~^, 
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cosmic expansion redshifts the energy density of the inhomogeneous modes faster than 
it increases due to parametric resonance. We find that this occurs when > Mp/300, 
for the model and couphng studied here. 

In many chaotic inflation scenarios, the mean-field amplitude at the end of the 
slow-roll period can be as large as Mp/3 [11, 147]. In light of our result, in such 
models, it is clearly essential to include the effect of spacetime dynamics in order 
to study mean-field dynamics and resonant particle production during reheating. In 
addition, our result indicates that for the case of a minimally coupled inflaton with 
unbroken symmetry, parametric amplification of its own quantum fluctuations is not a 
viable mechanism for reheating, unless the coupling is significantly strengthened (see 
[54, 55]). Parametric amplification of quantum fluctuations may still play a dominant 
role in the reheating of chaotic inflaton models with an inflaton coupled to other fields, 
e.g., a (t?x^ model [87,88,95, 116, 166]. 

For more moderate cosmic expansion, where > 100 ti, parametric amplifi- 
cation of quantum fluctuations is an efficient mechanism of energy transfer to the in- 
homogeneous modes, and the asymptotic effective equation of state is found to agree 
with the prediction of a two-fluid model consisting of the elliptically oscillating mean 
field and relativistic energy density contained in the inhomogeneous mode occupations. 
In a collisionless approximation, the mean field eventually decouples from the mean- 
squared fluctuations (variance) and at late times undergoes asymptotic oscillations 
which arc damped solely by cosmic expansion [147]. For the case when cosmic expan- 
sion is subdominant, 3> ri, the mean-field dynamics and the growth of quantum 
fluctuations are in agreement with results of studies of preheating in Minkowski space 
[99]. In particular, the total adiabatically regularized energy density is found to be 
constant (to within the limits of numerical precision) for the case of ^ oo, in 
agreement with the predictions of field theory in Minkowski space. 

While there has been a large volume of work on understanding the preheating 



61 



period and parametric particle creation, the thermalization of inflationary models has 
not yet been understood from first principles [167]. Because of the absence of a sepa- 
ration of microscopic and macroscopic time scales at the end of the preheating stage, 
the Boltzmann equation is inadequate for studying collisional thermalization in most 
inflationary models [99]. In particular, in the leading-order, large- approximation 
employed here (and in the one- loop approximation which it contains) , this model also 
does not thermalize. However, it still may approach a radiation-dominated effective 
equation of state (as found in [99] for the case of Minkowski space) . Clearly, a first- 
principles analysis of thermalization is necessary. Continuing the early work of kinetic 
field theory [68] , and the recent work on correlation hierarchy [82] , we know that such 
a first-principles analysis should involve at a minimum the full two-loop, two-particle- 
irreducible effective action (or alternatively, next-to-leading order in the large- A?^ ap- 
proximation). Since it represents a rigorous truncation of the full Schwinger-Dyson 
hierarchy, in this sense it is a natural generalization of the coUisionless approxima- 
tions used previously to study reheating. However, the equations derived from it for 
the mean-field and two-point function are nonlocal and hence difficult to solve even 
numerically [117]. This Chapter is, therefore, concerned only with preheating via para- 
metric resonance particle creation. A possible approach to the thermalization problem 
is described in Chapter 5. 

This Chapter is organized as follows. In Sec. 3.2 we present the general theory 
of nonequilibrium dynamics of a scalar field in curved spacctimc, including a sum- 
mary discussion of reheating in inflationary cosmology. In Sec. 3.3 we specialize to the 
case of spatially flat FRW spacetime, and derive the dynamical equations. The ini- 
tial conditions used and the results obtained from numerically solving the dynamical 
equations are described in Sec. 3.4. Discussion and conclusions follow in Sec. 3.5. 
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3.2 Inflaton Dynamics in FRW spacetime 

3.2.1 quantum fields in curved spacetime 

As a simple model of inflation, let us consider a scalar field in semiclassical gravity, 
where the matter field is quantized on a classical, dynamical background spacetime. 

The classical action has the form 

where is the matter field action defined in Eq. (2.64) and 5° is the gravity action 
defined in Eq. (2.65). 

The inflaton field (p is then quantized on the classical background spacetime; we de- 
note the Heisenberg field operator by $h, and the quantum state by Of particular 
importance in a study of inflaton dynamics are the mean field 

4>{x) = {n\^»{x)\n), (3.2) 

the fluctuation field 

(/Ph(x) = $h(x) - 0(x), (3.3) 
and the mean-squared fluctuations, or variance 

{n\^l{x)\n) = {n\^l{x)\n) - {n\^^ix)\nf. (3.4) 

In Chapter 2, a systematic procedure was presented for deriving real and causal evo- 
lution equations for the mean field, two-point function, and the metric tensor in semi- 
classical gravity. Assuming a globally hyperbolic spacetime, one can evolve the coupled 
evolution equations forward from initial data specified at an initial Cauchy hypersur- 
face. The evolution equations follow from functional differentiation (and subsequent 
field identifications) of the closed-time-path (CTP) two-particle-irreducible (2PI) effec- 
tive action, r[(f), G,g'^'']. The CTP-2PI effective action is a functional of the mean field 
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0, two-point function G, and metric tensor g^^'^ , which now carry not only spacetime 
labels but also time branch labels, which have an index set {+, — }. The evolution equa- 
tions for (f>, (<^|), and g^i, then follow from Eqs. (2.82), (2.83), and (2.81), respectively. 
The variance ((^|) is related to the coincidence limit of any of the four components 
of the CTP two-point function through Eq. (2.48). The energy-momentum tensor 
(T^jy) is defined by Eq. (2.85), and it is this quantum expectation value which (after 
renormalization) enters as the source of the semiclassical Einstein field equation (2.84). 
Eqs. (2.81)-(2.83) constitute a set of coupled, nonlocal, nonlinear equations for the 
mean field, two-point function, and metric tensor. The renormalizcd versions are what 
enter into the description of inflaton dynamics. The CTP-2PI effective action can be 
computed using diagrammatic methods described in Chapter 2, where a covariant ex- 
pression for T[(l),G,g^''] was computed in a general curved spacetime (truncated at 
two loops). 

3.2.2 Inflaton dynamics in FRW spacetime 

We now consider a spatially flat Priedmann-Robertson- Walker (FRW) spacetime, which 
is spatially homogeneous, isotropic, and conformally flat. Its line element can be writ- 
ten in the form 



ds'^ = a{rif 



i=l 



(3.5) 



where a is the scale factor, {i € {1,2,3}) are the physical position coordinates on 
the spatial hypersurfaces of constant conformal time rj (related to the cosmic time t 
by = / dt/a). The Hubble parameter, which measures the rate of cosmic expansion, 
is 

H{r,) = ^, (3.6) 

where the over-dot denotes differentiation with respect to cosmic time t. Given our 
choice of sign convention and metric signature, the Ricci tensor in the FRW coordinates 
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is given by 



i?oo = 3 
Rij = — 



a" {a'f 
a 

a 



(3.7) 
(3.8) 



where the prime denotes differentiation with respeet to 77, and Rqq is the component 
of the Ricci tensor proportional to dr] ® drj. The scalar curvature is 

6a" 



and the Einstein tensor is 

Goo 

Gij 

Finally, the volume form on M is 



R 



a a? 



dij. 



Cm = a^{dr] A dx^ A dx^ A dx^). 



(3.9) 

(3.10) 
(3.11) 

(3.12) 



The higher-order (e.g., R^) geometric terms in the geometrodynamical field equation 
are not shown because the renormalized constants b and c are set to zero in Sec. 3.3.4. 
In restricting the spacetime to be a spatially flat FRW, we are reducing the number 

of degrees of freedom in the metric: 



(3.13) 



This reduction should not be carried out in the 2PI generating functional r[<^, G,g*^"\, 
but only in the equations of motion (2.81)-(2.83). This is because functional differen- 
tiation of r[(/), G, a~'^r|^^^] with respect to the scale factor a gives only the trace of the 
energy-momentum tensor, a~'^ri^^ [T^j^^,) , and not the additional constraint equation 
which the initial data must satisfy. 

The spatial homogeneity and isotropy of FRW spacetime permits only two alge- 
braically independent components of the energy-momentum tensor, which in the FRW 
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coordinates of Eq. (3.5) arc given by (loo) and (Tu); all other components are zero. 
These must be functions of rj only (due to spatial homogeneity). For the purpose of 
numerically solving the semiclassical Einstein equation, it is convenient to work with 
the trace 

r = g^''{T^,) = a-V{T^.u), (3.14) 

instead of {Ta). The trace T enters into the dynamical equation for 0(77), and (Too) 
enters into the constraint equation. 

Another consequence of the spatial symmetries of FRW spacetime is the restriction 
on the generality with which we may specify initial data for dynamical evolution. Let 
us choose to specify initial data on a Cauchy hypersurface of constant conformal 
time 770 • In the Heisenbcrg picture,^ for consistency with spatial homogeneity, the 
quantum state |0) must satisfy 

{n\^nir]o,x)\n) = ^{r]o), (3.15) 

{n\^'^{rjo,x)\n) = 4>'{no), (3.i6) 

for all X G M^, where $h is the Heisenberg field operator for the scalar field. The 
values of 4>{vo) and 4> (vo) constitute initial data for the mean field. In addition, the 
quantum state must satisfy 

{n\ip^{rjo,x)ip^{rio,x')\n) = F(r/o, \x-x'\), (3.17) 

— {^\ip^^{'n,x)ip^{'r],x )\n) = F'{r)o,\x - x\), (3.18) 

in terms of an equal-time correlation function F{riQ, \x — x'\) which is invariant under 

simultaneous translations and rotations of x and x' . As defined in Eq. (3.3), ip^ 

denotes the Heisenberg field operator for the fluctuation field. The spatial Fourier 

transform of F is related to the power spectrum of quantum fluctuations at rjo for the 

'^As discussed in Sec. 3.2.3, for our purposes it is sufficient to consider only the case of a pure state. 
The analysis can, however, be easily extended to encompass a mixed state with density matrix p. 
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quantum state Alternatively, we may say that F{riQ,r) and F'{r]Q,r) give initial 
data for the evolution of the two-point function via the gap equation (2.83). The 
symmetry conditions (3.15), (3.16), (3.17), (3.18), along with the spatial symmetries 
of the classical action in FRW spacetime, guarantee that the mean field and two-point 
function satisfy spatial homogeneity and isotropy for all time, i.e., 

(<1>h(x)) = 4>{v), (3.19) 
G++{x,x') = G++{r],r]',\x-x'\), (3.20) 

for all X £ M. The conditions (3.19), (3.20) permit a formal solution of the gap equa- 
tion (2.83) for G++ in terms of homogeneous mode functions, via a Fourier transform 
in comoving momentum k, as shown in Sec. 3.3.2. By rotational invariance, the Fourier 
transform depends only on the magnitude k = ■ k. Of course, the quantum state 

is not uniquely defined by the spatial symmetries; a unique choice of the initial 
conditions for <j) and Gab at S^,, is (in the Gaussian wave-functional approximation) 
equivalent to choosing The choice of quantum state depends on the physics of 
the problem we wish to study. 

As a consequence of covariant conservation of the energy-momentum tensor 

V'^(r^.) = 0, (3.21) 
the functions (Too(77)) and {Ta^r])) satisfy 

which comes from taking the z/ = component of Eq. (3.21). In analogy with the 
continuity relation for a classical perfect fluid in FRW spacetime, 

^(«V) = -P^M, (3.23) 
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we may define the energy density p and pressure by 



(3.24) 



(3.25) 



However, the quantity p should not be interpreted as the true hydrodynamic pressure 
until a perfect-fluid condition is shown to exist; otherwise, bulk viscosity corrections 
can enter into Eq. (3.25) [168]. The effective equation of state is defined as a time aver- 
age (over the time scale ri for the matter field dynamics, to be discussed in Sec. 3.2.4) 
of the ratio p/p, 



The effective equation of state 7 (where the bar denotes a time average) is an important 
quantity in differentiating between the various stages of inflationary cosmology. 

Several solutions to the semiclassical Einstein equation (2.87) for idealized equa- 
tions of state are of particular interest in cosmology. The effective equation of state 
7 = — 1 (eternally "vacuum dominated") leads to a solution 0(77) = ~l/{Hri), for 
—00 < ?7 < 0, where H = -ySvrGp/S and p is a constant. This solution corresponds 
to the "steady-state" coordinatization covering one-half of the de Sitter manifold [17]. 
The effective equation of state 7 = corresponds to nonrelativistic matter, in which 
case the scale factor conformal-time dependence is a a ry^. The effective equation of 
state 7 = 1/3 corresponds to relativistic matter, and its scale factor conformal-time 
dependence is a oc rj. 

3.2.3 Initial conditions for post-inflation dynamics 

In most realizations of inflationary cosmology, the Universe evolves through a period 
in which a dominant portion of the energy density p comes from a quantum field 
$H, the inflaton field, whose effective equation of state [defined as in Eq. (3.26)] is 

7 ~ — 1. In chaotic inflation, this condition is due to the fact that the inflaton fleld 




(3.26) 
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is in a quantum state in which the Heisenberg field operator <J>h acquires a large 
(approximately spatially homogeneous) expectation value, defined by Eq. (3.19). A 
requirement for chaotic inflation is that the potential energy F(0) of the expectation 
value (f) dominates over both the spatial gradient energy [coming from ((Vyjjj)^)] ^^^d 
kinetic energy for the inflaton field, and the energy density of all other quantum fields 
coupled to the inflaton. The potential energy V{^) gives a contribution to the energy- 
momentum tensor satisfying precisely 7 = —1. During inflation, the scale factor 
grows by a factor of approximately exp{HAt), where At is the interval of inflation 
in cosmic time, typically larger than 60H^^. While the Universe is inflating, the 
expectation value ($h) is slowly rolling toward the true minimum of the effective 
potential. (In reality, the situation is much more complicated than this. The effective 
potential is an inadequate tool for studying out-of-equilibrium mean-field dynamics [15, 
169].) Assuming the Universe was in local thermal equilibrium prior to inflation, the 
temperature during inflation decreases in proportion to 1/a. The energy density of any 
relativistic (nonrelativistic) flelds coupled to the inflaton is proportional to 1/a^ (1/a^). 
The contribution to the quantum energy density from spatial gradients of fluctuations 
about the inflaton fleld is proportional to (see Sec. 3.3 below). Most importantly, 
any inhomogeneous modes 6(1)^. of the mean field which might exist at the onset of 
inflation are redshifted. The physical momentum of a quantum mode, /Cphys = k/a, 
decreases as 1/a relative to the comoving momentum k. The quantum state of any 
fleld coupled to the inflaton at the end of inflation is, therefore, approximately given by 
the vacuum state. The inflaton fleld is well approximated by a spatially homogeneous 
mean field, with vacuum fluctuations around the mean-field configuration. The mean 
field can be thought of as representing the coherent oscillations of a condensate of 
zero-momentum inflaton particles. 

Let us consider the case of inflation driven by a single self-interacting scalar field 
(f) (with unbroken symmetry) in spatially flat FRW spacetime. The above arguments 
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imply that one can model post-inflationary physics with a quantum state which at 
rjo corresponds to a coherent state for the field operator $h [in which (ri|$H(a;)|fi) = 
(^(77)], and the fluctuation field tp^^ is very nearly in the vacuum state. ^ Then for 
rj < rjo, (Too) is dominated by the classical energy density of the mean field (f). The 00 
component of the Einstein equation then yields 



a 



af_ ISttGpc 

~2 



(3.27) 



where pc is the classical energy density of the mean field, defined by 

Pc = ^i^'? + V{4>). (3.28) 
The mean field 4> satisfies the classical equation 

<^" + — ^' + aV(<^) = 0, (3.29) 



where V{(p) denotes the classical potential. For the A^'* theory, the potential is [from 
the Minkowski-space limit of Eq. (3.1)] 

V{^) = \m^' + (3.30) 

The assumption that the Universe is inflating (i.e., 7 ~ —1) for 77 < 770 requires that 
the energy density pc be potential dominated, 

F(<^)» 2^(0V, (3.31) 

and that the mean field satisfies the slow-roll condition, 

7// 2a' -1 

(p < — (p . (3.32) 
a 



^Though this is a pure quantum state, the methods employed in this study can be used to treat a 
quantum field theory in a mixed state (for example, a system initially in thermal equilibrium with a 
heat bath). 
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Given Eqs. (3.31) and (3.32), an approximate "0th adiabatic order" solution to the 
Einstein equation can be obtained [normahzed to a{r]o) = 1], 

a{ri) ~ — -^7 r, (3-33) 

where H is & slowly varying function of ij, given by 



H(V) = (3.34) 

Prom Eq. (3.34), we can evaluate the expansion rate nonadiabaticity parameter 
[41] for r) < r)o using Eq. (3.34). During slow-roll it follows from conditions (3.31) and 
(3.32) that 

Oh = = , = < 1. (3.35) 

The solution (3.33) for a(ry) is exact in the limit of constant H (corresponding to a 
constant (p at the tree level). For simplicity, let us assume that (p goes to a constant 
value > Mp in the asymptotic past, 77 — — oo. The spacetime is then asymptoti- 
cally de Sitter, with the scale factor having an asymptotic cosmic-time dependence 
a{t) ~ ex.p{Ht). Because the enormous cosmic expansion during the slow-roll period 
redshifts away all nonvacuum energy in the Universe, it is reasonable to assume that 
the quantum state would register no particles for a comoving detector coupled to 
the fluctuation field (p at conformal-past infinity; i.e., that the fluctuation field ip is in 
the vacuum state at rj ^ —00. This would mean that a ~ l/{Hr]) at r] ^ —00. This 
spacetime is not asymptotically static in the past, but it is conformally static with a 
conformal factor whose nonadiabaticity parameter vanishes at conformal-past infinity. 
Therefore, the best approximation to a "no-particle" state for a comoving detector in 
the asymptotic past is given by the adiabatic vacuum [100] constructed via matching 
at ?7 — > — 00. All higher-order adiabatic vacua will in this case agree at conformal past 
infinity. 

To construct the nth-order adiabatic vacuum matched at an equal-time hyper- 
surface S,,^, one first exactly solves the conformal-mode function equation for the 
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quantum field [see Eq. (3.51) below]. Since the mode-function equation is second or- 
der, each k mode has two independent solutions, which can be represented as Uk and 
Uf.. A particular solution consists of a linear combination of Uk and u^. The adiabatic 
vacuum is constructed by choosing (for each k) a linear combination which smoothly 
matches the nth-order positive frequency WKB mode function at S^^. The resulting 
orthonormal basis of mode functions is then used to expand the Heisenberg field oper- 
ator ^nix) in terms of and a|.. The vacuum state is defined by afe|vac) = for all k, 
which can be shown to correspond (in the rfm — — oo limit) to the de Sitter-invariant, 
adiabatic (Bunch-Davies) vacuum. 

3.2.4 Post-inflation preheating 

Inflation ends when the mean field has rolled down to the point where condition 
(3.32) ceases to be valid, which we assume occurs at conformal time The inflaton 
mean field then begins to oscillate about the true minimum of the effective potential, 
leading to a change in the effective equation of state. A harmonically oscillating 
scalar mean field (m^ » A(/> /6) has an effective equation of state 7 = 0, and a 
scalar inflaton undergoing extreme large-amplitude oscillations {X(j) /6 ^ m^) has an 
effective equation of state 7 = 1/3 [10]. In realistic models, the inflaton field is coupled 
to various lighter fields consisting of fermions and/or bosons. These couplings, as 
well as the inflaton's self-coupling, provide mechanisms for damping of the mean-field 
oscillations via back reaction from quantum particle production, and energy transfer 
to the lighter fields and the infiaton's inhomogeneous modes. 

Let us consider the scalar field theory with unbroken symmetry in Minkowski 
space [with classical action given by the Minkowski-space limit of Eq. (3.1)], and 
suppose that the mean field <p oscillates about the stable equilibrium configuration 
= with initial amplitude 0o- For the moment we are neglecting the effect of 
spacetime dynamics, i.e., assuming a{rf) = 1. The time scale for oscillations of the 
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mean field is given by [99] 



where / and k are defined by 



AKjk) 



f = J^-^ (3-37) 
V 6 m 

k = , ^ (3.38) 

and is the complete elliptic integral of the first kind [170]. For harmonic oscilla- 

tions where 

< rm?, (3.39) 

time-dependent perturbation theory was used in the group lA studies (see Sec. 3.1) to 
compute the damping rate V for the mean field in the model. At lowest order in A, 
the damping rate for the mean field (j) corresponds to the rate for four zero-momentum, 
free-field excitations of the inflaton to decay into a (fluctuation field) particle pair, 
due to the self-coupling [95], 

r^^O(l)^^, (3.40) 
47rm 

with vacuum initial state for ip. The symbol 0(1) denotes a constant of order unity. In 
addition to the assumption (3.39), there is another crucial assumption in the deriva- 
tion of Eq. (3.40), namely, that the dominant contribution to the decay rate is given 
by the lowest order, jvac) — ^ jig, l_g) process, where the occupation numbers are 
for the fluctuation fleld ^p. It can be shown [22, 171] that for this (bosonic) case, the 
perturbative decay rate into the k momentum shell for the fluctuation field ip is en- 
hanced by l-|-2n when the occupation of the k shell is n. This is a stimulated emission 
effect due to Bose statistics.^ The use of Eq. (3.40) to estimate the damping rate thus 

^In contrast with the case with Bose fields, the use of time-dependent perturbation theory to study 
inflaton decay into fermions via a Yukawa coupling does not require the condition <C 1, because 
of the Pauli exclusion principle [87, 94, 95]. It is still necessary, however, to assume weak coupling (or 
small mean-field amplitude) in order to use perturbation theory [94,95]. 
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implicitly assumes that for all k, the fluctuation field occupation numbers are small, 
i.e., n-j <C 1. This is because time-dependent perturbation theory in terms of the 
interaction corresponds to an expansion of the field theory around the vacuum 
configuration. Equivalently, it corresponds to an amplitude expansion (in powers of 
the "classical field" (j)±) of the IPX closcd-time-path effective action r[^_|_,(^_], which 
is defined in Eq. (2.11) in Ref. [80]. When A^q sufficiently large at 770, or on a time 
scale for rij^ to grow to order unity, the perturbative expansion in A breaks down. 

In many inflationary scenarios, condition (3.39) does not hold at r/o. A correct 
analysis of the dynamics of the inflaton field must, therefore, be nonperturbativc, if 
the inflaton is self- interacting and/or coupled to Bose flelds. Again, of interest in 
"preheating" is the time scale for damping of the mean field due to back reaction 
from particle production into the inhomogeneous modes of the fluctuation fleld. This 
quantum particle production is known to occur by parametric amplification of quan- 
tum vacuum fluctuations, for the zero-temperature, unbroken symmetry system under 
study here. Boyanovsky et al. [99] have obtained an approximate analytic expres- 
sion (in Minkowski space) for the time scale ri for the variance ((p^) to grow to the 
point where X{ip'^)/2 is of the same order of magnitude as the tree-level effective mass 



+ Xcj) /6, 



r.^^mf ''IZ'/^' V (3.41) 

The function B(f) is of order unity, and in terms of the asymptotic value of / at 
ry — 00, B[f{rf — 00)] ~ 0.285 953. Their result is valid in flat space and based 
on a solution of the one-loop dynamics which neglects the back reaction of particle 
production on the mode functions. The essential feature of the time scale ri is that 
it depends on the ln(A~^). As a consequence of the analytic solution to the classical 
mean-field equation and the estimate for ri, it is possible to estimate (for the case of 
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Minkowski space) the effective equation of state 7c for the mean field [99] 




.A? [i-iA?] + 1(1 + ,/;?) 1 



E(k) 



(3.42) 



where E(k) is the complete elliptic integral of the second kind [170], pc is the pressure 
of the mean field, and pc is the energy density of the mean field, defined in Eq. (3.28). 
The late-time effective equation of state can be studied using an idealized two-fluid 
model consisting of the classical mean-field oscillations % and a relativistic component 
corresponding to the energy density in the quantum modes pq [defined in Eq. (3.61) 
below] . 

The physical processes discussed above neglect coUisional scattering of excitations 
of the inhomogeneous modes due to the self-interaction, for example, binary scat- 
tering. These scattering processes ultimately lead to thermalization of the system. A 
quantitative understanding of the time scales for such processes in the nonperturbative 
regime studied here within a rigorous field-theoretic framework is at present lacking. A 
perturbative treatment of collisional thermalization of the system using the Boltzmann 
equation assumes a separation of time scales for coUisionless processes (ri) and ther- 
malization. However, due to the nonperturbatively large occupation numbers which 
arise in the resonance band of the inhomogeneous field modes on the time scale ri, 
such a naive approach would predict that the time scale for thermalization is on the 
order of (or earlier than) the preheating time scale ri. A nonperturbative approach to 
studying the collisional thermalization of the system is, therefore, required. However, 
within the 1/N expansion (to be discussed in Sec. 3.3.1), the collisional scattering pro- 
cesses are subleading order in 1/N, and thus the separation of time scales is assured 
within this controlled expansion [99]. Let us denote the time scale for scattering by 

In typical inflationary scenarios, the self-coupling A of the inflaton is very weak 
[10], in the range lO^-^^-lO"-'^^ (see, however, [54,55]). In our numerical work, / is 



75 



initially unity, in which case the three time scales tq, ti and T2 separate dramatically, 

(3.43) 

T2/T1 ~ 0{N). (3.44) 

The period leading up to ri is called preheating, because (i) the energy transfer from 
the mean field is entirely nonequilibrium in origin, and (ii) the occupation numbers 
of the fluctuation field are extremely nonthermal. In this regime, since T2 ^ n, 
collisional effects can be neglected. In a collisionless approximation, the damping of 
the mean field is due to energy transfer into the inhomogeneous quantum modes, a 
process similar to Landau damping in plasma physics [117]. 

So far in this section we have not included the effect of spacetime dynamics on 
the particle production and back reaction processes. Cosmic expansion introduces an 
additional time scale H~^, where H is the Hubble parameter defined in Eq. (3.34). 
In typical chaotic inflation scenarios, the initial infiaton amplitude can be as large as 
Mp/3, leading to 

at the onset of reheating. In this case, <^ ri when A is very small. Clearly, 

for sufficiently large initial infiaton amplitude, it is necessary to include the effect of 
spacetime dynamics in a systematic study of preheating dynamics of the infiaton field. 

3.3 O(A^) inflaton dynamics in FRW spacetime 

In this section, we study the nonequilibrium dynamics of a quartically self-interacting, 
minimally coupled, O(A^) field theory (with unbroken symmetry) in spatially flat FRW 
spacetime. We use the covariant evolution equations derived in [80] , in order to study 
the dynamics of the mean field, variance, and the spacetime, at leading order in the 
1/7V expansion. 



n/ro ~ O In 



27r 
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3.3.1 The O(A^) model in the expansion 

The classical action for the unbroken symmetry 0(iV) model in a general curved 
spacetime is 



(3.46) 



where the 0(Ar) inner product is defined by^ 

i,4=c^^4P8ij. (3.47) 

As in Eq. (3.1), A is a coupling constant with dimensions of l/?i, and ^ is the dimen- 
sionless coupling to gravity (and is necessary in order for the quantized theory to be 
renormalizable) . 

In [80], the covariant mean-field equation, gap equation, and geometrodynamical 
field equation were computed for this model at leading order in the XjN expansion. 
The evolution equations follow from Eqs. (2.81)-(2.83), with the 2PI, CTP effective 
action truncated at leading order in the XjN expansion. At leading order in the XjN 
expansion, we need only keep track of one component of the CTP two-point function 
Gab{x, x'); we choose G+_|_(x, x'), which is the Green function with Feynman boundary 
conditions. The covariant gap equation for G++ at leading order in the 1/N expansion 
is 

[D^ + m'^ + ^R(x) + ^4>^(x) + ^G(x,x)] G++(x,x') = S'^ix - x')^l=, (3.48) 
V 2 2/ ^-g' 

plus terms of 0{1/N). The covariant 5 function is defined in Ref. [17]. The mean-field 
equation is, at leading order in 1/A/", 

U + m'^ + iR+^^^ + ^G{x, xU ^{x) = 0. (3.49) 



*In our index notation, the latin letters i,j,k,l,m,n are used to designate 0(iV) indices (with 
index set {1, . . . , A'^}), while the latin letters a,b,c,d,e, f are used below to designate CTP indices 
(with index set {+, —}). 
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Recall that all four components of Gab{x,x') arc the same in the coincidence limit, 
which we are denoting by G{x,x). The coincidence limit is divergent in four 

spacetime dimensions, and the regularization method is described in Sec. 2.5.4 below. 
The geometrodynamical field equation is given by Eq. (2.84). in terms of the (unrenor- 
malized) energy-momentum tensor computed at leading order in the expansion, 
which is shown in Eqs. (5.37) and (5.38) in Ref. [80]. 

3.3.2 Restriction to FRW spacetime 

Let us now specialize to the spatially flat FRW universe, with initial conditions ap- 
propriate to post-inflation dynamics of the inflaton field. As discussed in Sec. 3.2.3, 
initial Cauchy data for <j), and a are specified on a spacelike hypersurface S^,, 

(at conformal time 770). The spatial symmetries of and for a quantum state |0) 
consistent with a spatially homogeneous and isotropic cosmology are given in (3.19-b). 
As a consequence of these symmetries, both the mean field (j) and variance {^p\) are 
spatially homogeneous, i.e., functions of conformal time only. 

Eq. (3.48) for G++ in spatially flat FRW spacetime has the formal solution 

G++{x,x')=aiv)-'airj')-' J ^e^^^^"-") [G(V - 7?)%(r?)*^ifc(V) 

+ eiri-v')uk(.V>k{vr], (3.50) 

in terms of conformal-mode functions which satisfy a harmonic oscillator equation 
with conformal-time-dependent effective frequency, 

{£2+^Uv))u, = 0. (3.51) 

The fact that Uk (r?) depends only on rj and k (where k is comoving momentum) implies 
that G++ is invariant under simultaneous spatial translations and rotations of x and 
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x' . The effective frequency ^Ikiv) appearing in Eq. (3.51) is defined by 



(3.52) 



(3.53) 



(3.54) 



Initial conditions for the positive frequency conformal mode functions Ukiv) must 
be specified (for all k) at r^o- A choice of initial conditions corresponds to a choice 
of quantum state for the fluctuation field ip^^; initial conditions are discussed in 
Sec. 3.4.1 below. The (bare) variance {(p^) has a simple representation in terms of the 
conformal-mode functions: 



It should be noted that this expression is divergent, in consequence of our having 
computed the variance in terms of the bare (unrenormalized) constants of the theory. 
In terms of a physical upper momentum cutoff K, the variance G{x, x) diverges like K^. 
Even after removing the 0{K'^) divergence, there remains a logarithmic dependence 
on K which must be regularized. In addition, the mode functions Uk depend on ((^|) 
through Eq. (3.54). The leading-order, large- A/", mean-field equation in spatially fiat 
FRW spacetime becomes 



where the time-dependent bare effective mass M{r]) is given by Eq. (3.54). For sim- 
plicity of notation, we will henceforth write M instead of M{r]), and similarly for 



Finally, we can express the bare energy-momentum tensor in terms of the conformal- 
mode functions Uk{v)- As discussed in Sec. 3.2.2, it is convenient to work with the 
00 component and the trace of the energy-momentum tensor. The components of the 




(3.55) 



(3.56) 



a 



m{r]). 
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classical part of the energy-momentum tensor arc spatially homogeneous, and given 

by 



1 



3e 



2a' 



1 2 / 2,^1^, 3e(a' 



a 

r^ir,) = \ii6^-i)i4>'f + 6^(4>" + —^' 



/^2 



2a4 



+ 2 ( + ^(^% |i? 



(3.57) 

4^ ■ 2-- ■ ^'-''^ 
The quantum energy-momentum tensor components are also spatially homogeneous. 
We find for the 00 component, 

^'o(^) = i / (£f f + ^ ^ " '^^^) 



+ (6^-i)-[(4)*ufe + 4u^ 



, (3.59) 



and for the trace, 



(6^ - !)<; 141^ - {k' + a'M')\uk\' - - {u'kTuk + « 



+ 



{a') 



l\2 



(3.60) 



It can be shown by asymptotic analysis that, in terms of a physical upper momentum 
cutoff K, the bare is quartically divergent, i.e., 0{K^), and that (for minimal 
coupling) T*^ is quadratically divergent. In addition, the components of the bare 
energy-momentum tensor contain the effective mass M^, which contains the divergent 
variance The energy density pg of quantum modes of the (p field is defined in 

terms of by 

Pq = ^T^o - ^{flf- (3-61) 
We shall also refer to Pq as the energy density of the "fiuctuation field." 



3.3.3 Renormalization of the dynamical equations 

The variance ((^|) and quantum energy-momentum tensor components Tqq and T'^ are 
divergent in four spacetime dimensions, and must be regularized within the context 
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of a systematic, covariant renormalization procedure. In the "in-out" formulation of 
quantum field theory, renormalization may be carried out via addition of countert- 
erms to the effective action, which amounts to renormalization of the constants in 
the classical action [172]. The closed-time-path formulation of the effective dynamics 
is renormalizable provided the theory is rcnormalizable in the "in-out" formulation 
[31,67], as is the case with the 0{N) field theory in curved spacetime [80, 150, 154]. 
For our purposes it is convenient (in this model) to carry out renormalization in the 
leading-order, large-iV, evolution equations, rather than in the CTP effective action 
[66]. 

In this study we employ the adiabatic regularization method of Parker, Pulling, 
and Hu [100, 101]. The idea is to define an adiabatic approximation to the conformal 
mode function, and then to construct a regulator for the integrands of the bare energy- 
momentum tensor and variance from the adiabatic mode functions [130]. Renormaliza- 
tion occurs when we define the renormalized variance and energy-momentum tensor to 
be the difference between the bare expressions and the regulators and simultaneously 
replace the bare quantities m, A, G, b, c. A, and ^ by their renormalized counterparts. 
The equivalence of this procedure to other manifestly covariant methods (such as di- 
mensional continuation) is well established [173]. We implement renormalization as 
a two-step process: Pirst, we adiabatically regularize the variance and renormalize 
^, m, and A; Second, we adiabatically regularize the energy-momentum tensor and 
renormalize the scmiclassical geometrodynamical field equation. 

We define the adiabatic order of a conformal mode function as follows: let i^kiv) ~^ 
Uk{v/T), where T is introduced as a time scale which is formally taken to be unity 
at the end of the calculation. Then the adiabatic order of an expression involving 
derivatives of fl^ is simply the inverse power of T, of the leading-order term in an 
asymptotic expansion about T ^ oo. However, in order for the adiabatically reg- 
ulated energy-momentum tensor for an interacting scalar field theory to agree with 
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the renormalized energy-momentum tensor obtained by manifestly covariant methods 
(e.g., covariant point splitting [164]), it is necessary to define the adiabatic order of 
expressions involving A and derivatives with respect to 77, such as A(^ )", as the sum of 
the exponent of (f) and the number of conformal time differentiations [160] . Therefore, 
M^h)" is considered fourth adiabatic order, as is X{cl) )" ■ 

Having defined adiabatic order, we now construct the adiabatic mode functions. 
It is well known that the WKB Ansatz 

Ukiv) = ^7f^ f dHW{ri')^ (3.62) 

turns the harmonic oscillator equation (with time-dependent frequency fife) into a 
nonlinear differential equation for W, 

^^"^^ = ^^^^^ + - m^- ^^-^^^ 

Starting with the lowest-order Ansatz W^^\r]) = Q,k{jl)i can iterate this equation; 
the nth-order iteration yields the nth-order WKB approximation for u^. For the free 
field theory, the nth-order WKB approximation gives an expression for uj^ which is of 
adiabatic order 2n. In the interacting case, the above definition of adiabatic order calls 
for removing terms such as A(^ )"" at 4th adiabatic order. Thus we have a method 
of deriving expressions for T^, T*^, and (<^|) at fourth, fourth, and second adiabatic 
orders, respectively. One then sets T = 1 in the truncated expression. We can thus 
obtain a fourth-order adiabatic approximation to the quantum energy-momentum ten- 
sor (T|2y)ad4) and a second-order adiabatic approximation to the variance {(p^ad2- By 
subtracting {T^i,)^4, from the divergent T^j^ and (<^|)ad2 from the divergent (y'l), finite 
expressions for the renormalized energy-momentum tensor and variance are obtained. 

First we regularize the variance ((/?|), and carry out a renormalization of A, m, and 
^. In the leading-order, large- A'^ approximation, no terms appear in the mode-function 
equation (3.51) which would permit addition of counterterms; therefore, $7^ must be 
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finite [174]. The effective frequency 0,k which appears in Eq. (3.52) is the "bare" 
effective frequency, which we denote by {'^h)B- In conjunction with the adiabatic 
regularization procedure, we fix the renormahzation scheme by demanding equivalence 
of the bare and renormahzed effective mass [117], 

i^Dn = {^Ih, (3.64) 

where the "R" subscripted quantities are renormahzed. Using Eqs. (3.52) and (3.53), 
we have 

UR + Ml = ^^R + Mi, (3.65) 
where M| is defined in Eq. (3.54), 

M| = m| + ^<^%^(<^|)e, (3.66) 

and is defined similarly, 

Ml = ml + ^^' + ^{<plU. (3.67) 

Now, Ab, tub, and suce the bare constants of the theory which appeared (without 
B's) in the classical action (3.46). The renormahzed quantities in Eq. (3.52) are defined 
below. The bare ((^|)b is a conformal-time-dependent function defined by Eq. (3.55), 

{'Pliri))s = ^\Mn)\', (3.68) 

where the conformal-mode functions Uk{il) obey Eq. (3.51). Now we demand that the 
renormahzed theory be minimally coupled, i.e., we set = 0. Because of Eq. (3.64), 
we can formally use in computing the adiabatic regulator for the variance (¥'|)b- 

Computing the asymptotic series (in l/T) of the quantity \uk{v)\^ to 0{1/T'^), where 
^\{r]/T) is the effective frequency, we obtain (after setting T = 1) 



h f d^k 



2C J (27r)2 



1 {C'f-2CC" MlC" 5M^(C') 



+ 



(3.69) 
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in terms of an auxiliary function 



(3.70) 



In Eq. (3.69) the symbol Uk is defined as follows 



(3.71) 



In the adiabatic prescription, the renormalized variance (¥'h)r appearing in Eq. (3.67) 



is defined by 



(3.72) 



where the first term on the right-hand side is given by Eq. (3.68), and the second term 
on the right-hand side is given by Eq. (3.69). Everything on the right hand side can 

be expressed in terms of renormalized quantities, so this procedure is well defined. 
Written out explicitly, the renormalized variance satisfies the equation 



1 

c 



(2^ 



2uk 



1 {C'f-2CC" MlC" 5M^(C" 



i2n 



+ 



(3.73) 



One can use the WKB approximation for Uk{rj) to compute the asymptotic series 
for the integrand in Eq. (3.73) in the limit k — ^ oo, and show that the integral is 
convergent. Since M| is contained in the integrand above, Eq. (3.73) leads to an 
integral equation for the renormalized effective mass Mr, 



Mi 



d^k 



2C J (27r) 



2Cbk 



2CC" 



+ 



MjC" 



5M|(C- 



(3.74) 



Eqs. (3.73) and (3.65) together define Ar and m^. All physical quantities should now 
be expressed in terms of the renormalized parameters and Ar of the theory. The 
renormalized mean-field equation now reads 



2a' 7' 
+ — <j) + a'M, 
a 



2 ;\/r2^ 



0, 



(3.75) 
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where is given by Eq. (3.74), and the mode functions in Eq. (3.74) obey the 
homogeneous equation, 

0. (3.76) 



dr] 



The initial conditions for the conformal-mode functions at r]o are discussed in Sec. 3.4.1 
below. 

Having obtained a renormalized mean-field equation, we now turn our attention to 
regularizing the quantum energy- momentum tensor. As a consequence of Eq. (3.64), 
we can substitute M and ^ in the equations for the components of the 

quantum energy-momentum tensor, Eqs. (3.59, 3.60). Since we wish to study the 
minimal coupling case, we set = 0. To avoid confusion we denote the bare energy- 
momentum tensor components (3.59) and (3.60) by {Tqq)b and {'T'^)b, respectively. Let 
us also relabel the bare constants b, c, G, and A appearing in the bare semiclassical 
Einstein equation (2.84) as 6b, Cb, Gb, and Ab. Applying the method described above 
to construct the adiabatic regulator, for we find 



(^oo)ad4 ~ 4C J 



(2^ 



Ml(C'f 9(C" 



/\4 



AC 



64C4 



C'C" 



+ 



4 



4C3 
t2 



16C2 8^2 



A3 



2y 



+ 10C{C'fC" + 2C'^{C"f - AC'^C'C'"^ 
( Mj \ 



J ^k 



^ V128C2 ) 
AC^C'C" 



1 



+ 



5((7')^ + 40(C')^C" + 2C^{Cy 
7Ml{Cr 



128C 



■(-5(C")2 + 2CC")i 



1024 



CO 



11 1 ' 



(3.77) 



where Uk is defined in Eq. (3.71). For T, we find 
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2C2 2CJQk 
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+ 



^ {-3{Cr + ^CC") + ^ (9(C')^ + AC'C'iMl)' 



2 

2lC{C'fC" + 6C2(C")^ + 4C^(m2)" + SC^C'C" 
1 



+ 



^(-3(C'f + CC7") 

+ ^^^(87(C')^ - UC^C'{Ml)' - 2QSC{C'fC" 
+ mC'^{C"f + 16C^(M|)" + SOC^C'C" - 16C^C"") 

A<dC{C'fC" + lh{C' f{C" f + 20C2C"C'" - C^C"") 

(l5(C")' - SQC{C'fC" + 6C2(C")' + SC^C'C") 



+ 
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+ 
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(i5(c")^-iicc")4r + 



2048 



10 



(3.78) 

In the free-field limit (Ar = 0), the regulators (3.77) and (3.78) agree with the minimal- 
coupling, spatially fiat limit of the adiabatic regulators obtained by Bunch [173]. 

The renormalization procedure for the semiclassical Einstein equation (2.84) can 
now be precisely stated. According to the adiabatic prescription, we define the quan- 
tum energy-momentum tensor by 



(3.79) 



It can be checked that the momentum-integral expressions for the two independent 
components of Eq. (3.79) are convergent. In terms of {T^i,)^-, the total energy- 
momentum tensor (after renormalization) is 



(3.80) 



where (T'^i,)r stands for T^^^, and renormalized quantities are substituted for bare 
quantities. The bare quantities Gbj Ab, 6b) and Cb are now replaced by Ga, A^, 6^, 
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and Cr in the renormalized semiclassical geometrodynamical field equation, 



(3.81) 



3.3.4 Renormalized semiclassical Einstein equation 

Using semiclassical methods to study the dynamics of the inflaton field in FRW space- 
time requires that the Hubble parameter be much less than the Planck mass, H ^ Mp. 
On dimensional grounds, Cr and 6r are likely to be of order h'^Mp^'^, in which case 
R S> Cj^R^, and R 2> bj^R'^^Rap, provided Rap 0. Let us, therefore, set 6r = 
and Cr = 0, and additionally, let us choose Ar = 0, so that Eq. (3.81) becomes the 
renormalized semiclassical Einstein equation (without cosmological constant). 



TOa + TOH-^((¥'l)R 



(3.82) 



Taking the trace of Eq. (3.82) in spatially flat FRW spacetime, we find 



K = SvtGr 



Cr-)H + cr«)R-^((^^)R)^ 



(3.83) 



Recalling that ^r = 0, and using Eq. (3.58), the classical part of the trace of the 
renormalized energy-momentum tensor is given by 



-(</.)2 + 2 m| + 



Ar 12 



(3.84) 



and the quantum trace of the renormalized energy-momentum tensor is given by 



^ a^J (27r)3 



('^*^)ad4 5 

(3.85) 

where (T'^)ad4 is defined in Eq. (3.78). As discussed in Sec. 3.2.2, the 00 component 
of the semiclassical Einstein equation is a constraint, and is given by 



^ } = 87rGR 



(3.86) 
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From Eq. (3.57), the expression for the classical part of the 00 component of the 
renormalized energy-momentum tensor is given by 



and the quantum part of the 00 component of the renormalized energy-momentum 
tensor is given by 



■-00^(14; 

(3.88) 

where (ro^o)ad4 is defined in Eq. (3.77). 

Eqs. (3.83) and (3.75) are coupled differential equations for a and 0, involving 
complex homogeneous conformal-mode functions Uk which satisfy Eq. (3.76). The 
conformal mode functions are related to the variance (<^|)r by Eq. (3.73). This is 
a closed, time-reversal-invariant system of equations. The initial data at 770 must 
satisfy the constraint equation (3.86). We now drop all "R" subscripts, because we 
will henceforth work only with renormalized quantities. 

3.3.5 Reduction of derivative orders 

The adiabatic regulators (3.69), (3.77), (3.78) for the variance and energy-momentum 

-2 

tensor contain derivatives of up to fourth order in a and up to second order in cp and 
((/^l). The presence of the former can be understood as resulting in part from the well- 
known trace anomaly for a quantum field in curved spacetime [175], which contains 
higher-derivative local geometric terms, e.g., CR. In addition, there are nonanomalous 
finite terms which result from the renormalization of the energy-momentum tensor and 
the choice of minimal coupling. 

The effect of higher derivatives in the semiclassical Einstein equation has been 
much studied in the literature [176-180]. The higher-derivative evolution equations 
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for a and (p have a much larger solution space than the classical Einstein and mean- 
field equations, and in general, the higher-derivative semiclassical Einstein equation is 
expected to have many solutions which are unphysical. In addition, the semiclassical 
Einstein equation (which is fourth order in a) requires more initial data than the 
classical Einstein equation in order to uniquely specify a solution. However, Simon 
and Parker [179,180], following the methods of Jaen, Llosa, and Molina [181], have 
shown that in one-loop semiclassical gravity, there exists a procedure for consistently 
removing the unphysical solutions within the perturbative (H) expansion in which the 
equations arc derived. The procedure corresponds to the addition of perturbative 
constraints, thereby yielding second-order equations which require the same amount 
of initial data as does the classical Einstein equation. Their method involves reducing 
the order of the a'" and a"" terms in the semiclassical Einstein equation using strict 
perturbation theory in h. 

In this study we follow the approach of Simon and Parker to reduce the order of 
the equations for 0, a, and {<^|) to second order. We replace all expressions involving 
a'" and a"" with expressions a'" and a'"' obtained from the classical Einstein equation, 
i.e., Eq. (3.83) with h = 0. This procedure is physically justifiable in this model for the 
following reason: At early times, the dominant contribution to the energy-momentum 
tensor is classical, T^^. Therefore, the deviations a'" — a'" and a"" — a'"', which are 
entirely quantum in origin and oc h, are at early times expected to be very small. In 
addition, at late times the Universe is expected to become asymptotically radiation 
dominated, in which case a'" = a"" = 0. The classical approximations to the late- 
time behavior of a'" and a"" should also have this property, regardless of whether the 
mean-field oscillations are harmonic or elliptic. This procedure is, therefore, physically 
justifiable in the system studied here. 
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3.4 Analysis 



Having derived coupled dynamical equations (3.75), (3.83), (3.76) for the mean field 
0, scale factor a, and conformal-mode functions Uk, respectively, we now proceed to 
solve them. 

3.4.1 Initial conditions 

At the Cauchy hypersurface at r/o, we specify initial conditions on the conformal-mode 
functions which correspond to a choice of quantum state for the fiuctuation field 
<^jj. Based on the analysis in Sec. 3.2.3, we choose boundary conditions at 770 which 
correpsond to the adiabatic vacuum state for at 77 ^ — 00. Prom the semiclassi- 
cal Einstein equation (2.87), the slow-roll condition (3.32), the potential-dominated 
condition (3.31), and assuming that the variance ((/?|) satisfies 



for r] < rjQ, it follows that the spacetime is asymptotically de Sitter at conformal-past 
infinity. Using the approximate solution (3.33) for the scale factor for 77 < r/o, we can 
solve the mode function equation (3.76) for rj < rjo at the same (0th) adiabatic order. 
The general solution is of the form 



(3.89) 
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+ clHi^^{k{r,-H-\r,)-r]o)} 



where H^^^ and H^^^ are the Hankel functions of first and second kind, respectively 



[170], and v is defined by 
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(3.91) 



The function H{ri) is defined as in Eq. (3.34), 




(3.92) 
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where now pc = o?Tqq. The Hubble parameter must be slowly varying for this ap- 
proximation to hold, i.e., the expansion rate nonadiabaticity parameter [41] 

= ^ < 1. (3.93) 

The Wronskian condition on the mode functions (which comes from the canonical 
commutation relations for the fluctuation field operator) requires that 

141' + 141' = 1- (3.94) 

By choosing and c|, different vacua are obtained. The Oth-order adiabatic vacuum 
(matched aX rj = — oo) is constructed by choosing c\, and so that % smoothly 
matches the positive-frequency Oth-order WKB mode function at r] = — oo. This 
corresponds to = 1 and = 0, for all k. Using the asymptotic properties of the 
Hankel function, the adiabatic limit A;, \r]\ — > oo, can be derived, and verified to have 
the correct form, 

lim % ~ ^6"'*=". (3.95) 
fc,|T)|-»oo v2A; 

In addition, the high-momentum, flat-space limit {k, — > oo) gives the same result. 
The initial conditions for the Uk at rjQ are then defined by demanding that the Uk 
functions smoothly match the approximate adiabatic mode function solutions (for 
rf < rfo) at r) = rjQ. This leads to the following initial conditions for the confer mal- 
mode functions: 



4(r/o) = ^ 



(3.97) 



where Hq = H(r]o). The above conditions are valid only at 0th order in the above- 
defined adiabatic approximation, where terms of order H'/H are discarded. It is 
straightforward to show that Eq. (3.93) is valid given the slow-roll (3.32) and inflation 
(3.31) assumptions. In addition to the initial conditions for Uk at 770, we may freely 
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choose initial values for (p{r]o) and (j) (rjo), subject to the constraint that (j) must be 
small enough that conditions (3.31) and (3.32) are valid. We are already assuming 
that a{r]o) = 1. The initial value of a'{r]o) is fixed by the constraint equation (3.86). 

3.4.2 Numerical solution 

In this section we describe the methods we used to solve the coupled evolution equa- 
tions for 4> [Eq. (3.75)], a [Eq. (3.83)], and % [Eq. (3.76)] numerically.^ We evolved 
a representative sampling of mode functions Uk for the region of momentum space 
< A; < Ka, where X is a physical upper momentum cutoff.^ Employing a physical 
[183], as opposed to comoving, momentum cutoff is necessary because a comoving 
cutoff would require the use of the renormalization group equation to track how the 
renormalized parameters flow as the scale factor a increases at each time step. For 
a comoving cutoff the qiiadratic divergence in the variance would be proportional to 
1/a^, requiring a time-dependent renormalization (see [120], for example). The use of 
a physical upper momentum cutoff yields a quadratic divergence which can be removed 
by a non-time-dependent mass renormalization [183]. 

We chose a variety of values of K/m between 50 and 70. The sampling of momentum- 
space is carried out with a uniform binning, with total number of bins A^bins- Various 
values of A^^^ins were used, all greater than 10^. Eq. (3.74) was solved by iteration, and 
the momentum space integrations were performed numerically using the 0{1/N^.^J 
extended Simpson rule. The differential equations (3.83) and (3.75) were evolved us- 
ing 4th-order Runge-Kutta with adaptive step-size control; the target precision for the 
time steps varied between 10~^ and 10~^. Cutoff independence was verified a poste- 
riori by explicitly checking that the results of the numerical solution were insensitive 

^Henceforth, we set /i = 1 and work in units of energy where m = 1. 

finite momentum cutoff is necessary in any case due to the crossing of the Landau point (in the 
large- iV approximation to the theory) when the cutoff is taken to infinity [117, 182]. 
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to a rescaling of K/m. The solutions were computed to a conformal-time scale of 
400 m~^. A typical solution computed according to the above methods required on 
the order of 300 h of CPU time on a modern workstation. 

3.4.3 Results 

A primary goal of this work is the quantitative study of the effect of spacetime dy- 
namics on the parametric resonance energy-transfer mechanism in nonequilibrium 
zero-mode oscillations of a quantum field. As discussed in Sec. 3.2.4, this energy 
transfer, and the corresponding damping of the mean field due to back reaction, occur 
on a time scale of order t\ defined in Eq. (3.41). We numerically evolved the evo- 
lution equations for a, 0, and [ip^ for various values of Mp/m, ranging from very 
large values (corresponding to Minkowski space), to small values (corresponding to 
a strong-curvature, rapid-expansion regime). Figs. 3.1-3.19 show the resulting time 
dependences for the mean field ^, the scale factor a, the variance, A(v?|)/2, the energy 
density p, the energy density in quantum modes pg [defined in Eq. (3.61)], and the 
pressure-to-energy-density ratio 7. The different solutions plotted correspond to 
different values of Mp/m, with A = 10^^"^, K/m = 50, and (j){r]o)/m = 2.0 x 10^. As 
discussed in Sec. 3.4.2, a physical momentum cutoff K was used. The values chosen 
for Mp/m were 10^^, 10^^, 6 x 10^°, and 6 x 10^. The choice of 0(r?o) and A fixes rjo by 
Eq. (3.33) and Hq by Eq. (3.34). Table 3.1 shows (in units where m = 1) the values 
of Mp, the inverse Hubble constant H^^{riQ), and the figure numbers in which the 
corresponding solutions are plotted. The H~^{'r]o) column is the initial inverse Hubble 
constant, which gives the initial time scale for cosmic expansion. Figs. 3.1-3.19 plot 
the resulting solutions. 

The time scales defined in Sec. 3.2.4 can now be explicitly computed. Using Eqs. 
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Figure 3.1: Plot of 4> vs 77, with Mp/m = 1.0 x 10^^. 
(3.37) and (3.34), we have /(770) = \/2, po = 0O) and 

Using Eq. (3.36), we find tq ~ 4.11832 m^^ The value of ti predicted by Eq. (3.41) 
is 132.624 m~^, which is very close to the value predicted by Eq. (3.43), 132.759 mT^. 
For the cases Mp/m = 10^^ and 10^^, it is clear from Table 3.1 that Hq^ S> ti, so 
that the effect of cosmic expansion is expected to be insignificant on the preheating 
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1 X 10" 


-14 
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50.0 


103.65 



Table 3.1: Values of parameters for numerical solutions of Eqs. (3.75), (3.83), (3.76). 
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Figure 3.2: Plot of A((^2)/2 vs r?, with Mp/m = 1.0 x 10^^. 

time scale ri. For the case Mp = 6 x 10^°, 1/{Hqti) ~ 7.8, so that the effect of 
cosmic expansion should be apparent and non-negligible. For the case Mp = 6 x 10^, 
l/{HoTi) ~ 0.78, and cosmic expansion should have a significant effect on parametric 
amplification of quantum fiuctuations. 

Figs. 3.1-3.8 show the dynamics of the mean field and variance in the regime of 
very weak cosmic expansion, ^ n. As expected, under the infiuence of the 
elliptically oscillating mean field, the variance (v^h) gi'ows exponentially in time until 
X{ipl)/2 is of the same order as + /2, at which point back reaction shuts off 
the resonant transfer of energy to the inhomogeneous modes. The time scale for the 
variance to become of order unity can be clearly seen to be ~ ri. As seen previously 
in studies of preheating dynamics in Minkowski space [99] , on the time scale ~ ri , the 
mean field decouples from its own fiuctuations and oscillates with an asymptotically 
finite amplitude, given by [99] A0^(2m2 ) = 0.914. In the Minkowski space limit, 
corresponding to Mp/m — ^ oo, covariant conservation of the energy-momentum tensor 
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Figure 3.3: Plot of 7 vs 77, with Mp/m = 1.0 x 10^^. 

implies that dp/dt = 0. This was verified for the case of Mp/m = 10^^, where no change 
in p was detected to within the numerical precision of our algorithm, as expected from 
dimensional analysis of Eq. (3.83). The increase in the scale factor for these cases was 
within a few parts in 10^ of the initial value 0(770) = 1. The asymptotic equation of 
state plotted in Fig 3.13 is observed to be 7 ~ 0.18. This is exactly what would be 
predicted for a two-fiuid model consisting of a mean field with equation of state given 
by Eq. (3.42), 7c ^ 0.0288, and a relativistic gas corresponding to the energy density 
of the (f field, with 7q ~ 0.3333. The average 7q + 7c = 0.182. 

For the case Mp/m = 6 x lO^'^, the effect of cosmic expansion is clearly visible in 
Figs. 3.9-3.14. In Fig. 3.9, the coherent oscillations of the mean field for the time 
period < r/ — ryo < ~ 27ro are clearly seen to be redshifted by the usual 1/a factor 
expected from the Hubble damping term in Eq. (3.75). The expected asymptotic 
equation of state (taking into account cosmic expansion) computed from a simple 
two-fiuid model is ~ 0.133, in agreement with Fig. 3.13. 
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Figure 3.4: Plot of pq vs 77, with Mp/m = 1.0 x 10^^. 

Figs. 3.15-3.19 show the solution for Mp/m = 6 x 10^. In this case, 1/{HqTi) ~ 
0.781. From Fig. 3.17, we clearly see that cosmic expansion renders parametric am- 
plification of quantum fluctuations an inefficient mechanism of energy transfer to the 
inhomogeneous modes. The very rapid oscillations of the mean field at late times are 
due to the conformal time scale used here, in which the oscillation period of the mean 
field decreases inversely with a. Damping of the mean field due to cosmic expansion is 
the dominant effect in Fig. 3.15. The power-law decrease in energy density consistent 
with matter having an effective equation of state 7 ~ 0.0288 can be seen in Fig. 3.18. 
At rj = 300 m~^, the ratio Pq/p ~ 0.0002, so the fraction of energy density in the inho- 
mogeneous modes is negligible in comparison to the classical, mean-field contribution. 
Since the variance is never large enough that it dominates the effective mass M, 
the mode functions approximately obey the one-loop equation, in which the effective 
frequency is k"^ + a'^{m'^ + X(f) ), neglecting the a" /a term. The width of the resonance 
can then be shown to be approximately given by k"^ < X(j)o/2. The variance is damped 
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Figure 3.5: Plot of 4> vs 77, with Mp/m = 1.0 x 10^^. 

by 1/a^ due to cosmic redshift, so when ~ ri, the variance never grows to be of 
order unity. 

In addition to varying Mp, the coupling A was varied, with results in agreement 
with Eq. (3.41), showing a logarithmic dependence of ri on A~^. 

3.5 Summary 

In this Chapter we use the minimally coupled, quartically self- interacting scalar 0{N) 
field theory as a model for the inflaton field, and study its nonequilibrium dynam- 
ics nonperturbatively in a spatially flat FRW spacetime whose evolution is driven by 
the quantum field. We solve the coupled, self-consistent semiclassical Einstein equa- 
tion, mean-field equation, and conformal-mode-function equations numerically. Our 
goal in this Chapter is to study the effects of spacetime dynamics on the mean field, 
and parametric amplification of quantum fluctuations. This process of energy trans- 
fer from the mean field to the inhomogeneous modes is inherently nonperturbative 
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Figure 3.6: Plot of A((^2)/2 vs r?, with Mp/m = 1.0 x 10^^. 

and nonequilibrium. It requires the use of the closed-time-path formalism and the 
two-particle-irreducible effective action. As our focus in this Chapter is on the para- 
metric amplification of quantum fluctuations, we assume unbroken symmetry. Our 
analysis is, therefore, most relevant to reheating in chaotic inflation scenarios. We 
use the two-loop, covariant equations for the mean field and the two-point function 
for the fluctuation field derived in Chapter 2 and study the case of leading order in 
the 1 /N expansion, an approximation which is valid on time scales much shorter than 
the mean-free time for multiparticle scattering (T2). For FRW spacetimes, we use the 
well-established adiabatic regularization procedure to obtain finite expressions for the 
renormalized variance and energy-momentum tensor which enter into the mean-field 
equation, conformal-mode function equations, and the semiclassical Einstein equation. 
In our approach, covariant conservation of the energy-momentum tensor is preserved 
at all times, as it should be. (It should not and need not be put in by hand, as 
was done in a recent study of reheating in a fixed background FRW spacetime [162].) 
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Figure 3.7: Plot of 7 vs ij, with Mp/m = 1.0 x 10^^. 

We use the adiabatic vacuum construction (matched at conformal-past infinity) to 
define the quantum state for the fluctuation field in FRW spacetime with asymptotic 
de Sitter initial conditions; this is the most physical vacuum construction given the 
decidedly nonadiabatic conditions which prevail at the end of inflation. The instan- 
taneous Hamiltonian diagonalization constructions used in earlier studies of reheating 
in curved spacetime [116, 124] are known to be problematic [130]. 

We evolved the coupled dynamical equations for the mean field, variance, and scale 
factor using standard numerical methods, for time scales of 400 m~^, where the initial 
period of mean-field oscillations is 4.11832 m"^. Several regimes for the parameters 
of the system were investigated. From the solutions of the dynamical equations we 
studied the behavior of the scale factor a, the mean field cj), the energy density p, 
pressure-to-energy-density ratio 7, and the inhomogeneous-mode (fluctuation-field) 
energy density pq. The solutions of the dynamical equations were analyzed for a 
variety of values for MpTo, the parameter which controls the rate of cosmic expansion 
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Figure 3.8: Plot of pq vs 77, with Mp/m = 1.0 x 10^^. 

relative to the time scale for mean-field oscillations in the model. In the case of 
negligible cosmic expansion, corresponding to very small initial inflaton amplitude, 
the dynamics is identical to that seen in the group 2B (see Sec. 1.2) studies of O(A^) 
preheating in Minkowski space [99]. In particular, the conservation of energy and 
logarithmic dependence of the preheating time scale ri on the inverse coupling 
[as shown in Eq. (3.41)] are confirmed. For the case of moderate cosmic expansion, 
H{rio)Ti ~ 10 [where H{rjo) is the Hubble parameter at the initial time 770], energy 
transfer via parametric amplification of quantum fluctuations is still efficient, and the 
dynamics can be understood using the analytic results of [99] (for Minkowski space), 
in terms of the conformally transformed mean-field amplitude ^ ^/a and oscillation 
period tq = To/a. The asymptotic effective equation of state is found to be consistent 
with the prediction of a simple two-fluid description of the late-time behavior of the 
system. 

The most significant physical result concerns the case of rapid cosmic expansion. 
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Figure 3.9: Plot of (f) vs r), with Mp/m = 6.0 x 10^°. 

where H~^{r]Q) ~ ri. In this case we find that parametric ampHfication of quantum 
fluctuations (via parametric resonance) is an inefficient mechanism of energy transfer 
to the inhomogeneous modes of the inflaton, because the parametric resonance effect 
is inhibited both by redshifting of the mean-field amplitude and by the redshifting of 
the physical momenta of the modes out of the resonance band. The energy density of 
particles produced through parametric resonance is in this case redshifted so rapidly 
that, in our model, the term A((/?|)/2 never grows to be of the order of the tree-level 
effective mass, m? + Xcp /2. As the mean-field amplitude is damped (oc 1/a) due to 
cosmic expansion, eventually the resonant particle production ceases, and the mean 
field oscillates with a damped envelope at late times. This leads us to the following 
conclusions: (i) On the physical level, in chaotic inflation scenarios with a inflaton 
minimally coupled to gravity and with a large initial inflaton amplitude at the end of 
slow roll, parametric amplification of the infiaton's own quantum fiuctuations is not 
a viable mechanism for reheating the Universe, unless the self-coupling is significantly 
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Figure 3.10: Plot of a vs r?, with Mp/m = 6.0 x 10^°. 

increased.^ This does not imply that the phenomenon of parametric amplification 
of quantum fluctuations does not play a vital role in the "preheating" period of in- 
flationary cosmology, for different models and/or couplings. The interesting case of 
a (fy^x^ model is currently under investigation, (ii) On a more methodological level, 
we conclude that a correct study of the reheating period in a chaotic inflation model 
with large inflaton amplitude at the onset of reheating must take into account the 
effects of spacetime dynamics. This should be carried out self- consistently using the 
coupled semiclassical Einstein equation and matter-field equations, so that no ad hoc 
assumptions need be made about the effective equation of state and/or the relevant 
time scales involved. 

A full two-loop treatment of the unbroken symmetry mean-field dynamics of the 
0{N) field theory [which involves solving the nonlocal, integro-differential equations 
(2.116) and (2.117)] includes multiparticle scattering processes, which provide a mech- 
anism for reheating; but they are of a qualitatively different nature than the parametric 

^In recent work on galaxy formation from quantum fluctuations, Calzetta, Hu, and Matacz [54, 55] 
report that A can be as high as ~ 10~^. 
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Figure 3.11: Plot of X{^'^)/2 vs r/, with Mp/m = 6.0 x 10^°. 

resonance energy-transfer mechanism studied here. In addition, the nonlocal nature of 
the gap equation in the full two-loop analysis makes numerical solution of the coupled 
Einstein and matter equations difficult. In our model, scattering occurs on a time scale 
T2 which is significantly longer than the time scale ti for parametric amplification of 
quantum fluctuations. Therefore, in this model the leading-order, large-iV (collision- 
less) approximation is sufficient for a study of parametric amplification of quantum 
fluctuations. In addition, realistic models of inflation invariably involve couplings of 
the inflaton to other fields, which provides additional mechanisms of energy transfer 
away from the inflaton mean field, and into its (or other fields') quantum modes. 

The issues involved in a systematic study of the thermalization stage of post- 
inflationary physics are more complex. A quantum kinetic field theory treatment 
taking into account multiparticle scattering is required. The two-loop, 2PI effective 
action is the simplest and most generally applicable rigorous formalism which contains 
the leading-order multiparticle scattering processes. The leading-order, 1/N approxi- 
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Figure 3.12: Plot of p vs r], with Mp/m = 6.0 x 10^°. 

mation is a coUisionless subcase of the two-loop, 2PI effective action; it is employed in 
this study in order to obtain local dynamical equations which can be solved numeri- 
cally, and is adequate for a study of parametric amplification of quantum fluctuations. 
In addition, the growth of entropy must be understood within the context of a physi- 
cally meaningful coarse graining of the full time-reversal-invariant quantum dynamics 
of the field theory. This is discussed in Chapter 5 below. A first-principles analysis of 
the thermalization stage is currently underway [184]. 
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CHAPTER 4 



Fermion production, noise, and stochasticity 

4.1 Introduction 

4.1.1 Background and issues 

In Chapter 3 we studied the effect of parametric particle creation on the inflaton dy- 
namics in the post-inflation, "preheating" stage. It was shown that analysis of the 
dominant physical processes during the early stages of reheating necessitates consid- 
eration of self-consistent back reaction of the inflaton variance on the mean field and 
inhomogeneous modes. This is a nonperturbative effect, and its description requires 
a consistent truncation of the Schwinger-Dyson equations [82,99, 125]. Assuming ini- 
tial conditions conducive to efficient parametric particle creation, the end state of the 
regime of parametric particle creation consists of a large inflaton variance (i.e., on the 
order of the tree level terms in the inflaton's effective mass). Thus, both the inflaton 
mean field and variance should be treated on an equal footing. This requires a two- 
particle-irreducible (2P1) formulation of the effective action which is a subclass of the 
master effective action [82]. 

During the later stages of reheating, the dynamics of the inflaton field is thought, 
in the case of unbroken symmetry, to be dominated by fermionic particle creation 
[166]. This stage of infiaton dynamics is the subject of this chapter. We consider 
a model consisting of a scalar inflaton field (with A^^ self-coupling) coupled to 
a spinor field ip via a Yukawa interaction and we attempt to present as complete 
and rigorous a treatment as mandated by the self-consistency of the formalism and 
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the actual solvability of the equations. Thus, we adopt a closed-time-path (CTP), 
two-particle-irreducible (2PI), coarse-grained effective action (CGEA) to derive the 
dynamics of the inflaton field. We have explained the significance of CTP and 2PI in 
Chapter 2 [67-69, 80, 82, 84, 185], and their relevance to the study of inflaton dynamics 
in Chapter 3. Here, an added feature of an open system is introduced: we wish to 
include the averaged effect of an environment on the system, and a useful method is 
via the coarse-grained effective action [186-190]. Let us explain the rationale for this. 

4.1.2 Coarse- Grained Effective Action 

In inflationary cosmology at the onset of the reheating period, the inflaton field's zero 
mode typically has a large expectation value, whereas all other fields coupled to the 
inflaton, as well as inflaton modes with momenta greater than the Hubble constant, 
are to a good approximation in a vacuum state [8]. This suggests imposing a physical 
coarse graining in which one regards the inflaton field as the system, and the various 
quantum fields coupled to it as the environment. Prom the closed-time-path, coarse- 
grained effective action (CTP-CGEA) [186-190] derived in Sec. 4.2 below, one obtains 
effective dynamical equations for the inflaton field, taking into account its effect on 
the environment, and back reaction therefrom. For the present problem, the system 
consists of the inflaton mean field and variance, and the environment consists of the 
spinor field(s) coupled to the system via a Yukawa interaction. 

We wish to emphasize here a subtle yet important distinction between the system- 
environment division in nonequilibrium statistical mechanics and the system-bath di- 
vision assumed in thermal field theory. In the latter, one assumes that the propagators 
for the bath degrees of freedom are fixed, finite-temperature equilibrium Green func- 
tions, whereas in the case of the CTP-CGEA, the environmental propagators arc slaved 
(in the sense of [82]) to the dynamics of the system degrees of freedom, and are not 
fixed a priori to be equilibrium Green functions for all time. This distinction is im- 
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portant for discussions of fermion particle production during reheating, because it is 
only when the inflaton mean field amplitude is small enough for the use of perturba- 
tion theory, that the system-bath split implicit in thermal field theory can be used. 
Otherwise, one must take into account the effect of the inflaton mean field on the bath 

(spinor) propagators. 

4.1.3 Earlier work 

The first studies of particle production during reheating in inflationary cosmology 
were [91,93], where reasonable estimates of particle production were made, but back 
reaction effects were not addressed. The earliest studies of fermionic particle creation 
during reheating used time-dependent perturbation theory to compute the imaginary 
part of the self-energy for the zero mode of the inflaton field, which was related to the 
damping parameter in a friction-type phcnomcnological term in the equation of motion 
for the inflaton zero mode [91,92,94,95, 191]. In these studies, it was assumed that 
the effect of fermionic particle production on the dynamics of the inflaton zero mode 
is that of a T4> friction term. However, it has been shown [69, 87, 90, 106-108, 112, 192] 
that this assumption is not tenable for a wide variety of field-theoretic interactions and 
initial conditions. Rather, the effect of back reaction from particle production must be 
accounted for by deriving the effective evolution equation for the quantum expectation 
value of the inflaton zero mode, where the dynamics of the degrees of freedom of the 
produced particles are either coarse-grained (as in Sec. 4.3), or accounted for through 
self-consistent coupled equations (as in Sec. 4.2). In general, particle creation leads to 
a nonlocal dissipation term in the inflaton mean field equation, and it is only under 
rather idealized conditions and specialized cases that one can expect the dissipation 
kernel to approach a delta function (i.e., a friction term) [53,56,74,75]. Therefore, 
[91,94,95] missed the time-nonlocal nature of fermion particle production and its 
effect on the dynamics of the inflaton field. In addition, these studies computed the 
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self-energy in flat space, thereby neglecting the effect of curved spacetime on fermion 
production, and did not examine the stochastic noise arising from the coarse graining 
of the fermion field. 

In addition, most early studies of fermion production during reheating, in using 
time-dependent perturbation theory to compute the vacuum particle production rate 
[91,94,95], did not include the effect of back reaction of the produced fermion parti- 
cles on the particle production process itself. In [112], the effect of a thermal initial 
fermion distribution on the particle production process was investigated (and a Pauli 
blocking effect was shown), but their analytic derivation of the Pauli blocking effect 
involves the same perturbative expansion (i.e., system-bath split) described above, 
and therefore does not incorporate the effect of the produced fermion quanta on the 
particle production process. In order to take this effect into account, it is necessary 
to include the effect of the time- varying inflaton mean field in the equation of motion 
for the spinor propagator, which amounts to a coarse graining of the fermion field, 
in the system-environment sense, as described above. The perturbative amplitude 
expansion of the effective inflaton dynamics, in contrast, amounts to a system-bath 
coarse graining which does not include this back reaction effect. 

More recent studies of fermion production during reheating [110-112] obtained 
dynamical equations for the inflaton mean field using a one-loop factorization of the 
Lagrangian, and solved them numerically. However, these studies were carried out 
in flat space, and because they studied only the dynamics of the inflaton mean field 
(and at one loop), their analysis did not take into account the back reaction of the 
inflaton variance on the fermion quantum modes, nor the back reaction of particle 
production on the dynamics of the inflaton two-point function. The importance of 
the curved spacetime effect was addressed in Chapter 3, and we will discuss below the 
importance of treating the inflaton quantum variance on equal footing with the mean 
field. 
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4.1.4 Present work 

In this study, we wish to describe the late stages of the reheating period, in which the 
damping of the inflaton mean field is dominated by fermionic particle production; our 
starting point is the end of the preheating period (in which the inflaton dynamics was 
dominated by back reaction from parametric particle creation, as discussed in Chap- 
ter 3. Because the inflaton variance {ip^) can (for sufficiently strong self-coupling) be 
on the order of the square of the amplitude of mean field oscillations at the end of 
preheating (in the case of unbroken symmetry) [125], it is necessary to treat the infla- 
ton mean field and variance on an equal footing in a study of the subsequent effective 
dynamics of the inflaton field. This requires a two-loop, two-particle-irreducible for- 
mulation of the coarse-grained effective action. At two loops, both the inflaton mean 
field and the inflaton variance couple to the spinor degrees of freedom, and are damped 
by back reaction from fermion particle production; all the previous studies mentioned 
above, in using the IPI effective action, missed this possibly important effect. 

In addition to having a large variance, the inflaton amplitude at the end of the 
preheating period may be large enough that the usual perturbative expansion in powers 
of the Yukawa coupling constant is not valid [see Eq. (4.49) below], in which case a 
nonperturbative derivation of the inflaton dynamics is required. In the construction of 
the CTP-2PI, coarse-grained effective action below, the spinor propagators obey one- 
loop dynamical equations in which the inflaton mean field acts as an external source. 
Studies of fermion particle production during reheating which relied on the use of 
perturbation theory in the Yukawa coupling constant [87, 91-95, 191] therefore do not 
apply to the case of fermion particle production at the end of preheating with unbroken 
symmetry, when the Yukawa coupling is sufficiently large. The dynamical equations 
derived in Sec. 4.2 below for the inflaton mean fleld and variance are applicable even 
when, as may be the case, the inflaton mean field amplitude is large enough that a 



114 



perturbative expansion in powers of the Yukawa coupling is not justified. 

Although, as discussed above, a proper treatment of the early stage of fermion 
production during reheating, starting at the end of preheating with unbroken symme- 
try, should in principle employ the CTP-2PI-CGEA to obtain coupled equations of 
motion for the inflaton mean field and variance, at very late times the inflaton mean 
field and variance will have been damped sufficiently (due to the dissipative mecha- 
nisms derived below in Sec. 4.2) that the perturbative IPI effective action will yield 
a qualitatively correct description of the inflaton mean-field dynamics. While curved 
spacetime effects should in principle be incorporated self-consistently for a quantita- 
tive calculation of the reheating temperature in a particular inflationary model (as 
discussed in [125]), for a general discussion of dissipative effective dynamics of the in- 
flaton mean field in the case of weak cosmic expansion (where the Hubble constant H 
is much smaller than the frequency of inflaton oscillations), it is reasonable to neglect 
curved spacetime effects in computing the spinor propagators. Therefore in Sec. 4.3, 
we derive the perturbative, flat-space CTP-IPI-CGEA to fourth order in the Yukawa 
coupling constant, and obtain an evolution equation for the inflaton mean field with 
nonlocal dissipation. 

Another new feature of our work obtainable only from the stochastic approach 
adopted here is the derivation, in Sec. 4.4, of a Langevin equation for the inflaton mean 
field, with clear identification of the dissipation and noise kernels from the CGEA. 
We have calculated the energy dissipated and the fluctuations in the energy. From 
the latter we obtain the range of parameters where the conventional "mean-field" 
approach breaks down. We believe the methodology presented here provides a better 
theoretical framework for the investigation of phase transitions in the early universe, 
as exemplified by our treatment of reheating in inflationary cosmology. 
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4.1.5 Organization 

Our work is organized as follows. In Sec. 4.2, we derive the coupled equations of motion 
for the inflaton mean field and variance, in a general curved spacetime, including 
diagrams in the coarse-grained CTP-2PI effective action of up to two-loop order. In 
Sec. 4.3, we specialize to Minkowski space and small mean-field amplitude, and obtain 
a perturbative mean field equation including dissipative effects up to fourth order 
in the Yukawa coupling constant. In Sec. 4.4, we examine the dissipation and noise 
kernels obtained in Sec. 4.3, and show that they obey a fluctuation-dissipation relation. 
We then derive a Langevin equation which self-consistently includes the effect of noise 
on the dynamics of the inflaton zero mode. We summarize our results in Sec. 4.5. 

4.2 Coarse-grained inflaton dynamics in curved spacetime 

In this section, we present a first-principles derivation of the nonequilibrium, nonper- 
turbative, effective dynamics of a scalar inflaton field (j) coupled to a spinor field ip via 
a Yukawa interaction, in a general, curved classical background spacetime. The use of 
the Schwinger-Keldysh closed-time-path (CTP) formalism allows formulation of the 
nonequilibrium dynamics of the inflaton from an appropriately defined initial quantum 
state. The evolution equations for the inflaton mean field and variance are derived from 
the two-loop, closed-time-path (CTP), two-particle-irreducible (2PI), coarse-grained 
effective action (CGEA). As the name suggests, there are two approximations of a 
statistical mechanical nature. One is the coarse graining of the environment — here 
the inflaton field is the system and the fermion field is the environment [186]. The 
other refers to a truncation of the correlation hierarchy for the inflaton field [82] — the 
two-particle-irreducible effective action. This formulation retains the inflaton mean 
field and variance as coupled dynamical degrees of freedom. Back Reaction of the 
scalar and spinor field dynamics on the spacetime is incorporated using the semiclas- 
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sical Einstein equation, which foUows from functional differentiation of the effective 
action with respect to the metric. It is shown that these dynamical equations are both 
real and causal, and the equations are cast in a form suitable for implementation of 
an explicit curved-spacetime renormalization procedure. 



4.2.1 The model 

We study a model consisting of a scalar field ^ (the inflaton field) which is Yukawa- 
coupled to a spinor field tp, in a curved, dynamical, classical background spacetime. 
The total action 

S[4>,i;,i;,gf^^] = S^[g'^''] + S^[<j>,i^,i;, g^, (4-1) 
consists of a part depicting classical gravity, 5° [5^^], and a part for the matter fields. 



5^[<^,v;,v,5n = s'^[ct>,gn+s'^[ip,i^,gn+s''[ct^,i^,i^,gn, 



(4.2) 



whose scalar (inflaton), spinor (fermion), and Yukawa- interaction parts are given by 



d ^x\J—g 



d x\/ —g 



0(n + m2 + ei?)</. + -</)^ 



d '^x^/—g(|)'^|J'tp. 



(4.3) 
(4.4) 
(4.5) 



For this theory to be renormalizable in semiclassical gravity, the bare gravity action 
S°[g^"'] of Eq. (4.1) should have the form given by Eq. (2.65) [17, 129]. In Eqs. (4.3)- 
(4.5), m is the scalar field "mass" (with dimensions of inverse length); ^ is the dimen- 
sionless coupling to gravity; /j, is the spinor field "mass," with dimensions of inverse 
length; □ is the Laplace-Beltrami operator in the curved background spacetime with 
metric tensor g^,^; is the covariant derivative compatible with the metric; y/—g 
is the square root of the absolute value of the determinant of the metric; A is the 
self-coupling of the infiaton field, with dimensions of 1/Vh; f is the Yukawa coupling 
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constant, which has dimensions of and R is the scalar curvature. The curved 

spacetime Dirac matrices satisfy the anticommutation relation 

{7^7"}+ = 25'^^l.p, (4.6) 

in terms of the contravariant metric tensor g^'^ . The symbol l^p denotes the identity 
element in the Dirac algebra. 

It is assumed that there is a definite separation of time scales between the stage of 
"preheating" discussed in Chapter 3, and fermionic particle production, which is our 
primary focus in this work. However, this does not imply that perturbation theory in 
the Yukawa coupling constant / is necessarily valid, which would require that condition 
(4.49) (defined in Sec. 4.3 below) be satisfied. In addition, the fermion field mass ^ 
is assumed to be much lighter than the inflaton field mass m, i.e., the renormalized 
parameters m and // satisfy m 3> 

4.2.2 Closed-time-path, coarse-grained effective action 

We denote the quantum Heisenberg field operators of the scalar field 4> and the spinor 

field 'ip by $h and ^'h, respectively, and the quantum statc^ by For consistency 

with the truncation of the correlation hierarchy at second order, we assume $h to 

have a Gaussian moment expansion in the position basis [118, 150], in which case the 

relevant observables are the scalar mean field 4> [defined in Eq. (3.2)], the mean-squared 

fluctuations, or variance [defined in Eq. (3.4)]. As discussed above, at the end of the 

preheating period, the inflaton variance can be as large as the square of the amplitude 

of mean-field oscillations. On the basis of our assumption of separation of time scales 

in Sec. 4.2.1, and the conditions which prevail at the onset of reheating, the initial 

quantum state |0) is assumed to be an appropriately defined vacuum state for the 

^Although in this case the particular initial conditions constitute a pure state, this formalism can 
encompass general mixed-state initial conditions [68]. 
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spinor field. 

The construction of the CTP-2PI-CGEA for the (p'ijjip theory in a general, curved, 
background spacetime closely parallels the construction of the CTP-2PI effective ac- 
tion for the 0{N) model in curved spacetime discussed in Chapter 2. Following the 
approach of Sec. 2.4, we define a "CTP manifold" Ai, and a volume form ejn on Ai, 
in terms of the discrete set {+, — } labeling the "time branch," and a Cauchy hyper- 
surface which is far to the future of the time scales in which we are interested. The 
spacetime manifold M is the past domain of dependence of Sy^. The restrictions of a 
function (p, defined on A^, to the + and — time branches are denoted by 0+ and 
respectively. We can then define a matter field action on M, 

(4.7) 

where the spacetime integrations in are now over M only. We use the symbol to 
distinguish it from the action on M. Let us also simplify notation by suppressing 
time branch indices in the argument of functionals on A^, i.e.. 

Let us also define the functional 5^ on M. by 

S^[ct>,iP,iP,gn = S-[,^+,V;+,V+,5r] -'5"['^-,VS-,V'-,5n, (4.9) 

in analogy with Eq. (4.7). For a function <^ on A^, the restrictions of cf) to the -|- and 
— time branches are subject to the boundary condition 

(</'+)|s. = (</'-)|s. (4.10) 

at the hypersurface S^. The gravity action 5°, promoted to a functional on M, takes 
the form 

= ^'^bf] - S°[gn, (4.11) 
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where the range of spacetime integration in S'^ on the right-hand side of Eq. (4.11) is 
understood to be over M. 

To formulate the CTP-2PI-CGEA, our first step is to define a generating functional 
for n-point functions of the scalar field, in terms of the initial quantum state which 
evolves under the influence of a local source J, and a non-local source K coupled to the 
scalar field (in the interaction picture with the external sources being treated as the 
"interaction"). This generating functional depends on both J and K, as well as the 
classical background metric g^'^. In the path integral representation, the generating 
functional Z[J, K, g^^'^] takes the form of a sum over scalar field configurations (p and 
complex Grassmann- valued configurations ip on the manifold M, 

\ I d^x^ I d^x'^c''^c'^Kac{x,x')Ux)cl>a{x') 



' ctp 

ab 



+ I d'^xV^d^'JaCt 
M 



(4.12) 

where Ja{x) is a local c-number source and Kab{x,x') is a nonlocal c-number source. 
The subscript CTP on the functional integral denotes a summation over field configu- 
rations (f)±, -ijj^, and 'tp± which satisfy the boundary condition (4.10). The latin indices 
a, 6, c, . . . , have the discrete index set {-|-, — }, and denote the time branch [67, 68]. 
The boundary conditions on the functional integral of Eq. (4.12) at the initial data 
surface determine the quantum state |ri). The CTP indices have been dropped from 
gi^" for ease of notation; it will be clear how to reinstate them [80] in the two-loop 
CTP-2PI effective action shown below in Sec. 4.3. The two-index symbol c"^ is de- 
fined by the n = 2 case of Eq. (2.27). The generating functional for normalized n-point 
functions is 

W[J,K,g''^] = -imnZ[J,K,g'"'], (4.13) 
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in terms of which we can define the classical scalar field on Ai, 

1 5W[J,K,gn^ ^^^^^ 



^ SJb{x) 

and the scalar two-point function on M in the presence of the sources Ja and Kab, 

hGab{x,x')jK = 2CacCbd—;= I T .jJ \ TV ~ 4>a{x)4>b{x')^ (4-15) 



where the JK subscripts indicate that ^ and G are functional of the Ja and Kab 
sources. In the limit Ja = K^b = 0, the classical field is the same on the two time 
branches and it is equivalent to the mean field 0, as shown in Eq. (2.43). In the same 
limit, Gab becomes the CTP propagator for the fluctuation field defined in Eq. (3.3), 
as shown in Eqs. (2.44)-(2.47). In the coincidence limit^ x' = x, all four components 
(2.44)-(2.47) arc equivalent to the variance (vh^) defined in Eq. (3.4). Provided we 
can invert Eqs. (4.14) and (4.15) to obtain Ja and Kab in terms of and Gab, the 
CTP-2PI eflFective action can be defined as the double Legendre transform (in both 
JaSMdKab) oiW[J,K,gi^% 



r[<^,G,3^'^] = I^[J,K,5n- / d''xV^c'''Jaix)Ux) 
-If d^x^ [ dVV=?c«V'^i^ac(a;,xO [^G6d(x,xO + 06(x)0d(x')l . (4.16) 

The inverses of Eqs. (4.14) and (4.15) can be obtained by functional differentiation of 
Eq. (4.16) with respect to 

1 S^[4>,G,g^^•^] 



-9 54>a{x) 



= -c'^''Jb{x)^a - if'c^" 



and with respect to Gab-, 



M 



d [Kbdix,x')^G + Kdb(.x',x)^G] <Pcix), (4.17) 



/^v^ 5Gab{x,x') - 2" ' ^c.(x,x),„ 



(4.18) 



^The variance {ip(x)^) is divergent in four spacetime dimensions, and should be regularized using 
a covariant procedure [17, 193]. 
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where the <pG subscript indicates that Kab and J a are functionals of (p^^ and Gab- In- 
serting Eqs. (4.17) and (4.18) into Eq. (4.16) yields a functional integro-differential 
equation for the CTP-2PI effective action in terms of and G only, so the JK sub- 
scripts can be dropped. It is useful to change the variable of functional integration 
to be the fluctuation field about defined by Eq. (2.33). Performing the change-of- 
variables Dcf) Df, the equation for F is 

T[lG,gn= [ d'xf dV-^Pf^G„,(x,x') 

Jm Jm dGba{x',x) 



-ihlni / Dip+Dil;^D'tp+Dip^D'tp_D'tp^ex.p 



JM 



4 5r[,A,G,<7n 



which has the formal solution 



ih. 



" JM JM 



X . . -ipaKx)ipb[X , 



5Gba{x',x) 



(4.19) 



r[<^, = S'^[4>\ - ylndetGafe - iMndetF„b + PsI^, G] 



ih 



+ — I d^xy/^ f d^x'^f^A"-^{x',x)Gab{x,x), (4.20) 

JM 

where A°''^ is the second functional derivative of the scalar part of the classical action 



S"**", evaluated at 0, 
1 



iA^^x^x') 



g \6(pa{x)S(j)b{x') 



1 



c"^(n, + m' + ^R{x)) + c«''^'*^0,(x)^,(x 



abed 



A 



S{x — x) 



1 



(4.21) 



The symbol F^b denotes the one-loop CTP spinor propagator, which is defined by 

Fab{x,x') = B-^{x,x'), (4.22) 

where we are suppressing spinor indices, and the inverse spinor propagator B"'^ is 
defined by 



iB^{x,x^) 



1 



1 



6'llJa{x)Stpi,{x') 



c'^{il''V^-fi)-c'^"'mx') S{x'-x) 



(4.23) 



'-9 
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It is clear from Eq. (4.23) that the use of the one-loop spinor propagators in the 
construction of the CTP-2PI-CGEA represents a nonperturbative resummation in the 
Yukawa coupling constant, which (as discussed above) goes beyond the standard time- 
dependent perturbation theory. The boundary conditions which define the inverses of 
Eqs. (4.21) and (4.23) are the boundary conditions at the initial data surface in the 
functional integral in Eq. (4.12), which in turn, define the initial quantum state 
The one-loop spinor propagator F^b is related to the expectation values of the spinor 
Heisenberg field operators for a spinor field in the presence of the c-number background 
field (p, 

nF++(x,x%^^^_^^ = {n\Ti^,ix)^,{x'))\n), (4.24) 
hF__{x,x')\^^^^_^^ = {n\f{^,{x)^,{x'm), (4.25) 
hF+_{x,x')\^^^^_^^ = -(0|*H(x')*H(x)|n), (4.26) 
hF_+{x,x\^^^_^^ = {n\^^{x)^^{x')\n), (4.27) 

where the spinor Heisenberg field operators obey the equations 

{i-f^V^-fi-f4>)^ = 0, (4.28) 
(-^7^V^ -fi- f^)^ = 0. (4.29) 

The CTP spinor propagator components satisfy the relations (valid only when (f>+ = 

F++{x,x'y = F_^{x',x), (4.30) 

F__{x,x'y = F++(x',x), (4.31) 

F^+{x,x'y = F_+(x',x), (4.32) 

F+_{x,x'Y = F+_{x',x). (4.33) 

The functional T2[4>, G] is defined as —ih times the sum of all vacuum diagrams drawn 
according to the following rules: 
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1. Vertices carry spacetime (x G M) and time branch (a G {+, — }) labels. 

2. Scalar field lines denote hGab{x,x'). 

3. Spinor lines denote the one-loop CTP spinor propagator hFab{x,x') (spinor in- 
dices are suppressed), defined in Eq. (4.22). 

4. There are three interaction vertices, given by iS^/h, which is defined by 

S%ip,i>,2l;] = S' , (^+ , V+ , V'+] - ^' [4>- ,^-,4^-,^^-], (4.34) 

(4.35) 



S^[4>,(p,ip,tl)] = - J d^xy^ 
where we have followed the notation of Eq. (4.9). 



5. Only diagrams which are two-particle- irreducible with respect to cuts of scalar 
lines contribute to 

The distinction between the CTP-2PI, coarse-grained effective action which arises 
here, and the fully two-particle-irreducible effective action (2PI with respect to scalar 
and spinor cuts) , is due to the fact that we only Legendre-transformed sources coupled 
to 0; i.e., the spinor field is treated as the environment. In Eq. (4.21), the curved- 
spacetime Dirac S function is defined as in [17]. Comparison of Eq. (4.20) above with 
Eq. (4.13) of Ref. [80] [which was computed for the 0{N) model] shows that the 
TrlnFab in Eq. (4.20) differs from the usual one-loop term by a factor of 2, owing to 
the difference (in the exponent) between the Gaussian integral formulas for real and 
complex fields [79]. 

The functional T2[ct),G,g^^'^] can be evaluated in a loop expansion, which corre- 
sponds to an expansion in powers of h, 

oo 

r2[4>,G,gn = ^/i'rW[<^,G,5n, (4.36) 

1=2 

starting with the two-loop term, F^^^. The functional P^^^ has a diagrammatic expan- 
sion shown in Fig. 4.1, where solid lines represent the spinor propagator F (as defined 
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Figure 4.1: Diagrammatic expansion for F^^^ , the two-loop part of the CTP-2PI coarse- 
grained effective action. 

in Sec. 4.2), and dotted hnes represent the scalar propagator G. The vertices termi- 
nating three lines are proportional to the scalar mean field (/). Each vertex carries 
spacetime (x) and CTP — ) labels. The A(/>'' self-interaction leads to two terms 
in the two-loop part of the effective action, the second and third graphs of Fig. 4.1. 
They are the "setting sun" diagram, which is O(A^), and the "double bubble," which 
is 0(A), respectively. The Yukawa interaction leads to only one diagram in F^^^, the 
first diagram of Fig. 4.1, 

^_f_^aa'a"^WV' j ^4^^ J ^ Q ab{x , x' )T, [F,> (x , x' ) F,. {x' , x)] , 

(4.37) 

where the trace is understood to be over the spinor indices which are not shown, 
and the three-index symbol c"^'^ is defined by the n = 3 case of Eq. (2.27). Here, 
we treat the A self-interaction using the time-dependent Hartree-Fock approximation 
[79], which is equivalent to retaining the 0(A) (double-bubble) graph and dropping 
the O(A^) (setting sun) graph. We assume for the present study that the coupling A is 
sufficiently small that the O(A^) diagram is unimportant on the time scales of interest 
in the fermion production regime of the inflaton dynamics. The mean-field and gap 
equations including both the 0(A) and the O(A^) diagrams were derived for a general 
curved spacetime in Chapter 2. 



125 



4.2.3 Evolution equations for (p and G in curved spacetime 

The (bare) semiclassical field equations for the two-point function, mean field, and met- 
ric can be obtained from the CTP-2PI-CGEA by functional differentiation with respect 
Gab-> 'i>a-> 5^*^, followed by identifications of 4> and g^'^ on the two time branches, 
as shown in Eqs. (2.81)-(2.83). These equations constitute the semiclassical approxi- 
mation to the full quantum dynamics for the system described by the classical action 
(4.1). Equation (2.81) should be understood as following after time branch indices have 
been reinstated on the metric tensor in the CTP-2PI-CGEA [80]. The field equation of 
semiclassical gravity (with bare parameters) is obtained directly from Eq. (2.83), and 
given by Eq. (2.84). The right-hand side of Eq. (2.84) is the (unrenormalized) quantum 
energy-momentum tensor defined by Eq. (2.85). The energy-momentum tensor (T^jy) 
is divergent in four spacetime dimensions, and must be regularized via a covariant 
procedure [17], as discussed in Sec. 3.3.3 above. 

Making the two-loop approximation to the CTP-2PI-CGEA, where we take V2 — 
h'^r^'^\ and dropping the O(A^) diagram from the mean-field equation becomes 

(D + m'^ + (R{x) + ^4>\x) + Y^{x, x)^ 4> + hfTijFabix, x)] - h^g^T.{x) = 0, 

(4.38) 

where G{x,x) is the coincidence limit of Gabix,x'), and in terms of a function S(y) 
defined by 

S(y) = J d^x^g J dVy=7|G++(x,x')Tr.p [F++ix,y)F++{y,x')F++{x' ,x)] 

-G-+{x, x')Tr,p [F_+(x, y)F++{y, x)F+^ {x , x)] 
-G+_ (x, x')TV,p [F++ (x, y)F+_ (y, x')F_+ {x' , x)] 

(x, xOTr,p [F„+ (x, y)F+_ (y, x')F__ (x', x)] | . 

(4.39) 

Making use of the curved spacetime definitions of the scalar and spinor field Hadamard 
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kernels [17] 



G^^Hx,x') = {n\{ip^{x),ip^{x')}\n), 

fW(x,x') = (0|[^h(x),^h(x')]|0), 



and retarded propagators 



GR{x,x') = ie{x,x'){n\[ip^{x),<p»{x')]\n), 

Fr{x,x') = ie{x,x'){n\{^^ix),^^ix')}\n), 



(4.40) 
(4.41) 

(4.42) 
(4.43) 



the function S(y) can be recast in a manifestly real and causal form, 

S(y) = -2 j df^x^ y^dV^^ReTr^p (o{x,x')G^^\x' ,x)F^^\x,x') 

- Gr{x, x'^Fnix, x')) Fniy, x'yFniy, x)] , (4.44) 

from which it is clear that the integrand vanishes whenever x or x' is to the future of 
y. The "gap" equation for Gab is obtained from Eq. (2.83), 

(G-y^ix, x') = A^^ix, x') + — c'«G(x, x)5{x - x')^= 

+ n/2c-''^"c^'''^"'IV.p [Fa'b'ix,x')F,.a"ix',x)] . (4.45) 

Multiplying Eq. (4.45) through by Gab, performing a spacetime integration, and taking 
the ++ component, we obtain 



{D + m'' + ^R+-^ +—G{x,x)\g++{x,x') 



+ nf [ dx"J^ICix,x")G++(x",x') = -iS(x - x')-^, (4.46) 

J V-g 



in terms of a kernel JC{x,x") defined by 

K:{x,x') = -zTr,p [F++{x,x')F++{x',x) - F+^{x,x')F^+{x' ,x)] . (4.47) 



Making use of Eqs. (4.41) and (4.43), this kernel takes the form 
/C(x, x') = ReTr^p Fji{x,x')F^^\x' ,x) , 



(4.48) 
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which shows that the gap equation (4.46) is manifestly real and causal. As will be 
shown below in Sec. 4.3 (in a perturbative limit), the kernel IC{x,x') is dissipative, 
and it reflects back reaction from fermionic particle production induced by the time- 
dependence of the inflaton variance. The gap equation (4.46) is therefore damped for 
modes above threshold,^ and this damping is not accounted for in the IPX treatments 
of inflaton dynamics (where only the inflaton mean field is dynamical). In contrast 
to previous studies [108, 110, 112] which assumed a local equation of motion for the 
inflaton propagator, the two-loop gap equation obtained from the CTP-2PI-CGEA 
includes a nonlocal kernel, which is a generic feature of back reaction from particle 
production. As stressed above, the dissipative dynamics of the inflaton two-point 
function can be important when the inflaton variance is on the order of the square of 
the inflaton mean-field amplitude; such conditions may exist at the end of preheating. 

The set of evolution equations (4.38) for (p and (4.46) for G, is formally complete 
to two loops. Dissipation arises due to the coarse graining of the spinor degrees of 
freedom. These dynamical equations are useful for general purposes, and are valid in 
a general background spacetime. However, in order to get explicit results, one needs 
to introduce approximations, as we now do. 

4.3 Dynamics of small-amplitude inflaton oscillations 

The effective evolution equations for the inflaton mean fleld (f> and variance {ip^) de- 
rived in the previous section are useful for studying fermion production when ^q, 
the amplitude of the spatially homogeneous inflaton mean-field oscillations, is large, 
and the inflaton variance is of the same order-of-magnitude as (^o)^- discussed 
in Sec. 4.2 above, such conditions can prevail at the end of the preheating period in 
chaotic inflation with unbroken symmetry [99, 125]. Because of the dissipative kernel 
^See [84] for a similar discussion in the context of spinodal decomposition in quantum field theory. 
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lC{x,x') in the gap equation (4.46), which damps the evolution of G, and the back re- 
action terms in the mean-field equation, which damp the oscillations of 4>, eventually 
the condition, 

fk « m (4.49) 

will hold, at which point it is justifiable to follow the mean-field dynamics using the 
perturbative, IPI, coarse-grained effective action [95] . Although in principle one should 
study this process in a general curved spacetime, for simplicity we assume spatial 
homogeneity, and that the inflaton mass is much greater than the Hubble constant, 
m ^ H. While this condition alone is in general not sufficient to ensure that curved 
spacetime effects are negligible during reheating (see, for example, [125], where cosmic 
expansion does affect preheating dynamics even though m':$> H), with the additional 
assumption of condition (4.49) it is reasonable to neglect the effect of cosmic expansion 
in the spinor propagators [111, 112]. In this and the following section, we also neglect 
the self-coupling A, because for the case of unbroken symmetry, the lowest-order A- 
dependent contribution to the perturbative inflaton self-energy is 0(^2) [108], and we 
are only concerned with one-loop dynamics in this section. 

Let us therefore specialize to Minkowski space, and implement a perturbative ex- 
pansion of the CTP effective action in powers of the mean field 0. This formally entails 
a solution of the gap equation (4.46) for G, a back-substitution of the solution into the 
CTP-2PI coarse-grained effective action, and a subsequent expansion of this expression 
(now a functional of (f) only) in powers of (f). The resulting perturbative expansion for 
the effective action contains only free-field propagators. For consistency, one should 
use an initial density matrix for the spinor degrees of freedom which corresponds to the 
end-state particle occupation numbers of the nonperturbative dynamics of Sec. 4.2. 
For simplicity, however, we assume the initial quantum state for the spinor field is the 
vacuum state. Hereafter, F^b denotes the free-field, Minkowski-space, vacuum spinor 
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CTP propagator, whose components are given by [63, 107, 117] 

F^Ux x') = [ ^e-'P^--')Si±^ (4 50) 

F-M = - I ^e-(-')4^±!!!L, (4.51) 

F-+{x,x') = J -^e-^P^---'hn{^ + m)d{p^ - fx^)e{p% (4.52) 

F+-ix,x') = J ^e-^f(---')27r(i, + m)5(/ -/z2)e(-/). (4.53) 

The ++ and propagators admit a representation in terms of a time-ordering 

function 9{x,y) = 6{x^ — %p), 

F++(x,x') = e{x,x')F_^{x,x') + 9{x' ,x)F^-{x,x'), (4.54) 
F__{x,x') = e{x,x')F+_{x,x') + 9{x' , x)F__|_(x, x'). (4.55) 

The CTP effective action can be expanded in powers of and we find 

r[0] = <S'^[0] - ylndet(^«^)-^ - mndetFab + Ti[4>], (4.56) 

The CTP effective action can be expanded in powers of and we find 

T[$\ = S'^[(i)] - ^lndet(^'''')-^ - imndetF„6 + Ti[4>], (4.57) 

where the kernel A is defined by 

iA''\x, x') = -c''\n^ + m^)6{x - x'), (4.58) 

and Fi is defined as —ih times the sum of all one-particle-irreducible diagrams con- 
structed with lines given by hA~^ and hFab, and vertices given by S^[(l),tlj,'ip]/h and 
S^[(p,'tp,'ip]/h. Because the free-field propagators .4"^ and Fab do not depend on (p, 
the log(det) do not contribute to the variation of r[0] with respect to (f>, and therefore, 
they can be dropped. The functional ri[0] can be expanded in powers of h, 

oo 

ri[<^] = 5;n'r«[0], (4.59) 
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where the term r(')[0] is the sum of all IPI /-loop graphs. Order by order in the loop 
expansion and the coupling constant, the OTP IPI effective action must satisfy the 
unitarity condition 

which has been verified to two-loop order in the case of scalar field theory [66]. 
The one-loop term in the loop expansion of the CTP effective action, P'^^^f^], can be 
further expanded in powers of 



(4.61) 



n=l 



which corresponds to the usual amplitude expansion of the CTP effective action [112]. 
Figure 4.2 shows the diagrammatic expansion of P^^^ , where solid lines represent the 



.(1) 




Figure 4.2: Diagrammatic expansion for p(^), the one-loop part of the CTP-lPI coarse- 
grained effective action. 

spinor propagator F (as defined in Sec. 4.3), and dotted lines represent multiplication 
by the scalar mean field cf). Each vertex carries spacetime (x) and CTP (-I-, — ) labels. 
The terms P2n['^] generally divergent, but since the theory is renormalizable in 
the standard "in-out" formulation, it is renormalizable in the closed-time-path, "in-in" 
formulation [66,67]. 



131 



4.3.1 One-loop perturbative effective action at 0{p) 

The 0{p) term in the expansion of the one-loop CTP effective action, which is the 



first term in Fig. 4.2, takes the form^ 

r?^M = -|c«''^c«'^'<='y" d^xd^x'Ux)^a'{^')T,^, [F,^{x,x')F,,,{x',x)\ 
Making use of sum and difference variables 

nx) = \[M^)+L{x') 

A{x) = 4>+{x)-^M), 
the functional f!^^ [^] can be recast in the form 

d^xd^x' 

in terms of manifestly real kernels V2{x,x') and Af2{x,x') defined by 



(4.62) 



(4.63) 
(4.64) 



J:{x)A{x')V2{x,x') + ^A{x)A{x')J\f2{x,x') 



(4.65) 



'D2{x,x') = lmTr^p[F++{x,x')F++{x',x) + F+^{x,x')F^+{x' ,x)], 
M2{x,x') = - ReTrjF++{x,x')F++{x',x)]. 



(4.66) 
(4.67) 



Only the kernel T>2{x, x') contributes to the mean-field equation of motion. The kernel 
Af2ix,x') constitutes a correlator for noise, and will be discussed in Sec. 4.4. The 
unitarity condition (4.60) requires that the sum of diagrams proportional to 'E{x)T,{x') 
vanish identically. With the definitions of the retarded spinor propagator, Eq. (4.41), 
and the spinor Hadamard kernel, Eq. (4.43), which in Minkowski space take the form 



Fji{x,x') = i9{x,x')[F-^{x,x') — F^ {x,x')], 

F^^\x,x) = F_+{x,x') + F+_{x,x'), 

the kernel V2{x,x') can be written in a manifestly causal form, 

V2{x,x') = ^ReTrjFR{x,x')F^^Hx',x)]. 



(4.68) 
(4.69) 



(4.70) 



*Note that there are no nonzero graphs with an odd number of vertices in this model. 
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Using Eq. (4.70), one can compute T)2{x, x') in an arbitrary curved background space- 
time. It should be noted that 'D2{x,x') is just the lowest-order term in the series 
expansion of K{x,x') [defined in Eq. (4.48)] in powers of the coupling constant /. The 
appearance of the retarded propagator in Eq. (4.70) guarantees that the contribution 
of r2^'* to the mean-field equation of motion is causal. 

Let us now evaluate V2{x,x') using dimensional regularization and the modified 
minimal subtraction (MS) renormalization prescription [153,182,194]. Dimensional 
regularization requires changing the coupling constant so that the interaction has 
the correct dimensions in n spacetime dimensions, 

f^f^"^, (4.71) 

where we have introduced a parameter A, the renormalization scale, which has di- 
mensions of mass. By Lorentz invariance and causality, the product of Feynman 
propagators can be written in terms of an amplitude A2 [153], 

TV.,[F++(x,x')F++(x',x)] =i J ^e-"'^---")A2{e + te), (4.72) 

and with this choice of renormalization prescription, the amplitude A2{k'^) takes the 
form 

A2ik') = £ daE{a-k')log (^^^^ , (4-73) 

where E{a;k'^) is defined by^ 

E{a;k^) = n"^ -a{l-a)k'^. (4.74) 

Note that in Eq. (4.73), the a integration shows up via the Feynman identity [153] 



— ^—pr- = (-^ - 1)' / dai-- - [ daNS{ai -\ h aiv-i) ["iCi -I h a^vCAr] ^ 

(^l---(^N Jo Jo 

(4.75) 



®The notation E{a; k^) used here should not be confused with E{k), the complete elliptic integral 
of second kind. 
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The ie appearing in Eq. (4.72) ensures that the ampUtude A2 is evaluated on the 
physical sheet [153,195]. The logarithm in Eq. (4.73) has a negative real argument 
when the two conditions > Afi^ and |2a — 1| < a/1 — Ajj? /hP' are both satisfied. 
When |2q; — l| < ^1 — Ajj.'^/k'^, the amplitude ^2(^^) has a branch cut (considered 
as an analytically continued function of k^) for (/c*^)^ > fc^ + Afp. The discontinuity 
across the branch cut is related to the "cut" version of the diagram [the second term 
in Eq. (4.66)] via the Cutosky rules [172, 195-199], 

TrjF+_{x,x')F^+{x',x)] = -i J il^e-M--a:')Disc[^2(A:2)]^(/c°). (4.76) 

From Eqs. (4.72), (4.73), and (4.76), it is straightforward to obtain an expression for 
the dissipation kernel. 



^2(x,x') = i/(0e-^(-')^^aE(a;/c 



log 



\E(a:k 



2) 



A2 / 



, (4.77) 



where we have now taken the limit e — > 0+. One can verify by inspection that this 
kernel is real. However, the second term in Eq. (4.77) breaks time-reversal invariance 
and leads to dissipative mean field dynamics. The one-loop Fourier-transformed mean- 
field equation is (dropping the caret from (f>) 

'\E{a-k'^) 



e - m'^ + ik^j2ik) - t daE{a-k^)\og 



m = -J{k), (4.78) 



A2 

where 72 (A;) is the dissipation function, defined as —ihplk^ times the Fourier trans- 
form of the second term in Eq. (4.77), 

Uk) = ¥lm[P2(A;)] = (l - ^V^%(/c^ - V). (4.79) 



k^ ' " Sir V J 
The one-loop 0{p) dissipation kernel agrees with previous calculations of the proba- 
bility to produce a fermion particle pair of momentum A; [91, 94, 198, 200]. In Eq. (4.78), 

J{k) is an external c-number source. The imaginary term ik^^2^ in Eq- (4.78) breaks 
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time-reversal invariance and acts as a /c-dependent dissipative force in the mean field 
equation. The 6 function enforces the energy threshold for the virtual fermion pair in 
the one-loop 0{p) diagram to go on-shell. The dissipative mean-field equation (4.78) 
is essentially the linear-response approximation to the effective inflaton dynamics. It 
should be noted that the dissipation kernel T>2{x,x') is nonlocal, in contrast with the 
local friction- type dissipation assumed in earlier studies of post-inflation reheating [10]. 
However, in the limit ji^ 0, the dissipation kernel does become time-local, as there 
is no longer a length scale in the expression for 'D2{x,x') which could define a time 
scale for nonlocal dissipation [68]. 



4.3.2 One-loop perturbative effective action at 

The 0(/^) term in the one-loop CTP effective action consists of the "square" diagram, 
which is the second term in Fig. 4.2, 



r«[<^] = y^bc^a'b'c'^def^d'e'f j ^^d" x' y' y 



Tr,p [Fw{x, x')Fc'f{x', y')Fe>e{y', y)Ffc{y, x)] 



X k{x)k'{x')4>d{y)4>d'{y') 



(4.80) 



Expanding out the contracted CTP indices, we obtain 



■.(l)r 



j d^^xd^x'd^d^y 
ix)4>+ {x')4>+ {y)4>+ (y')Tr.p { F++ (x, x')F++ (x', y')F++ {y\ y)F++ (y, x) } 

+ 4>- {x)4>- (x')0- {y)^- (y')iv.p (i^— (x, x')i^— (x', y')^— {y', y)F— (y, ^) } 

-40+(x)0_(x')<^-(y')'^-(?/)Tr.p{F+_(x,x')F__(x',y')i^— 

-40_(x)0+(xO<^+(2/')'^+(y)^.p{^-+(^,^O^++(^',yO^++(y',2/)^+-(y,^)} 

+ 40+(x)0+(xO<^-(yO0-(2/)'I^.p{^++(a^,aj')^+-(^',y')^--(2/',y)^-+(2/,^)} 
+ 20+ (x)0_ (x')<^+ (y')0- (y)TV.p (x, xOF_+ (x', y')F+_ {y' , y)F_+ {y, x) ] 

(4.81) 
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When r4^^ is expressed in terms of A and S [defined in Eqs. (4.63) and (4.64)], only 
those terms with one factor of A and three factors of S contribute to the mean field 
equation of motion. As a consequence of the unitarity condition (4.60), the sum of 
terms proportional to four factors of S must vanish. Keeping only those terms in the 
effective action which contribute to the mean field equation or are quadratic in A, we 
find 



+ - A(x) A(x')S(y')S(y)M(x, x',y', y) 



in terms of a kernel V^i^x, x' ,y' , y) defined by 



Vi{x,x' ,y\y) = -ImTr, 



F++{x,x 



+ F++{x,x')F+^ 
+ F+_(x,x')F— 
+ F+_(x,x')i^-+ 

-F+_{x,x')F^+ 

-F++(x,x')F+- 
-F++{x,x')F++ 



F++ix',y' 

x',y')F^- 

x',y')F-+ 
x',y')F+_ 
x',y')F++ 
x',y')F^+ 

x',y')F+- 



F+_{x,x')F_ 



and a "noise" kernel M^lx , x' , y' , y) defined by 



F++iy',y)F++{y,x) 

y',y)F-+{y, x) 
y\y)F++{y,x) 
y',y)F-+{y, x) 
y',y)F++{y,x) 
y',y)F++iy,x) 
y',y)F-+{y, x) 

x',y')F_.iy',y)F_+iy,x) 



Af4{x,x',y',y) = ReTr,. 



F++{x,x')F++{x',y')F++{y',y)F++{y,x) 



+ F++ {x, x')F+^ (x', y')^— (?/', y)F-+ (y, x) 
-F+_{x, x')F__{x', y')F_+{y', y)F++{y, x) 
- F+_ (x, x')F_+ (x', y')F+^ (y' , y)F^+ (y, x) 
+ F+_(x, x')F_+{x', y')F++{y', y)F++(y, x) 



(4.82) 



, (4.83) 
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- F++{x, x')F+_{x', y')F_+{y', y)F++{y, x 

- F++{x, x')F++{x', y')F+_{y', y)F_+iy, x 
+ F+_ (x, x')F__ {x\ y')F__ {y' , y)F_+ {y, x 
+ {F++ (x, y')F++ {y', x')F++ {x' , y)F++ {y, x) 

- F+_ (x, y')F__ {y', x')F_+{x' , y)F++{y, x 

- F++ix, y')F+^iy', y)F^+iy, x 
+ F+_ (x, y')F_+ iy', x')F+_ {x' , y)F_+ (y, x 

- F+_{x, y')F_+{y', x')F++{x', y)F++{y, x 
+ F++(x, y')F^-{y\ x')F-+(x', x 

- F++ (x, y')F++ {y', x')F+_ (x', y)F_+ {y, x 

+ F+_ (x, y')i^— (y', ^')^— y)F-+iy, x)}/2 



, (4.84) 



The noise kernel Ma does not contribute to the mean field equation of motion. There 
are, of course, terms in r4^^[0] which are higher order in A, for example, O(A^), but 
in passing over to a stochastic equation for (f> in Sec. 4.4, we will be assuming that A 
is small, so that higher-order terms in powers of A can be ignored. Such terms will in 
general contribute to non-Gaussian noise, which will be studied in an upcoming paper 
[201]. 

Let us evaluate the first term of Eq. (4.83), which consists of only Feynman propa- 
gators. The term is logarithmically divergent, and as in Sec. 4.3.1, we use dimensional 
continuation and the modified minimal subtraction (MS) renormalization scheme. Be- 
cause we are only interested in deriving the dissipative terms in the mean-field equation 
coming from this diagram, and because we are assuming m 3> we include only the 
one-loop logarithm. We find 



Tr 



F++{x, x')F++{x', y')F++{y', y)F++{y, x) 



log only 
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J [ZttP [Zirp [Zirp 



where the amplitude A4{ki,k2, ks) is defined by 



^4(^1,^2,^:3 



27r2 



daida2da3 log 



-£'4(0:1, «2, 0:3; ^1,^2,^3) 
A2 



(4.86) 



n2 



in terms of a function E4 defined by 

E4{ai,a2,a3;ki,k2,k3) = (1 - q;i)A;i + (1 - ai - 02)^:2 + (1 - oi - 0:2 - q;3)A;3 

- (1 - ai)kl - (1 - ai - a2)(fei + 2fei • ^2) 

- (1 - ai - 0:2 - a3)(2A;i • ^3 + 2A;2 • ^3 + A;|) + /x^. 

(4.87) 

As in Sec. 4.3.1, the cut diagrams in Eq. (4.83) are related to the discontinuities in 
-£4(01, 02, 03; ki, k2, ks) via the Cutosky rules. The details are shown in the appendix. 
We can then express the dissipation kernel V4{x,x' ,y' ,y) as a Fourier transform over 
external momenta, 



V4{x,x',y',y) 



d'^ki d'^k2 d'^ks i\_(ki+k2+k3)-x+ki-x'+k2-y'+k3 
(27r)4 (27r)4 (27r)4 



y^v^ikiMM), 

(4.88) 



in terms of a function 'D4{ki,k2,ks) defined by 



'D4{ki,k2,k3) = — 



1 /" , , , 1 /|-E4(ai,a2,Q;3; A:i, A:2,fc3)| 
— / daio(a2aa3log I 



+ isgn(A;^ + ki)h[{k2 + kaf] + isgn(/c? + k^)h[{ki + fcs)^] 
+ isgn(A;? + k^ + k^)h[k^] + isgn{k^)h[{ki + k2 + ksf] 

+ isgn{k^s)h[kl] + isgnik'Dhikl] - iH{ki,k2, k^) 



(4.89) 
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and the functions h{s) and H{ki, k2, k^) are defined by 



Ms) = \/i-^^(s- V), 



(4.90) 



H{ki,k2,k3) = / {9[-E4^{ai,a2,a3;ki,k2,k3)] 

J ai,a2,a3>0 



(4.91) 

Equation (4.89) leads to the foUowing mean-field equation at 0{f"^), 



f d^q dH 



47r2 



/ daE{a]k'^)\og 
Jo 



A2 



2n 



J daida2da3 log ^- — ^ — — — — 



+ sgn(g° + f)h[{q + Z)2] + sgn(fe'^)/i(g2) 
+ sgn(A;° - l^)h[{k - if] + sgn(g°)/i(A:^) 
+ sgn(A;° - g° - /°)/i(/2) 

+ sgn{f )h[{k -q- If] 



+ H{l + q-k,-q,-l) 



-J{k). 



(4.92) 



The presence of terms of the form isgn{p^)h{p'^ ) in Eq. (4.92) clearly signifies dissipative 
dynamics. The ^-function in Eq. (4.90) enforces the energy threshold for the virtual 
fermion quanta created at a particular vertex to go on-shell. Comparing Eq. (4.92) 
and Eq. (4.78), and assuming spatial homogeneity, we see that the 0{f^) dissipation 
kernel must be taken into account whenever the condition (4.49) fails to hold for the 
solution (p{t) to Eq. (4.78). 

At the end of the regime of parametric resonance in chaotic inflaton, i.e. the 
"preheating" regime, the inflaton mean field may oscillate with an amplitude as large 
as ^ rn/g^^, where g^^^ is the coupling to another scalar field typically on the order 
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of 10 ^ [166]. Condition (4.49) would then be violated if / > gtj,^. In this case it would 
be necessary, at a minimum, to take into account higher order terms (such as P4) in 



point where Eq. (4.49) is satisfied. 

4.4 Noise kernel and stochastic inflaton dynamics 

Although the kernels N'2{x, x') and N^ix, x' ,y' , y) do not contribute to the mean field 
equation, i.e., the equation of motion for 0, they contain information about stochas- 
ticity in a quasi-classical description of the effective dynamics of the inflaton field 
[53, 54, 56, 74-76, 188, 202, 203]. In this section, we study the effect of stochasticity on 
the dynamics of the inflaton mean field, within the perturbative framework established 
above. 

4.4.1 Langevin equation and fluctuation-dissipation relation at 0{p) 

In this section we show how to obtain a classical stochastic equation for the infiaton 
field from the 0{f^) perturbative CTP effective action. Prom Eq. (4.65), it follows 
that the 0{f) one-loop perturbative CTP effective action has the form 



the mean-field equation, until such time as the amplitude i?!>o(i) has decreased to the 




where for simplicity we have defined 



l'2ix,x') = hf 



'^M2{X,X'), 



(4.94) 



H2{x,x') = hp'D2{x,x'). 



(4.95) 
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In order to extract the stochastic noise arising from the kernel M2{x,x'), we use the 
Gaussian identity [151] 

1 



exp 



2h 

D^2 exp 



d*xd^x'A{x)A{x')iy2{x,x' 



, (4.96) 



where iV is a normahzation factor which does not depend on A, and ^2 is a c-number 
functional integration variable. Following [202] , we now define a functional 

1 



^[6] = A^exp 



h 



d^xd'^x'C2{x)i^2Hx,x')C2{x') 



(4.97) 



and it follows from Eq. (4.96) that P[C2] is normalized in the sense of 

J DC2PM = 1. (4.98) 
Using Eq. (4.96), we can rewrite the 0{p) one-loop OTP effective action, Eq. (4.93), 



as 



T[4>] = -iMog / D^2PMexp 



-(^5'^[(^]+ J d^xdVE(x)A(x')M2(a^,a;') 

+ J d^x^2ix)A{x)^ . (4.99) 

This suggests defining a new effective action which depends on both ^2 and (j)^ (drop- 
ping the carat from 0), 

T[(j),^2]=S'^[(t)]+ j d'^xd^x'^{x)A{x')n2{x,x') + J d'^x^2{x)A{x). (4.100) 

Let us define a type of ensemble average 

{{A}) = I L>6^'[6]^(6), (4.101) 

and note that Eqs. (4.97) and (4.101) imply that 



mx))) = 0, 

{{Ux)Hx'))) = nMx,x'). 



(4.102) 
(4.103) 
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Clearly then, 



-amm 



4>+=<t>- 



-rf 



5+ 



<l>+=4>-=<t> 



Taking the variation of T[^, ^2] with respect to and setting 0+ 
(after a Fourier transform) 



(4.104) 
(f), we obtain 



A2 



where ^2(fc) is defined by 



4>{k)+Uk) = -J{k), 
(4.105) 

(4.106) 



We now interpret Eq. (4.105) as a Langevin equation with stochastic force ^2- The in- 
fiaton Fourier mode <^ appearing in Eq. (4.105) should be viewed as a c-number stochas- 
tic variable, and the presence of the stochastic force ^2 indicates spontaneous breaking 
of spatial translation invariance by a Gaussian (but not white) noise source ^2 [56]. 
Moreover, this stochastic equation obeys a zero-temperature fluctuation-dissipation 
relation, as we now show. First, let us calculate the one-loop 0{p) noise kernel, 
A/2 (a;, a;') [defined in Eq. (4.67) above], using dimensional regularization and modified 
minimal subtraction. 



in terms of the Fourier-transformed noise kernel 



-ik{x—x')~ 



(4.107) 



Hk) = ^k^{i-'-fyeik^-4, 



k^ 



2x 3 



(4.108) 



The noise kernel V2{x,x') is colored; colored noise has been observed in other inter- 
acting field theories [53,56,202]. By inspection of Eqs. (4.79) and (4.108), it follows 
that 

\k%2{k) = V2{k), (4.109) 
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which leads to the zero-temperature fluctuation-dissipation relation [202], 

/oo 
dt'K{t-t'h2it',k), (4.110) 
-oo 

in terms of the distribution- valued kernel K{t) defined by 

roo J 

K{t) = / —ucos{ut), (4.111) 

and the spatially Fourier-transformed dissipation function and noise kernel, 

Mt,k) = j_ ^e-^^''Uk), (4.112) 
I2{t,k) = y^^_e-^°*72(fc). (4.113) 

This shows the physical significance of the noise kernel i'2{x, x') in an effective descrip- 
tion of the dynamics of the scalar mean field. 

4.4.2 Langevin equation and fluctuation-dissipation relation at 0(/^) 

In this section we consider the 0{f^) one-loop noise kernel, 7V4. The non-normal- 
threshold singularities in A4 lead to a noise kernel which depends on S, which is known 
to lead to ambiguities in the resulting Langevin equation [204, 205] . The meaning and 
interpretation of the non-normal-threshold parts of A/4 and P4 will be the subject of 
a future study [201]. Here, we focus on the effect of the normal-threshold singularities 
of ^4, which for the noise kernel, contribute a term 

dWx'A(x)A(x')S(x)S(x')i^4(a;,x') (4.114) 

to the CTP effective action, where the kernel i'4{x,x') is defined by 

u,{x,x') = J ^e-^-(— ')M,^), (4.115) 

and the function h{s) was defined in Eq. (4.90) above. The normal-threshold singu- 
larities of the dissipation kernel, P4, [the second and third terms of Eq. (4.89)], lead 
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to the following contribution to the CTP effective action, 

J dWx'A(a;)S(x) [E(x')]%4(a;,x'), (4.116) 
where the kernel fj,4{x,x') is defined by 

M4(x,x') = J ^e-^«-(-'-)sgn(g'^)Mg^). (4.117) 
With the definitions 

/X4(a:, x') = il ^e-^^<^'-^k'UQ), (4.118) 

'^4ix,x') = J -0^e-"^<-'--^h{q), (4.119) 

it follows immediately that the normal-threshold parts of V4 and A/4 obey a fluctuation- 
dissipation relation identical in form to Eq. (4.109), 

\q°\%{q) = Ml)- (4.120) 

Making use of Eq. (4.96), the 0{f^) effective action (including only normal-threshold 
contributions) can be written in the form 

r[<^,6,?4] =5*H+ J d'^xd'^x'A{x)I:{x')^^2{x,x') + J d^xUx)^{x) 

+ J cZ^xdVA(x)S(x) [S(x')]V4(a;,x') + J d^x^4{x)A{x)^{x), (4.121) 

where the stochastic noise source ^4 satisfies the conditions 

{{U{x)))=0 (4.122) 

{{Ux)Ux'))) = hu4{x,x'). (4.123) 
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Taking the functional derivative of Eq. (4.121) and making the usual identification, 
we obtain a Langevin equation with an additive noise ^2 and a multiplicative noise ^4, 

p2 /■! /IJ?/'^,. 



q -m +<72(g)-^^ daE{a;q ) log I ) 

/d^ k dH 



(4.124) 

where d^a = daida2da3, and Q{1, q, k) is defined as 

Q{l,q,k) = trj4{l} - d-'alog I 1 . (4.125) 

The stochastic force ^4 is clearly seen to contribute multiplicatively to the Langevin 
equation for cj). 

4.4.3 Homogeneous mean field dynamics at 0{p) 

To make connection with post-inflationary reheating, it is customary to assume that 
the mean field (j) is spatially homogeneous [90,95,166]. In this case, the Langevin 
equation (4.105) takes the form 

[lo^ -m^ + iL0/3{uj) + ri{u)] ^{uj) + ^(a;) = -J(a;), (4.126) 

where we have defined 



"(•^'^ 8^Ri'-^J ('(-'-V), (4.127) 
,H =-^jf'a„£(„;.^)l„g(S^) . (4.128) 

The total energy dissipated to the fermion field over the history of the dynamical 
evolution of the mean field is given [at 0{p)] by 

£ = - dtFS)^, (4.129) 
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where the friction force Fy(t) is the Fourier transform of if3{uj)iO(l){io). After a bit of 
Fourier algebra, we obtain an expression for the ensemble-averaged, total dissipated 
energy, 

f2 foo / A,,2\ 3/2 I f/ \|2 



It is straightforward to compute the variance in the total dissipated energy. Making 
use of Eq. (4.103), we find 

(4.131) 

where the function is defined by 



j-^^Lo'^-Ati^ / 42 \3/2 

Following [91], we assume that the infiaton field is held fixed via an external, constant 
c-number source J for t < 0, and that the source is removed for t > 0, 

J(t) = Je{-t). (4.133) 

Setting ^ = 1, assuming that m ^ fj,, and expanding to lowest order in /, we obtain 
for the ensemble averaged dissipated energy. 



{{S)) = ^S^. (4.134) 



Let us now compute the variance in the total dissipated energy. Performing a regu- 
larization via dimensional continuation, we obtain 

where (5 is a constant of order unity defined by (5^ = |119/60— 7em— log(47rm2/A2)|, and 
7em is the Euler-Mascheroni constant, 0.5772. Taking the ratio of the square root of 
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Eq. (4.135) and Eq. (4.134), we obtain the relative strength of the RMS fluctuations 
in the total dissipated energy density, S^^^, 

The parameter J is related to the initial inflaton amplitude ^o(^o) by --^ = 4>o{to)m'^ /2, 



which leads to 



^^"^^ 0.390-;^. (4.137) 



{{£)) Vl5TT^o{to) Mto) 

The fundamental assumption which justified the pcrturbativc expansion in /, Eq. (4.49), 
is seen to be independent of Eq. (4.137). Therefore, the ratio i'rms/^' is not required to 
be small by consistency with perturbation theory. As the initial inflaton amplitude 0o 
is made larger, the relative strength of the rms fluctuations of £ is seen to decrease, in 
accordance with the correspondence principle. It has been shown that the fluctuations 
in the total dissipated energy density are related to the fluctuations in the occupation 
numbers of modes [76]. 

Let us now examine whether the rms fluctuations in the total dissipated energy, as 
given by Eq. (4.137), is significant, given a reasonable value for the inflaton amplitude 
at the end of the preheating regime (the period of parametric resonance-induced par- 
ticle production). In chaotic inflaton with a scalar held % coupled to the inflaton field 
via a coupling constant 5^^, the typical inflaton amplitude at the end of the preheating 
regime is on the order of m/gcj,-^ [166]. In this case, we would find ^rms/((^)) — /9</>x5 
from which it is clear that fluctuations in the total dissipated energy are not signif- 
icant relative to the ensemble-averaged total dissipated energy, and therefore should 
not appreciably affect the reheating temperature. However, in new inflation scenarios 
where the inflaton amplitude 0o can be on the order of m at the onset of reheating, 
the ratio m/0o can be of order unity [10]. In this case, the ratio S^^s/{{£)) ~ /, which 
may not be a negligible effect. 

Although as shown above, stochasticity does not dramatically affect the total en- 
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ergy dissipated via fermion production in chaotic inflation, we may inquire whether 
the noise term in the Langevin equation for the inflaton zero-mode, Eq. (4.126), may 
nonetheless be non-negligible during the reheating period. Let us compute the rms 
fluctuations in the inflaton zero- mode, 4>tuis- Following methods described above, we 
find that the rms fluctuations of the inflaton zero mode are given, to 0(/^), by 

^ -.a, (4.138) 



where ct^ = |61/30 ~7bm — log(47rm^/A^)|. Equation (4.138) is seen to be independent 
of the inflaton zero-mode amplitude 0o- 

In order to determine the relative importance of fluctuations in the inflaton zero- 
mode amplitude 4>q during and at the end of the reheating period, we must introduce 
curved spacetime arguments. This is because the end of the reheating period is deter- 
mined by the time tend at which the Hubble constant becomes of the order of 3/3(m) for 
the case of A = being discussed in this section [10, 95]. Starting with the semiclas- 
sical Einstein equation (2.84) for spatially flat Friedmann-Robertson- Walker (FRW) 
cosmology, setting 6 = c = Ac = (following arguments similar to those of Sec. Ill D 
of Ref. [125]), retaining only the inflaton zero mode as the dynamical degree of freedom 
(for consistency with FRW), and retaining both the ^ field energy density and the 
classical, stochastic energy density of the inflaton zero mode, we have 

^' = ^ = ^fp(?)+P^V (4.139) 



p 



where a is the scale factor, the dot denotes a derivative with respect to cosmic time, 

and p{(j) ) is the energy density as a function of the time-average (over one period of 
oscillation) of 4> , which is given by the virial theorem, 

p{4>) ~ m^4>^ = ^m'^i^of. (4.140) 

Making use of Eqs. (4.126), (4.133), and (4.139), we obtain an approximate expression 
for the (ensemble-averaged) energy density of the inflaton zero-mode at the end of the 
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reheating period^ (to lowest order in /), 

p{t..a) - (g^j^se ' (4-141) 

where e is the base of the natural logarithm. Note that this expression is independent 
of the initial inflaton amplitude [166]. Equation (4.141) allows us to solve for the value 
of 00 the end of reheating. We find 

Mte..) - ^^^f - (4-142) 

The rms fluctuations in the inflaton zero mode, (p^^^, can only play a role in the 

inflaton zero-mode dynamics during reheating if the ratio (p^^^/cpQ is not small relative 

to higher-order [e.g., 0(/^)] processes which we are neglecting. In light of the minimum 

inflaton zero- mode amplitude attained during reheating, Eq. (4.142), we find that the 

ratio of fluctuations in the inflaton zero-mode to the zero-mode amplitude is given by 

Sam He' ^ m , . ^ .^^ 

'— ~ 2.37— — . (4.143) 



<^o(iend) /^p V 45 /Mp 
We estimate the ratio of the mean-squared inflaton amplitude fluctuations to the shift 

in the inflaton mass to be 

™^ 0.01. (4.144) 



ri(m) 

If, prior to the end of reheating at tend, 0rms/'/'o(^) becomes larger than higher-order 
terms which are neglected in our perturbative expansion, then fluctuations in the 
inflaton zero mode are a non-negligible effect. This will happen when (f^^^^/cpoit^nd), 
given by Eq. (4.143), is not <C 1. This shows that stochasticity must be taken into 
account in the dynamics of the inflaton zero mode, during the late stages of the 
rdioating period. 

®We wish to emphasize, however, that this expression does not take into account the regime of 
nonperturbative dynamics discussed in Sec. 4.2, and therefore should not be expected to yield a correct 
reheating temperature in a realistic inflationary scenario. However, it suffices for the present discussion 
of rms fluctuations of the inflaton amplitude, where we assume an idealized case similar to Eq. (106) 
of Ref. [95]. 
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4.5 Summary 

In this chapter, we present the results of a study of (unbroken-symmetry) inflaton 
dynamics during the late stages of reheating, which is dominated by fermion particle 
production to a light spinor field coupled to the inflaton field via a Yukawa coupling. 
We derived coupled nonperturbative equations for the inflaton mean field and two- 
point function, in a general curved spacetime, and showed that, in addition to the 
dissipativc mean-field equation, the gap equation for the two-point function is also 
dissipative, due to fermion particle production. Simultaneous evolution of the inflaton 
mean-field and two-point function is necessary for correctly following the inflaton dy- 
namics after the end of the preheating period, because the large value of the variance 
invalidates use of the ordinary perturbative, IPI effective action. 

We also derived the dissipation and noise kernels for the small-amplitude dynamics 
of the inflaton field, valid in the late stages of reheating when the inflaton mean-field 
amplitude is very small. The 0{p) noise and dissipation kernels, as well as the 
normal-threshold parts of the 0{f^) noise and dissipation kernels, are shown to obey 
a zero-temperature fluctuation-dissipation relation. With the noise and dissipation 
kernels, a Langevin equation for the inflaton zero mode is derived, and it is shown 
that the noise leads to a variance for the inflaton amplitude which is non-negligible 
before the end of reheating. 
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CHAPTER 5 



Correlation entropy in effectively open systems 
5.1 Introduction 

As emphasized in Chapter 3, an important issue within inflationary cosmology is how, 
and to what temperature, the Universe reheats following the period of profuse par- 
ticle production. Knowledge of this temperature is important for consistency of the 
inflationary Universe picture with the standard big bang cosmology. Because of the 
self-coupling of the inflaton fleld and its coupling to other quantum flelds, it is expected 
that the inflaton fleld and those coupled to it will eventually come to local thermal 
equilibrium in the expanding background spacetime, with the usual equilibrium equa- 
tion of state depending on the masses of the particle species relative to the temperature 
of the plasma [10]. The quantitative description of how the system reaches local ther- 
mal equilibrium, and at what temperature, is known as the thermalization problem. 
Though there have been some previous studies of this problem, most either employed 
classical arguments in an essential way [116], or utilized initial conditions which are 
not appropriate to the end of preheating in realistic inflation scenarios [206]. In light 
of the recent, newer understanding of the role that parametric resonance effects play 
during reheating, the thermalization problem is at present unsolved [207] . Because of 
the tremendous variety of different inflationary models, and the fact that the details of 
preheating dynamics of the inflaton fleld is severely model-dependent, any discussion 
of the thermalization period will of necessity be rather general.^ However, it is useful 

^In fact, in some cases, the preheating and thermahzation stages do not separate at all, and particle 
production and thermalization must be considered simultaneously [147], though in this discussion we 
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to explicate some of the challenges, both conceptual and technical, inherent in the 
thermalization problem. 

First and most important, it is not clear that a separation of macroscopic and 
microscopic time scales exists at the end of the preheating period [99]. In the two- 
particle- irreducible language of Chapter 3, a separation of macroscopic and microscopic 
time scales would require that the initial conditions be such that one can assume a 
quasilocal expansion for the time dependence of the effective mass. In fact, preliminary 
evidence suggests that for many realistic inflation scenarios, such a separation does 
not exist at the end of the preheating stage, as can be seen in Chapter 3, where the 
time scale for variations of the effective mass is on the order of the microscopic time 
scale for the theory. Such a separation constitutes the basis for passing over from the 
full, closed-system dynamics of quantum fields (given by the functional Schrodinger 
or quantum Liouville equation) to a quantum kinetic theory description in which one 
obtains separate equations for quantum field-theoretic processes on microscopic time 
scales, and relaxation phenomena on macroscopic time scales [68]. By quantum kinetic 
field theory, we are referring to the hierarchy of coupled equations for the relativistic 
Wigner function and its higher-correlation analogues, which are obtained by a Fourier 
transform of the relative coordinates in the Schwinger-Dyson equations for correlation 
functions (or alternatively, in the master effective action whose variation yields the 
Schwinger-Dyson equations). This is a quantum analogue of the BBGKY hierarchy 
[208], expressed in a representation convenient for distinguishing between microscopic 
(quantum field-theoretic) and macroscopic (transport and relaxation) phenomena. As 
such, it does not require near-equilibrium conditions, and in fact, is applicable for 
a rather general moment expansion of the initial density matrix [68]. It should be 
pointed out that in order to identify the relativistic Wigner function with a distribution 

will assume that such a separation exists, and will think of the end state of the preheating stage as 
constituting the initial conditions for thermalization. 
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function for quasiparticles, one must show that the density matrix has decohered, and 
this is neither guaranteed nor required by the existence of a separation of macroscopic 
and microscopic time scales [209]. 

Let us briefly comment on the important relation between quantum kinetic field 
theory in its full generality, and an effective rclativistic Boltzmann description of re- 
laxation phenomena for the one-particle distribution function of quasiparticles. In 
nonequilibrium statistical mechanics, as is well known [208, 210], the act of truncating 
the BBGKY hierarchy does not in itself lead to irreversibility and an iJ-theorem. One 
must further perform a type of coarse graining of the truncated, coupled equations 
for n-particle distribution functions. For example, if one truncates the hierarchy to 
include only the one-particle and two-particle distribution functions, it is the sub- 
sequent assumption that the two-particle distribution function at some initial time 
factorizes in terms of a product of single-particle distribution functions (which is re- 
lated to the assumption of molecular chaos) , and leads to the (irreversible) Boltzmann 
equation [211]. The assumption that the two-particle distribution function factorizes 
is an example of a type of coarse graining called slaving of the two-particle distribution 
function to the single-particle distribution function, in the language of Calzetta and 
Hu [82]. The situation in quantum kinetic field theory is completely analogous. One 
may choose to work with a truncation of the hierarchy of the Wigner function and its 
higher correlation analogues, or one may instead perform a slaving of, for example, 
the Wigner-transformed four-point function, which leads (within the context of per- 
turbation theory) directly to the relativistic Boltzmann equation [68] and the usual 
i?-theorem. Typically this slaving of the higher correlation function(s) involves impos- 
ing causal boundary conditions to obtain a particular solution for the higher correlation 
function(s) in terms of the lower order correlation functions [68,82]. The truncation 
and subsequent slaving of the hierarchy within quantum kinetic field theory can be 
carried out at any desired order, as dictated by the initial conditions and relevant 
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interactions. As with any coarse graining procedure, in implementing the slaving of 
a higher correlation/distribution function to lower correlation/distribution functions, 
one is going over from a closed system to an effectively open system, the hallmarks of 
which are the emergence of dissipation [68, 185] and noise/fluctuations [56,82]. This 
fact has led some to search for stochastic generalizations of the Boltzmann equation 
[212, 213], motivated by that fact that systems in thermal equilibrium always manifest 
fluctuations, as embodied in the fluctuation-dissipation relation [214,215]. 

The essential point about the process of slaving of higher correlation (or distribu- 
tion) functions is that it is a step which is wholly independent of the assumption of 
macroscopic and microscopic time scales. In fact, a completely analogous procedure 
exists at the level of the Schwinger- Dyson equations (i.e., without Wigner transfor- 
mation) for correlation functions in an interacting quantum field theory [82]. Recall 
that the Schwinger-Dyson equations arc, in the context of nonequilibrium field theory 
formulated using the Schwinger-Keldysh closed-time-path, an infinite chain of cou- 
pled dynamical equations for all order correlation functions of the quantum field. The 
importance of the closed-time-path formalism in nonequilibrium situations is that it 
ensures that the equations are causal and that the correlation functions are "in-in" 
expectation values in the appropriate initial quantum state or density matrix. As with 
the BBGKY hierarchy in nonequilibrium statistical mechanics, the general strategy 
is usually to truncate the hierarchy of correlation functions at some finite order. A 
general procedure has been presented for obtaining coupled equations for correlation 
functions at any order / in the correlation hierarchy, which involves a truncation of 
the master effective action at a finite order in the loop expansion [82]. By working 
with an I loop-order truncation of the master effective action, one obtains a closed, 
time-reversal invariant set of coupled equations for the first I + 1 correlation functions, 
(p, G, C3, . . . , C/_|_i. In general, the equation of motion for the highest order corre- 
lation function will be linear, and thus can be formally solved using Green's function 
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methods. The existence of a unique solution depends on supplying causal bound- 
ary conditions. When the resulting solution for the highest correlation function is 
then back-substituted into the evolution equations for the other lower-order correla- 
tion functions, the resulting dynamics is not time-reversal invariant, and generically 
dissipative. As with the slaving of the higher-order Wigner-transformed correlation 
function in quantum kinetic field theory, we have then gone over from a closed sys- 
tem (the truncated equations for correlation functions) to an effectively open system. 
In addition to dissipation, one expects that an effectively open system will manifest 
noise/fluctuations, as shown by Calzctta and Hu for the case of the slaving of the 
four-point function to the two-point function in the symmetry-unbroken field the- 
ory [82]. Thus a framework exists for exploring irreversibility and fluctuations within 
the context of a unitarily evolving quantum field theory, using the truncation and 
slaving of the correlation hierarchy. The effectively open system framework is useful 
for precisely those situations, such as thermalization in the post infiationary Universe, 
where a separation of macroscopic and microscopic time scales (which would permit 
an effective kinetic theory description) does not exist. ^ 

While it is certainly not the only coarse graining scheme which could be applied 
to an interacting quantum field,^ the slaving of higher correlation functions to lower- 
order conx^lalion fuuclions within a particular truncation of the correlation hierarchy, 

^At late times in the thermalization stage, when the quantum field is near equilibrium, an effec- 
tive kinetic description may be justified, but will likely require resummation of hard thermal loops 
[206]. Under such circumstances, even the evaluation of transport coefficients is nontrivial for high 

temperatures [216]. 

^Recently, a paper [206] claimed to study noriequilibrium relaxation without assuming a priori the 
existence of a separation of microscopic time scales and incorporating hard thermal loop resummation. 
However, this study assumed a decohered initial density matrix and near-equilibrium initial conditions, 
and thus did not encompass the most general initial conditions (and neither of the above-stated 
conditions holds true at the end of preheating in realistic inflation scenarios). 
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as a particular coarse graining method, has several important benefits. First, it can be 
implemented in a truly nonperturbative fashion, which is essential for post-inflation 
reheating, where the inflaton variance can be on the order of the tree-level effective 
mass at the end of preheating [99, 125]. This necessitates a nonperturbative resum- 
mation of daisy graphs, which can be incorporated in the truncation/slaving of the 
correlation hierarchy in a natural way. Second, the truncation of the correlation hi- 
erarchy accords with our intuition that the degrees of freedom readily accessible to 
physical measurements are often limited to the mean field and two-point function. For 
example, the transition rate of a particle detector coupled via a m(r)^[a;^(T)] inter- 
action to a quantum field [where m(T) is the detector's monopole moment operator] 
depends on the field's positive-frequency Wightman function [17]. 

Related to the non-existence of a separation of microscopic and macroscopic time 
scales in the conditions which prevail at the end of preheating (where one cannot treat 
the time- varying effective mass in the quasilocal approximation), is the fact that for 
any collection of parametric oscillators with time dependent frequency, the notion of 
vacuum state, and hence, that of particle, is ambiguous [17, 130]. While the growth of 
the variance of the inhomogeneous modes of the infiaton field during preheating can 
be attributed to parametric particle creation [95] , it is unlikely that a well-defined and 
unique particle concept for the inflaton field exists at the end of the preheating pe- 
riod. This is because the inflaton zero mode amplitude and variance contribute a time 
dependent term to the effective mass, which means that modes of the inflaton field 
(or other fields coupled to it) perceive a time-dependent effective frequency. In addi- 
tion, the expansion of the background spacetime also contributes a time-dependence 
to the effective frequency of quantum modes. The closest approximation to a "no- 
particle" state is obtained by defining an adiabatic vacuum, [22] which was discussed 
in Sec. 3.2.3, and in terms of this, one can define an adiabatic particle basis [217]. 
In this study, we perform a coarse graining which is independent of particle repre- 
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sentation, and thereby avoid the subtle issues involved in defining particles in time 
dependent background. 

A third challenge presented by the thermalization problem in post-inflationary re- 
heating is the issue of how to define the entropy of the inflaton field as it proceeds 
towards thermal equilibrium. This is an essential point for post-inflation reheating 
in cosmology. As discussed in Chapter 3, at the end of the slow roll period in infla- 
tionary cosmology, the quantum state for matter fields (other than the inflaton field, 
which has a large expectation value) is to a good approximation given by an adia- 
batic vacuum state [8]. This implies an extremely small entropy density per comoving 
volume element. At the end of the reheating period, the matter fields have reached 
a state of local thermal equilibrium, and the entropy per comoving volume element 
should be given very nearly by finite-temperature field theory calculations. Therefore, 
during the reheating period, the entropy per comoving volume grows by an enormous 
amount.^ An important criterion of the physicality of a particular coarse graining 
scheme is whether it predicts a monotonically increasing entropy during the thermal- 
ization stage. It should be pointed out that defining the entropy of a quantum field 
is an older problem than inflationary cosmology, and dates to early studies of particle 
production and vacuum viscosity in curved spacetime. In these studies, hydrodynamic 
transport coefficients such as bulk and shear viscosity were computed in finite tem- 
perature field theory, and were related to the rate of entropy growth through the first 
law of thermodynamics [33,168,191,218,219]. This approach required the assump- 
tion of near-equilibrium conditions, an "imperfect fluid" (hydrodynamic) form for the 
energy-momentum tensor, and a background temperature T. 

Because entropy is of such fundamental importance, it is useful to discuss how 

entropy can be defined for the various ways (outlined above) of approximating the 

''A rough calculation shows that the entropy per comoving volume element grows by a factor of 
about 10^^° during the reheating period [10] . 
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dynamics of a quantum field. First, for a unitarily evolving quantum field theory 
whose dynamics is a closed system and governed by the quantum Liouville equation, 
it is well known that the von Neumann entropy of the density matrix, 

5vN = -1V[p(i)lnp(t)], (5.1) 

is exactly conserved. If one can assume a separation of macroscopic and microscopic 
time scales, one can go over to the quantum kinetic field theory framework. However, in 
merely carrying out the Wigner transform, one has not sacrificed any information, and 
therefore, one should not expect any increase in entropy. Of course, if one additionally 
makes the assumption of factorization (equivalently, slaving of the Wigner-transformed 
four-point function), one indeed obtains the relativistic Boltzmann equation in the 
binary collision approximation. The Boltzmann entropy Sb defined in terms of the 
phase space distribution f{k, X) for quasiparticles can in this case be shown to satisfy 
a relativistic i7-theorem [68, 220]. However, in the case where there does not exist such 
a separation of time scales, how does one define the entropy of a quantum field? For 
nonperturbative truncations of the dynamics of interacting quantum fields, this is a 
nontrivial question [171]. Intuitively, one expects that any coarse graining which leads 
to an effectively open system with irreversible dynamics will also lead to the growth 
of entropy. This intuition is based on nonequilibrium quantum statistical mechanics, 
where if one has a specific projection operator P projects out the irrelevant degrees 
of freedom from the density operator and retains only the relevant degrees of freedom 
(thus going over to an open system), 

p^{t) = Pp{t), (5.2) 

there exists a formalism for deriving the equation of motion of the reduced density 
matrix p^, and in terms of it, the coarse-grained entropy 

5cG = -Tr[pj,(t)lnp^(i)], (5.3) 
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which win in general not be conserved [208,221]. Another equally powerful method 
adept to field theory is the Feynman- Vernon influence functional formalism [70] which 
has been used to treat open systems [223]. 

Let us give a brief depiction of entropy generated associated with parametric par- 
ticle creation for a free quantum field in an expanding Universe [22], or for an inter- 
acting field such as the A^"^ theory in the Hartree-Fock approximation or the 0(A'") 
field theory at leading order in the large- A/" expansion [117,217], where the dynam- 
ics of the quantum field reduces to a collection of parametric oscillators, each with a 
time-dependent frequency. Since the underlying dynamics is clearly unitary and time- 
reversal invariant in this suitable coarse graining leading to entropy growth 
is not trivially evident. Hu and Pavon [224] first made the observation that a coarse 
graining is implicitly incorporated when one chooses to depict particle numbers in the 
n-particle Fock (or "N") representation or to depict phase coherence in the coherent 
state (or "P") representation. Various proposals for coarse graining the dynamics of 
parametric oscillators have followed [118, 223-230]. The language of squeezed states is 
particularly useful for describing entropy growth due to parametric particle creation 
[222, 223, 226, 231]. For our purposes, the essential features of entropy growth due to 
parametric particle creation which distinguish it from correlational entropy growth to 
be discussed below, are that parametric particle creation involves a choice of repre- 
sentation for the state space of the parametric oscillators, and also usually involves 
an explicit coarse graining which can be expressed in terms of a projection operator 
acting on the density matrix. 

In contrast to entropy growth resulting from parametric particle creation, the 
coarse graining implicit in the slaving of the correlation hierarchy represents a choice 
of relevant correlations versus irrelevant correlations. In this sense, it accords with our 
intuitive notion that only the lower correlation functions are readily accessible to phys- 
ical measurement. However, it is not clear in what sense the nonperturbative slaving of 
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the correlation hierarchy corresponds to the projection of irrelevant variables from the 
full density matrix of the quantum field. ^ Nevertheless, it is clear that once one trun- 
cates the Schwinger-Dyson hierarchy, and carries out a subsequent slaving of a higher 
correlation function to the lower correlation functions, then the resulting effectively 
open system will manifest irreversibility and dissipation. As such, this coarse graining 
scheme should result in entropy growth. Because the equal-time correlation functions 
determine the moment expansion of the (Schrodinger-picture) density matrix, one may 
attempt to associate a reduced density matrix with the dynamical equations for the 
equal-time correlation functions in an effectively open system. Using Eq. 5.3, one 
may then define a coarse-grained entropy for the effectively open system, which we 
call the (Calzetta-Hu) correlation entropy to emphasize its origin in the slaving of the 
correlation hierarchy. This is the entropy which we endeavor to compute in this study. 

It is useful to compare the coarse graining scheme in the Calzetta-Hu correlation 
entropy with the coarse graining scheme proposed by Hu and Kandrup in their study 
of the entropy growth due to particle interactions in quantum field theory [171]. In the 
language of a collection of coupled parametric oscillators, the Hu-Kandrup proposal 
is to define a reduced density matrix by projecting the full density operator onto each 
oscillator's single-oscillator Hilbert space in turn, 

9{k) ^ Tt^^^^p, (5.4) 

and defining the reduced density operator as the tensor product of the projected single- 
oscillator density operators g{k), 

Pn = (S)9ik). (5.5) 

k 

^Recently there has been an interesting attempt to express the slaving of higher correlation func- 
tions in the correlation hierarchy in terms of the Zwanzig projection operator formalism [232], though 
so far it has only been implemented within the framework of perturbation theory. 
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The coarse-grained (Hu-Kandrup) entropy is then just given by Eq. (5.3), from which 
we obtain 

ScG = -Y,Tr[g{k)lng{k)]. (5.6) 

k 

It is interesting to observe that for a spatiahy translation-invariant density matrix 

for a quantum field theory which is Gaussian in the position basis, the Hu-Kandrup 

entropy is just the von Neumann entropy of the full density matrix, because the 

spatially translation- invariant, Gaussian density matrix separates into a product over 

density submatrices for each k oscillator. Like the Calzetta-Hu correlation-hierarchy 

coarse graining scheme, the Hu-Kandrup coarse graining does not choose or depend 

on a particular representation for the single oscillator Hilbert space. In this sense, 

it is not sensitive to parametric particle creation, but instead, it is sensitive to the 

establishment of correlations through the explicit couplings^ between the oscillators 

[171]. The Hu-Kandrup coarse graining also has a direct interpretation in terms of 

projection operator language, which is very convenient from the standpoint of defining 

a coarse-grained entropy. 

In this study, we are interested in the growth of entropy due to the coarse graining 

of the correlation hierarchy by slaving of a higher correlation function. The simplest 

nonperturbative truncation of the Schwinger-Dyson equations for the field theory 

which contains the time-dependent Hartrcc-Fock approximation (necessary for taking 

into account the large variance of the inhomogeneous modes of the inflaton field at 

the end of preheating) is the two-loop truncation of the master effective action, in 

which only the mean field ^, the two-point function G, the three-point function C3 are 

dynamical. All higher order correlation functions obey algebraic constraints, and can 

thus be expressed in terms of the three dynamical correlation functions. While this 

®The Hu-Kandrup proposal can also be applied to both free fields in anisotropic spacetimes, where 
the mode-mode interaction is called an intrinsic interaction, and a self-interacting quantum field, 
where the interaction is an extrinsic interaction [171]. 
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truncation of the Schwinger-Dyson equations is well-defined and could in principle be 
solved, it is disadvantageous for two reasons. First, as stated above, without some 
coarse graining, the system will not manifest irreversibility and will not equilibrate. 
Secondly, and on a more pragmatic level, it is much easier to work with a Gaussian 
density matrix than a non-Gaussian density matrix. Therefore we slave the three-point 
function to the mean field and two-point function, and thus arrive at an effectively 
open system. In principle, a systematic analysis of the coarse-grained dynamics of the 
mean field and two-point function should include stochasticity [82] , but we shall defer 
a study of stochasticity to a future investigation. 

This chapter is organized as follows. In Sec. 5.2, we show how the two-loop trun- 
cation of the correlation hierarchy leads to coupled, time-reversal invariant equations 
for the mean field ^, the two-point function G, and the three-point function C3. We 
then show that slaving of C3 to G and (p leads to an effective open system in which 
the dynamics for G and is irreversible and dissipative. In Sec. 5.3, we show how the 
two-loop truncation of the correlation hierarchy can be reformulated in local equations 
for equal-time correlation functions. In Sec. 5.4, we argue that slaving of the corre- 
lation hierarchy leads to the growth of correlation entropy. In Sec. 5.5 we summarize 
our results and discuss their physical significance, as well as possible extensions of this 
work. 

5.2 Correlation hierarchy and effectively open systems 

We seek a consistent truncation of the Schwinger-Dyson equations for the Minkowski- 
space field theory at third order in the correlation hierarchy, which means that the 
dynamical variables are the mean field (f), the two-point function G, and the three-point 
function C3. As mentioned above, this is the simplest (nonperturbative) truncation 
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of the correlation hierarchy which contains exphcit mode-mode interactions. In this 
truncation of the dynamics, higher correlation functions obey algebraic constraints 
which relate them to the dynamical correlations. Generalizing the CTP-2PI (closed- 
time-path, two-particle-irreducible) effective action to incorporate arbitrary-order non- 
local sources (e.g., three-point sources Jabc{x,y, z), four-point sources JabcdiWjX,!/, z), 
etc.) within the Schwinger-Keldysh framework for nonequilibrium quantum fields, 
one obtains the master effective action, which is a functional which, when properly 
truncated, yields an arbitrary-order truncation of the Schwinger-Dyson equations for 
correlation functions [68, 81, 82, 185]. In general, truncating the master effective action 
at a finite number of loops / yields a truncation of the Schwinger-Dyson equations in 
which only the first l + l correlation functions, ^, G, C3, . . . , C;+i are dynamical, and 
all higher correlation functions are constrained. Therefore the coupled equations for (j), 
G, and C3 which we sock correspond to the two-loop truncation of the master effective 
action. For the A^'^ field theory, the two-loop truncation of the master effective action 
(denoted by r;=2[<^, G, C3] in Chapter 2) has the form^ 

Ti=2[4>, G, Cs] = S[4>] - ^TVlnG + '-A^^'Gba - ^^^^^''GabGcd 

+ ^CABc{G-Y'''{G-Y'''iG-Y'''CA'B'C' - ^cj^^'^^'Cabc^d, (5-7) 

where we have (following [68]) introduced a compact notation using capital letters as 
indices to denote both spacetime and CTP labels [68], i.e., A = (a,x), B = {b,x'), etc. 
In terms of the CTP notation of the previous chapters, 

4>A = 4>ai^), (5.8) 

'^It should also be noted that a systematic study of the three-point function is useful in calculating 
the polarization tensor for the (p^ theory in six dimensions, which is a toy model for the three gluon 
interaction in QCD [233]. 

*In this chapter, we set 7i = 1 for notational simplicity. 
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and 

GAB = Gab{x,x'). (5.9) 
In this new notation, the four-index symbol a^^'-^^ denotes 



a 



ABCD _ abed 



and the two-index symbol a^^ denotes 

As in Chapter 2, A^^ denotes the second functional derivative of the classical action 
with respect to (f), evaluated at (p, 

'^^^ = A^lrl'^] = + ^>^'' - ^<'^''^''''^A'4>B'- (5.12) 
dcpAOcpB 2 

In the two-loop truncation of the master effective action, all n-point correlation func- 
tions for n > 3 are constrained, and can be expressed in terms of (/), G, and C3 [82]. 
The evolution equations for C3, G, and (j) are obtained by functional differentiation of 
the effective action with respect to C3, G, and 4>, respectively: 

liG-Y^iG-y^iG-Y^'CABC - ^a^^^^<^^ = 0, (5.13) 
i^G-y^-'-A^^ + ^a^^^^GAB 

+ lCABc{G-Y''iG-Y^'{G-Y^'{G-y^'CA'B'C' = 0, (5.14) 

BA/i-n_L^2\I .-^EABCl, X 1 A^^^ABEDri 1 '.'^^ABCEr^ n 

(7 [U + m )(t)A + g cr (Pa(Pb(Pc + 2 ^ Gba(Pd + Cabc = 0- 

(5.15) 

With a little bit of algebra, this reduces to 

i^G-y^'CABC - Aa^^^^c^z^G^sGGC = 0, (5.16) 
a^\n + m')4>A + ^cx^^^^ (14>a4>b + Gab) 4>d + ^cx^^^^^C^bc = 0, (5.17) 



Gap + -^o-^^'^^(i>DGFBC = ^^f- 



(5.18) 
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The above three equations are the coupled, time-reversal invariant, dynamical equa- 
tions for the mean field, the two-point function, and the three-point function. This 
truncation of the Schwinger-Dyson hierarchy does not lead to irreversibility or dis- 
sipation. However, let us now slave the three-point function to the mean-field and 
two-point function, which means solving Eq. (5.16) for the history of the three-point 
function with appropriate initial data for C3 [68,82,185]. Formally, the particular 
solution is 

Cabc = -iXa^'^^^'GEAGFBGGC^D, (5.19) 

where boundary conditions have to be imposed to get the homogeneous part of the 
solution. It is the causal boundary conditions introduced to solve Eq. (5.16) which 
introduces an "arrow of time" into the problem. If we view Eq. (5.19) as representing 
an approximation to the dynamics for C3, we may insert it into the two- loop effective 
action (5.7), obtaining 

r[<^, G] = S[4>] - ^IVlnG + '-A^^'Gba - ^(t^^'^^'GabGcd 

+ _LAV^^^^a^'^'^'^'0^<^^,G^A'GBB'Gcc', (5.20) 

which is precisely the two-loop, 2PI effective action for the theory (as discussed in 
Chapter 2). From previous studies [68, 82], it is known that the order term in this 
action leads to non-time-reversal-invariant dynamics for the mean field and two-point 
function. In slaving the three-point function to the mean field and two-point function, 
we have introduced a coarse graining which turns the closed system of (5.16)-(5.18) 
into an effective open system. To see this, let us compute the equation of motion 
for the mean-field (p by taking the functional derivative with respect to 0+ and then 
identifying 0+ = 0_ = 0, as described in Chapter 2. We obtain a real and causal 
equation of the form 

(d + m'^ + ^((^ + ^G^ 4>{x) + Y j d'^x'lm{G++{x,x'f) e{x,x') = 0. (5.21) 
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The theta function 6{x, x') appearing in the integrand dictates that the integral is over 
the past history of 4>; this term is clearly not time-reversal invariant. Furthermore, 
if one substitutes free-field propagators for the Gab in Eq. (5.21), [an approximation 
which corresponds to a quadratic expansion of r[^, G] in powers of (^], it follows that 
the equation for (p is dissipative for momentum modes of above the three-particle 
threshold Qm^ [68,82,234], 

Im(G++(x,xO^)0(x,x') = / ^,e^'-^^-^') H ds ^^^^ (5.22) 

where 

(A; + ie)2 = (fco + ie)2-P, (5.23) 
and the function h{s) is given by [82] 

h{s) = - / dty/{s + m?- tf - 4sm^\ 1 - -—9{s - 9m^). (5.24) 

S 74m2 V t 

Having illustrated how the slaving of higher correlation function (s) to the lower cor- 
relation functions leads to irreversible dynamics, let us momentarily return to the full 
two-loop truncation of the Schwinger-Dyson equations for (p, G, and C3. Although 
Eq. (5.16) readily leads to the time-oriented Ansatz (5.19) for slaving the three-point 
function to (f) and G, it is not convenient if, instead, we wish to simultaneously solve 
the coupled Eqs. (5.16)-(5.18). This is because one must have a closed expression for 
G~^ in order understand (5.16) as a differential equation. Making use of Eqs. (5.13) 
and (5.14), we have 

(G-y^ = z(n + m')a^^ + ':^a^^^^{4>Jn + Gab) 

- '-^a^^'^^iG-Y^hCABC, (5.25) 

from which we obtain a formal expression for 

(M-i)^p, (5.26) 



-1\DE 



(G-1) 



FE 



ABEF 



{4>a4>b + Gab) 



166 



in terms of the inverse of tlie "matrix" Mj^ defined by 

= 5/ + '-^a^'^^^^j.CAEC (5.27) 

Naturally, the inverse {M~^)^q has an infinite power series expansion in A, but for 
our purposes, it is sufficient to work at lowest order in A in the equation for CabCi so 
we will take as a first approximation ~ + 0(A), whereupon the equation for 
Cabc becomes linear, 

iA^'^Cabc = Aa^-^^^c^^GsFGcc, (5.28) 

where A^^ is defined above in Eq. (5.12). From Eqs. (5.17), (5.18), and (5.28), we 
see immediately that this approximation of the two-loop Schwinger-Dyson equations 
results in coupled dynamical equations for (p, G, and C3 which are manifestly time- 
reversal invariant. The essential features of Eq. (5.28) are that it is linear in C3 and 
that, for a spatially translation- invariant quantum state, it does not couple the various 
spatial momentum modes of the spatial Fourier representation of C3 to one another. 

5.3 Dynamics of equal-time correlation functions 

While the previous exposition of the functional integral approach to deriving the (trun- 
cated and coarse-grained) correlation dynamics is useful for a study of the origin of 
dissipation in an effective open system [68,185], it appears to be less convenient for 
computing the entropy of the quantum field, and in particular, making a connection 
between the slaving of C3 and the growth of entropy. This is because in the func- 
tional integral approach, the dynamical equations derived are time-nonlocal. Since in 
this chapter we are concerned with the problem with defining the entropy of an effec- 
tively open system, it would be useful to have a formulation of the dynamics which 
makes reference only to the equal-time expectation values of various products of the 
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Heisenberg field operators $h and which are the quantities which should actu- 
ally appear in a position-basis moment expansion of the Schrodinger-picture density 
operator [118,235]. Let us therefore reconsider the field theory in the Heisen- 
berg picture, and derive the evolution equations for equal-time correlation functions 
and their time derivatives.^ As in the previous section, we will first work with a 
two-loop truncation of the Schwinger-Dyson equations in the same approximation as 
Eqs. (5.28), (5.18), and (5.17). As in Chapter 2, the Heisenberg field operator is $h 
and its conjugate momentum operator, 6h. The Heisenberg equation of motion for 
the field operator has the form 

dH + m2$H + ^$| = 0. (5.29) 

Let us suppose that the quantum state of the system consists of a density operator 
Pjj which is invariant under spatial translations and rotations, and is Hermitian and 
has unit trace, Tr(pjj) = 1. In the Heisenberg picture, the density operator is time 
independent. The expectation value of an observable Oh is given by 

{0)=Tr{pM- (5-30) 

The mean field, defined in Eq. (2.2) above, is then given by 

0(xO) = ($)=TV(ph$h(x)), (5.31) 

which is spatially homogeneous due to the spatial translation and rotation invariance 
of the density matrix and the action for the theory. Following Eq. (2.3), we define 
the fluctuation field (^h by 

^h(x) = $h(x) -0(xO). (5.32) 

The expectation value of </Jh clearly vanishes as a consequence of Eqs. (2.2) and (2.3). 
Inserting Eq. (5.32) into Eq. (5.29), and taking the expectation value, we find 

^m^0+^ (^^ + [Ait)] + 3^[Git)]) = 0, (5.33) 

®This approach was also applied in [236] to the large- limit of the 0{N) model. 
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where we have defined 

= i^uixf) = TripMxf) (5.34) 
= ($h(x)2) = TV(p,$h(x)2), (5.35) 

which are spatially homogeneous because of the translation and rotation invariance of 
the density matrix.^° Making use of Eqs. (5.29) and (5.33), we can derive an equation 
of motion for (fj^, 

(p^ix) + mVH(x) + ^((^1 + 3(^Vh + ^<pI - [A] - 3^[G]) = 0. (5.36) 

Using Eqs. (5.29) and (5.36), we would like to derive coupled equations for the equal- 
time correlation functions and their derivatives. To simplify notation, we will use 
the latin indices i,j,k,l,m,n to denote spatial coordinates, so that Gij{t) stands for 
G{xi,t;xj,t). We will further simplify notation by dropping explicit notation of the 
time t, so that Gij{t) will be abbreviated as Gij. Finally, we drop the H subscript on 
Heisenberg field operators. In this notation, we seek coupled equations for 



Aijk = {ViVjVk)^ (5.37) 

Bijk = '^{'Pi^j^k + V>i'Pj'Pk + 'Pifj'Pk)^ (5.38) 

Cijk = ^{Vi<Pj<Pk + 'PiVj'fk + V'i'Pj'Pk)^ (5.39) 

Dijk = {'fi(Pj'fk), (5.40) 

Gij = i'Pi'Pj), (5.41) 

F^J = \{P^VJ+'Pi'Pj), (5.42) 

Eij = iViPj), (5.43) 



""^Naturally, [G] and [A] are divergent and must be regularized within a consistent rcnornializatiori 
procedure. The divergence structure for [G] for a spatially translation-invariant quantum state is well 
known (see, e.g., Chapter 3 above). We will not discuss the divergence structure of [A]. 
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along with (p and (f>= n. It should be emphasized that in the above definitions a 
limit process is understood which would avoid the appearance of Schwinger terms 
[206,237]. Clearly, Aijk as defined above is just the equal-time limit of C3 from the 
previous section, and Gij is just the equal-time limit of G from the previous section. 
In the two-loop truncation of the correlation hierarchy, the correlation functions C4 
and C5 are constrained [82]. The constraint equations for C4 and C5 are 

i'Pi'Pj'Pk'Pi) = GijGki + GikGji + GiiGjk (5.44) 

and 

{V'iV'jV'k'^lV'm) = GijAklm + GikAjim + GuAjkm + GimAjkl + GjkAum 

+ GjiAikm + GjmAiki + GkiAij^ + GkmAiji + GijnAj^jk- (5.45) 

We now differentiate each of Eqs. (5.37)-(5.43) (as well as tt and 0) with respect to 
time, and apply Eqs. (5.36), (5.44), and (5.45), and we find that within the approxi- 
mation where the equations of motion for spatial Fourier modes of C3 do not couple 
to one another [i.e., where we require agreement with Eq. (5.28)], we obtain a closed 
set of dynamical equations for the equal-time correlation functions and their time 
derivatives, Aijk, B^jk, Cijk, Dijk, Gij, Fij, Eij, (p, and vr, 

Aijk = ^Bijk (5.46) 

Bijk = "^Cijk — Tf'Aijk — —4>{GijGik + GjkGji + GkiGkj) (5-47) 

Cijk = B>ijk — '^SJf'Bijk — —4>{GikFij + GijFik + GijFjk) (5.48) 

bijk = -39Jl^ - \4>{FijFik + FjkFji + FkiFkj) (5.49) 

Gij = 2Fij (5.50) 

Fij = Eij - TfGij - j4>{Aiij + Aijj) (5.51) 

Eij = -2m^Fij (5.52) 

^ TT (5.53) 
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^ - (rn' + ^4>' + ^[G]^ - ^[A], (5.54) 

where Tl^ is the time-dependent Hartree-Fock effective mass defined by 

93t2=rn2 + ^<^' + ^[G]. (5.55) 

The appearance of A^/6 in Eq. (5.54), as opposed to the in Eq. (2.129), is because 
in this Chapter we are working in the time-dependent Hartree-Fock approximation for 
the A$*^ theory instead of the leading-order large- A'^ approximation in the 0{N) theory 
studied in Chapter 2. We remind the reader that [A] = Am and [G] = Ga. It can 
be verified that Eqs. (5.46)-(5.54) are time-reversal invariant, where under the time- 
reversal operator acts as follows 





(5.56) 


Q{Bijk) = -Bijk, 


(5.57) 


&{Cijk) = Cijk, 


(5.58) 


Q{Dijk) = -Dijk, 


(5.59) 




(5.60) 


0(-^ij) — ^ij^ 


(5.61) 




(5.62) 


= 0, 


(5.63) 


e(7r) = -TT. 


(5.64) 



Eqs. (5.46)-(5.54) represent a generalization of the equations derived in [236] to third 
order in the correlation hierarchy, for the case of the model. Because of the time- 
reversal invariance of the dynamical equations, it is expected that the density matrix 
whose moment expansion is given by the equal-time correlation functions (5.37)-(5.43) 

^^In [236], the quartic 0{N) model was studied in the large-AT Umit, which at leading order is 
structurally analogous to the time-dependent Hartree-Fock approximation. 
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would have conserved von Neumann entropy. However, this appears to be difficult to 
prove. 

It is interesting to note that the equations of motion for the time-dependent 
Hartree-Fock approximation are contained in the equations of motion (5.46)-(5.54), 
which is to be expected since this approximation is known to be a subcase of the two- 
loop truncation of the 2PI effective action, as pointed out in Chapter 2. In the present 
context, the time-dependent Hartree-Fock approximation consists of setting Aij^ = 
for all time, and in that case the resulting equations of motion for Gij, Fij and Eij 
take the simple form 

2Fij, (5.65) 
Eij -Tf Gij, (5.66) 
-2m^Fij, (5.67) 

which is seen to agree with the equations of motion in [236] . We note that in the time- 
dependent Hartree-Fock approximation, there is a first integral which can be obtained 
from Eqs. (5.65)-(5.67), 

±{EijGij-F^)=0. (5.68) 
Therefore it will be useful to define 

4 = 4(E,,Q,-4.), (5.69) 

where the aij function is a constant of the motion in the time-dependent Hartree-Fock 
approximation. It is not obvious from the definition (5.69) whether aij is necessarily 
real. However, it was shown in [118] that the spatial Fourier transform of aij, denoted 
by erg, is real and bounded from below, as a consequence of the uncertainty principle, 
that its spatial Fourier transform is indeed real, and bounded from below. We will 
also show in the next section that by including the "setting-sun" diagram which goes 



G 



Eij = 
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beyond the Hartree-Fock approximation (i.e., in the effectively open system discussed 
in Sec. 5.2), the aijk function defined by Eq. (5.69) is no longer constant. 

Let us therefore slave the three-point function to the mean field and two-point 
function as was done in Sec. 5.2, in which case the equal-time correlation function 
Aijk takes the form 

Ajkit) = 2XJ d^ye{t,y'')lm[G++{y,Xi)G++{y,xj)G++iy,Xk)] 4>{y), (5-70) 

where Xj, Xj and x^ are shorthand for {t,Xi), {t,Xj), etc., and the i,j,k label the 
spatial coordinate vectors but are not vector indices. The functions G+_|_(x,x') are the 
time-ordered Green functions satisfying the equation 

{x, x') = —iS{x — x'), (5-71) 

with appropriate "in-in" boundary conditions. The parameter 9Jl^ is the time-dependent 
effective mass defined in Eq. (5.55). In the effectively open system where the three- 
point function is slaved to the mean field and the two-point function, it is evidently 
not possible to represent the coupled equations for 4> and G in terms of completely 
time-local quantities, since the time-ordered propagators appear. This accords with 
intuition gained from the projection operator formalism in nonequilibrium statistical 
mechanics, where the resulting equation for the "relevant" density operator is time 
nonlocal [221]. However, the dynamical equations for cf) and G are still well-defined, 
provided boundary conditions are supplied to get a unique solution to Eq. (5.71). In 
this case, however, the dynamics becomes time-nonlocal and irreversible. 

5.4 Correlation entropy 

In the preceding sections we derived coupled equations for the equal-time correlation 
functions within a two-loop truncation of the correlation hierarchy for the field 

theory, and showed how the slaving of the three-point function to the mean field 
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and the two-point function leads to irreversibility and dissipation. In this section 
we attempt to define the entropy associated with this efi^ectively open system, which 
we call the (Calzetta-Hu) correlation entropy. Recall that the slaving of the three- 
point function leads to coupled equations (5.50), (5.51), (5.52), (5.53), and (5.54) for 
Gij, Fij, Eij, (j), and vr, where Aijk is expressed in terms of the two-point function and 
mean field by Eq. (5.70). Because we have slaved the three-point function, the reduced 
density matrix should have a Gaussian moment expansion in the position basis. In this 
section we investigate the consequences of this fact for the entropy of the effectively 
open system. 

For a calculation of the correlation entropy, we now go over to the Schrodinger 
picture, where and Ilj^ are the spatially Fourier transformed field operator and 
conjugate momentum, respectively. The most general Gaussian density operator sat- 
isfying spatial translation and rotation invariance has the position-basis matrix element 
[118,235] 



(0V|0) = n(2vr@-'/'exp 



i7r((?!>o - ^o) 



2% 



ill 



(4 - ' 



+ 



, (5.72) 



where erg, and ry-^ are all functions of time. The cr-j parameter controls the extent 
to which the density operator represents a mixed state. It has been called the "phase 
mixing" parameter by some authors [235]; we eschew this nomenclature because of 
the alternative meanings of "phase mixing" in statistical mechanics. The equal-time 
correlation functions can be computed directly from Eq. (5.72), 



Tr 



Pirn-, 



£2 



T^[p(*)(%-%o)']='?| + 



k ' 



Tr 



(5.73) 
(5.74) 

(5.75) 
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Then it is straightforward to determine the relations between the time-dependent func- 
tions in the moment expansion of the density operator and the equal-time correlation 
functions, 

Gkit) = ep (5.76) 

m = k^k^ (5-78) 
where, due to spatial translation invariance, the spatially Fourier-transformed two- 
point functions are defined as 

and similarly for Fij and Eij . In the approximation where the three-point function has 
been slaved to the mean field and two-point function, we can use Eqs. (5.50), (5.51), 
and (5.52) to obtain 

so that aij is no longer a constant of the motion. Note that the aj^ defined in Eq. (5.77) 
above is just the spatial Fourier transform of aij, 



''^^'^ = I ^s^'''^''-''^^kii^- (5.81) 



As stated above, due to the uncertainty principle, cr^ is bounded from below [118], 
cr^ > 1- We are now in a position to compute the coarse-grained entropy of the 
reduced density matrix in the effectively open system where the three-point function 
has been slaved, the correlation entropy. Making use of the calculation of [118,235], 
we find 

ScG = -Tr[p{t) In p{t)] (5.82) 



^fe(^) + l^ 1^ ( ^kit) + l \ _ / ag(t)-l \ fa^,{t) - 1 



, (5.83) 
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where cr^ is given by Eqs. (5.69) and (5.81). From Eq. (5.80), we see that S is given 



by 



k 



( 




) 



(5.84) 



where is given by Eqs. (5.80) and (5.81). It is not clear from Eq. (5.84) whether 



5.5 Summary 

In this chapter we have shown how the truncation of the master effective action in the 
Schwinger-Keldysh formahsm leads to coupled, nonperturbative, and causal dynamical 
equations for correlation functions in the field theory, and that the slaving of 
higher correlation functions to the lower correlation functions leads to irreversibility 
and dissipation. Wc then showed that a coupled set of equations can be derived 
for equal-time correlation functions which appears to be equivalent to the two-loop 
truncation of the master effective action. Finally, we showed that the slaving of the 
three-point function to the mean field and two-point function leads to nonconserved 
correlation entropy. The coupled equations derived in Sec. 5.3 may be useful in a 
further study of the thermalization stage in post-inflation reheating. 

There are several directions in which this study might be extended. First, we 
intend to investigate whether the two-loop coupled equations for equal-time correlation 
functions, which are manifestly time-reversal invariant, are also Hamiltonian. This 
would imply a conserved von Neumann entropy for the density operator in this closed 
truncation of the correlation hierarchy. Secondly, we would like to compute the Hu- 
Kandrup entropy for this truncation of the dynamics, and to compare it with the 
correlation entropy computed above in the case of the effectively open system (where 
the three-point function has been slaved to the mean field and the two-point function) . 
Finally, it would be useful to investigate under what conditions the correlation entropy 



ScG is positive definite. 
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is strictly increasing, as this would help to clarify the role that the slaving of the 
correlation hierarchy plays in describing the equilibration of nonequilibrium quantum 
fields. The dynamical equations derived for equal-time correlation functions in Sec. 5.3 
may prove useful in a systematic study of the thermalization stage of a realistic inflation 
scenario. 
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CHAPTER 6 



Conclusion 



In this dissertation, we have studied the nonequihbrium dynamics of quantum fields 
in both curved spacetime and Minkowski space, with particular emphasis on the re- 
heating problem in inflationary cosmology. The primary theoretical tools utilized (in 
various combinations) in this dissertation are the Schwinger-Keldysh closed-time-path 
formalism, the two-particle-irreducible and n-particle irreducible effective actions, and 
the coarse grained effective action. 

We first derived the coupled dynamical equations for the mean field and two-point 
function of a minimally coupled, quartically self-interacting 0(Ar)-invariant quantum 
field in a general curved, classical spacetime including diagrams up to two loops in 
the CTP-2PI effective action. The equations obtained are useful for the study of the 
dynamics of the inflaton fleld during the reheating period of inflationary cosmology, 
and, with a changeover to a tachyonic mass, would be useful for a study of the dynamics 
of a symmetry-breaking phase transition. Various subcases of the two-loop equations 
were discussed, including the leading order large-A'^ expansion, which is of particular 
use in a study of parametric particle creation including back reaction effects. 

Next we studied the dynamics of the 0{N) model in spatially flat FRW spacetime 
at leading order in the large- AT expansion, where the dynamics of the scale factor is de- 
termined self-consistently using the semiclassical Einstein equation. Initial conditions 
appropriate to the end-state of the slow-roll period in chaotic inflation were assumed. 
The coupled dynamical equations for the mean fleld, scale factor, and inhomogeneous 
modes of the inflaton fleld were solved numerically for different values of the ratio 
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of the inflaton mass rn to the Planck mass Mp. Nearly relativistic oscillations were 
assumed, where the initial inflaton amplitude satisfies ~ It was shown 

that for the case where ^ Mp/300, parametric resonance effects are not an efficient 
mechanism of energy transfer from the background field to the inhomogeneous modes 
because of cosmic expansion. This shows that cosmic expansion should be taken into 
account in a study of preheating dynamics in inflationary scenarios, such as chaotic 
inflation, where the inflaton amplitude is Mp/300 at the end of slow roll, and more 
generally, when the time scale for growth of the inflaton variance due to parametric 
resonance is on the order of the Hubble time. 

We then studied fermion particle production in a model consisting of a scalar 
inflaton field coupled via a Yukawa coupling / to a fermion field Fermion 
particle production is expected to be important at late stages during reheating in 
models with unbroken symmetry, after back reaction has caused parametric resonance 
effects to cease. We derived nonperturbative equations for the inflaton mean field and 
two-point function which are dissipative due to fermion particle production. In the 
small-amplitude limit where perturbation theory is valid, we showed that the effective 
dynamics of the inflaton zero mode, at order can be described by a stochastic 
equation. The dissipation and noise kernels were shown to satisfy a zero-temperature 
fluctuation-dissipation relation (FDR) . Furthermore, the normal threshold parts of the 
perturbative coarse-grained effective action at 0(/^) were shown to satisfy an FDR, 
and the noise kernel contributes multiplicatively to the effective stochastic equation 
for the zero mode. The stochastic variance of the inflaton zero mode was computed for 
the late stages of reheating, and it was shown that the rms fluctuations of the infiaton 
amplitude can be on the order of the inflaton amplitude before the end of reheating, 
and under such circumstances, the effect of stochasticity on the zero mode evolution 
should be taken into account. 

Finally, we discussed various proposals for defining the entropy of an interact- 
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ing quantum field, and in particular, the correlation entropy which arises when one 
performs a slaving of a higher correlation function within the context of a system- 
atic truncation of the correlation hierarchy. We showed how the slaving of a higher 
correlation function leads to an effectively open system where dissipation necessar- 
ily arises. We then presented a framework for deriving time-local, coupled equations 
for equal-time correlation functions within the context of a two-loop truncation of 
the Schwinger- Dyson hierarchy (in which the mean field, the two-point function, and 
the three-point function are dynamical). We then computed the correlation entropy 
arising from the slaving of the three-point function to the mean field and two-point 
function. 

One avenue of ongoing research is the application of the above-described methods 
to a study of the dynamics of a nonequilibrium phase transition in a field theory with a 
spontaneously broken symmetry [40]. While phase transitions have long been treated 
using phenomenological methods such as the time-dependent Landau-Ginzberg equa- 
tion [238, 239], a complete, first-principles picture of a nonequilibrium phase transition 
from the viewpoint of correlation dynamics, including domain growth and/or spinodal 
decomposition (the groundwork for which was set forth in [84]), has yet to emerge. 
The master effective action [82] , which allows systematic improvement over mean-field 
and Hartree-Fock calculations, is the preferred tool for this purpose. 

One of the more important cosmological applications of the study of nonequilib- 
rium phase transitions is the problem of computing the density of topological defects 
produced in a GUT-scale phase transition in the early Universe. In particular, GUT 
models with a vacuum manifold which has nontrivial first homotopy group may give 
rise to cosmic strings, which are one of the prevailing candidates for seeding large-scale 
structure [207]. Incorporating stochasticity arising from the slaving of the four-point 
function within the three-loop truncation of the Schwinger-Dyson equations (with un- 
broken symmetry) can in principle be used to compute the variance in the defect den- 
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sity, which would give an indication of the vahdity of most previous calculations [which 
focus on only on the ensemble-averaged defect density [240]] to observational limits 
on primordial fluctuations. In addition, there has recently been a great deal in inter- 
est in (and controversy regarding) the possibility of nonthermal symmetry restoration 
during the preheating dynamics of a quantum field with symmetry-breaking potential 
[97-99, 124, 163,241,242]. The 2PI effective action has already been applied to this 
problem [241]. The CTP-2PI effective action has also recently proven useful in the 
study of the dynamics of disoriented chiral condensates (DCCs) which might be pro- 
duced in relativistic heavy-ion collisions of sufficient energy to locally restore chiral 
symmetry [142-144, 217, 243, 244]. 

The use of the coarse grained effective action, in conjunction with the Schwinger- 
Keldysh closed-time-path formalism, as in Chapter 4 where we studied fermion particle 
production during reheating, may prove useful in several problems in cosmology and 
particle physics. In particular, the coarse grained effective action may be used to 
study the nonperturbative, effective dynamics of soft modes of the gauge field in 
the high temperature, symmetry-unbroken phase of the electroweak theory, which is 
necessary in order to correctly determine the hot electroweak baryon violation rate 
[205,245]. The coarse grained effective action has also been applied to a study of 
relaxation, transport, and thermalization phenomena [190, 206] in scalar field theory 
with a quartic interaction, and should be applicable to a study of effective dynamics 
of soft gluon modes in the quark-gluon plasma, where the microscopic theory is finite- 
temperature QCD. Of most direct importance for infiationary cosmology is the study 
of the coarse grained dynamics of super-horizon modes of the infiaton field during 
the slow roll period, where decoherence and stochasticity are directly related to the 
emergence of a classical picture of primordial density perturbations [53] . Work on this 
problem is ongoing [45, 54]. 

The problem of thermalization in post-infiation reheating is also the subject of 
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continued investigation [184]. As described in Chapter 5, a systematic study of the 
thermahzation stage will necessarily require inclusion of diagrams beyond the Hartree- 
Fock approximation (or equivalently, leading order in the large-iV expansion for an 
0{N) model). At this point, temporal non- localities enter into the effective dynamics 
of the two-point function and mean field, as derived from the 2PI effective action. 
Therefore, the method discussed in Sec. 5.3 of deriving coupled time-local equations 
for equal-time correlation functions should prove useful. More generally, the theoretical 
techniques involved in the study of thermahzation in post-inflationary reheating should 
prove useful for any problem which involves thermahzation of quantum fields in a 
dynamical background, such as GUT phase transitions in the early Universe. 
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APPENDIX A 



Discontinuities of the square diagram 



In this appendix, the seven terms of Eq. (4.83) involving cut propagators are exphc- 
itly evaluated using the Cutosky rules. The second and third terms of Eq. (4.83) 
correspond to normal-threshold singularities in the t and s channels, and are given by 



I 



(27r) 



-Tr 

4 ^'-sp 



F++{q)F+_iq + ki)F__{q + h + k2)F_+{q + fci + ^2 + h) 

= -iDisc[A4(A;i,A;2, A;3)u,=„3=o]^(^2 + ^3) (A-l) 



and 



/ 



(2vr) 



Tr 



F+_(g)F__(g + h)F^+{q + h + k2)F++{q + fci + ^2 + ^3) 

= -zDisc[A4(fcl, k2, fc3)|a2=0; a^+as=l]0{ki + k^), (A.2) 



respectively. The third term in Eq. (4.83) corresponds to the leading-order singular- 
ity of the square diagram [195] (i.e., the solution of the Landau equations in which 
ai,a2,a3 are all nonzero), 

d'^q 



I 



(27r) 



-Tr 

4 -^-^sp 



F+_{q)F_+{q + ki)F+_{q + ki + k2)F_+{q + ki + k2 + k^) 

= zDisc[A4(fel,fc2, A:3)|a,,a.,a3>0]^(fe?)^(-fc2)e(fe3)- (A.3) 



The last four terms in Eq. (4.83) correspond to the four remaining twice-contracted 
singularities, and are given by 

d^q 



(27r) 



-Tr 

4 -^^sp 



F++{q)F+_{q + ki)F_+{q + ki + k2)F++{q + ki+k2 + k^) 

= zDisc[^4(fci, k2, kz)\oi2=a3=Q]0ikl), (A.4) 
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/ 



(2^ 



Tr 



/ 



(27r) 



-Tr 
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(27r) 



Tr 

4 -^-^sp 



+ ki)F+_{q + h + k2)F_+{q + ki + k2 + k^) 

= zDisc[A4(A:i,/C2,fc3)|a3=0; ai+a2=l]'^(^3) > (^-5) 

F+^{q)F^+{q + + ki + k2)F++{q + fci + fcs + ^3) 

= iBisc[A4{ki, k2, k3)\a,=a2=o]diki), (A.6) 

F+_(g)F__(g + ki)F__{q + ki + k2)F_+{q + fei + ^2 + ^3) 

= ?Disc[^4(A;i, ^2, A;3)|„i=0; a^+az=Mkl + k^ + fcg). (A.7) 
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